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CONTROL POWER REQUIREMENTS OF VTOL AIRCRAFT 


PREFACE 


A vertical take-off and landing (VTOL) airplane is inherently difficult 
to handle since it has, among other problems, insufficient control power and 
imposes an extremely high task load on the pilot. Thus an automatic flight 
control system is desirable to improve the handling qualities of .the air- 
craft. A number of trends in guidance and control provide the necessary 
tools for the task. The foremost is probably the application of modern 
control theory to both guidance and control problems. Of equal importance 
is the application of high-speed, general purpose, vehicle-borne digital com- 
puters which make the use of modern control theory possible. The research 
reported here was initiated to investigate certain facets of modern control 
theory and their integration with a digital computer to provide a tractable 
flight control system for a VTOL aircraft. Since the hover mode is the most 
demanding phase in the operation of a VTOL aircraft, the research efforts 
were concentrated in this mode of aircraft operation. 

Subdivided into three sections, this report describes research work on 
three different aspects of the operation of the X-14B VTOL aircraft. In 
the first section, a general theory for optimal, prespecified, closed-loop 
control is developed. The ultimate goal was optimal decoupling of the modes 
of the VTOL aircraft to simplify the pilot's task of handling the aircraft. 
Modern control theory is used in the second section to design deterministic 
state estimators which provide state variables not measured directly, but 
which are needed for state variable feedback control. The third section 
examines the effect of atmospheric turbulence on the X-14B and determines 
a maximum magnitude gust envelope within which the aircraft could operate 
stably with the available control power. 
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Section I is part of a partially completed Doctoral thesis research; 
Sections II and III are Masters thesis research efforts. All sections are 


completely self-contained. Appropriate references and acknowledgements 
have been appended to the sections. Support of this research by the Nation- 
al Aeronautics and Space Administration, Grant NGR-05-004-051 , is gratefully 
acknowledged. 
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I. DESIGN AND APPLICATION OF SPECIFIED 
CLOSED-LOOP OPTIMAL CONTROL 
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DESIGN AND APPLICATION OF SPECIFIED 
CLOSED -LOOP OPTIMAL CONTROL 

NOMENCLATURE 


MATRICES 

A = system matrix 
A = desired system matrix 
= system input matrix 
C^ = system output matrix 
j< = feedback controller matrix 
£ = weighting matrix on the effort 

C[ = weighting matrix on the system states 

II = precompensating matrix 
S_ = Riccati matrix 

VECTORS 

i r = reference input vectors 
u. = input vector 

x = state vector 

~ output vector 

SCALARS 

p - roll rate 

q = pitch rate 

r = yaw rate 

u = forward velocity 

= side velocity 


v 



NOMENCLATURE 


SCALARS 

w = vertical velocity 
L = roll moment 
M = pitch moment 
N = yaw moment 
<f> = roll angle 
9 = pitch angle 
'E = yaw angle 



6 . 


CHAPTER 1 
INTRODUCTION 


In the control design of a vertical take-off and landing (VTOL) aircraft, 
a large number of problems arise. Some of these are general problems associ- 
ated with the complexity of multivariable systems and are related to the dim- 
ensionality of the subject systems and the interaction among their variables. 
The VTOL presents still additional problems to the designer. It has inher- 
ently poor stability and, due to its handling characteristics, requires the 
pilot to perform an unrealistically high number of tasks in order to maintain 
the aircraft at a desired attitude. 

The additional handling difficulties associated with a VTOL arise from 
the supplementary control forces and moments added to the aircraft in order 
to maintain the aircraft in a hover or near-hover condition. For the air- 
plane studied in the project, NASA's X-14B research aircraft, the control 
forces to provide lift and propulsion were obtained by vectoring the thrust 
of its two engines. The control moments for pitch and roll were obtained 
by bleeding air from the engines and ducting it to variable area nozzles 
in the wing tips and the tail of the airplane. 

For the design of an automatic flight control system which would im- 
prove the handling qualities of the X-14B aircraft, it was felt that the 
state variable decoupling techniques would yield the best results. The 
logic behind this decision was that if the pitch, roll, heave, and yaw modes 
were completely decoupled from each other, the pilot's task of flying the 
aircraft would be enormously simpler. This was the motivation for the op- 
timal control technique developed in this section. Although the results ob- 
tained may be used for arbitrary closed-loop system behavior specifications, 
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they are especially well suited for decoupling. 

As background for the control technique developed in this report, a 
few brief comments on decoupling and optimal control theories follow. 

The current literature on the subject on non-interacting, or decoupling, 
control is headed by Morgan [2], who arrived at sufficient conditions for 
decoupling and developed a procedure using a constant gain compensator with 
state variable feedback to obtain a decoupled controlled system. Rekasius 
[3] then obtained some useful results extending earlier work. A more gen- 
eral decoupling problem with static compensation, was solved subsequently 
by Falb and Wolovich [5]. Their result was the first complete solution to 
the decoupling problem for a significant class of linear systems. Gilbert 
[6] added to the results of Falb and Wolovich. The more general state feed- 
back decoupling problem, with dynamic compensation, was formulated and 
solved by Wonhatn and Morse. Additional work on system decoupling has been 
done by Silverman [9], Sato and Lopresti [11], Porter [34], and Yore [35]. 
The last author used a model reference approach to obtain an optimally de- 
coupled system. 

The foregoing methods, while useful in many applications, were only of 
limited use to the X-14B and hence a different approach was sought for the 
aircraft. Since most physical systems cannot be decoupled exactly, a tech- 
nique for approximate decoupling was sought. Thus research efforts were 
turned to the method of approximate decoupling with a specified index of 
performance. 

The theory of optimal control is well defined and conceptually not 
very difficult to apply. However, the methods of optimal control suffer 
from several drawbacks. The major difficulty arises from the fact that, 
even in well defined problems, it is difficult to specify a meaningful 
performance index that is analytically tractable. Furthermore, once a 
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performance index is defined, it is not a simple task to interpret the 
meaning of the elements of the performance index or to manipulate these 
to create desired changes in the dynamic behavior of the controlled sys- 
tem or the structure of the controller. 

In an effort to overcome some of the above difficulties, Solheim [15] 
presented a design of optimal control systems which minimized a quadratic 
performance index for a set of prescribed closed-loop eigenvalues. The 
procedure is as follows. For a performance index of the form 

J = / (x 1 + u] P_ uj dt (1-1) 

o 

The weighting matrix and the eigenvalues of the closed-loop system 
are selected by the designer. Then the 3. matrix, corresponding to the 
prescribed set of eigenvalues, is determined. Finally, the optimal feed- 
back gain matrix is determined from the steady state solution of the Ric- 
cati equation. 

The most exacting level of system design occurs when the engineer is 
given the task of designing a feedback system with specific input-output 
and dynamic behavior characteristics. This implies that a percompensator 
with a particular closed-loop system matrix must be available. Again, 
since precise agreement between specification and physical realization is 
not possible, some compromise result is usually accepted. The research 
reported herein considers an analytical approach based on the premise 
that exact system behavior cannot be dictated. The result is a design 
tool for finding an optimal control policy, with respect to some quadra- 
tic performance index, that yields a closed-loop system approximately 
equal to that desired by the designer. 

The method developed herein is completely general with respect to 
system structure and conceptually is very easy to apply. However, there 



are certain computational difficulties which are discussed in the main text 
of this report. Basically, it enables the control system designer to solve 
the following problem. A system is described by a set of input and output 
matrices and a system matrix A. The specification is a desired matrix A 
and the problem is to find a feedback control matrix K that will yield 
the desired closed-loop matrix and will minimize a performance index based 
on system error and effort. The procedure is to determine K_ by direct 
comparison of A and A , if possible, and then to determine P and Q 
from the result. This approach yields a closed-loop system behavior which 
is similar to the desired behavior. However, the magnitudes of state and 
input variables may be undesirably large. These characteristics may now 
be minimized, since and are known, while the closed-loop behavior 
is maintained at or near the desired condition, by varying the elements of 
P. and £. The final result is a compromise between the approximate equal- 
ity, between the desired and actual closed-loop system matrices and the 
magnitudes of the transients of the state input variables. 

Since A may not be arbitrarily specified, the above method is some- 
what restricted in application to decoupling problems. A more severe 
limitation to the method is that even when A can be calculated, there 
is no guarantee that corresponding weighting matrices ?_ and Q can be 
found. Another limitation of the method is that state variable feed- 
back must be used. Thus some form of state estimation has to be utilized. 

The general theory of specified closed-loop optimal control is devel- 
oped in Chapter 2. In Chapter 3, the results preapplied to controller 
design of the X-14B airplane and a successful design is shown. The 
Appendices include the derivation of the complete dynamic equations for 
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CHAPTER 2 

MATHEMATICAL DEVELOPMENT OF OPTIMAL 
CLOSED-LOOP SPECIFIED CONTROL 

2.1 Statement of the Problem and Assumptions 

In the design of a controller for a mul Invariable system, it is fre- 
quently necessary to impose restrictions on the behavior <. haractcristics 
of the closed-loop system. These restrictions include the location of the 
eigenvalues of the system as well as interaction between the inputs and 
the outputs. The most systematic way of specifying these restrictions is 
to prescribe a particular closed-loop system matrix and a precompensator 
matrix. Ordinarily this is not possible except for certain classes of 
systems. What is normally obtained is a system which resembles the desired 
system. It is not known however, what the degree of resemblance is, and 
what the penalty the designer must pay for the desired operation behavior. 
This chapter presents an analytical procedure that will systematically de- 
sign a controller such that a linear closed-loop system will behave close 
to a prescribed manner and is optimal for some quadratic criterion. 

The assumptions made for the analysis are the following: 

- the plant is represented by a set of linear, time-invariant 

set of state space equations, 

- the plant is completely observable and completely control- 

lable, 

- all states are directly available for measurement or can be 

generated by estimation techniques, 

- the control law is linear and constant, 

- the performance index is a quadratic with constant weighting 


coefficients. 



The class of systems considered consists of deterministic systems 
modeled mathematically by constant coefficient, linear differential equa- 
tions in the form 

i= 4 X + B u (1-2) 

x = c x 

The control policy is 

u[ = - K. x_ + r. (1-3) 

Here, x^ is an n-dimensional state vector, is a p-dimensional control 
vector, y_ is an m-dimensional output vector, and r. is a p-dimensional 
reference input vector. The matrices A, B_, C_, K, and £ are constant 
and have dimensions consistent with _x, and JO II Is a precompen- 

sator matrix which is necessary when additional requirements on input- 
output relations, such as decoupling, are made. 

A desired dynamic behavior may be described by the linear differential 


equations given by 

x = Ax+ ^ r . _ (1-4) 

Substitution of Equation (1-3) into (II -2) yields 

i=(A-BK)x + BJRr . (1-5) 

Thus if a controlled system is to behave according to specifications, it 
must obey the Equations 

A-BK-A (1-6) 

l« = L ■ (1-7) 


Mote that is not affected by the feedback controller K and is 

affected only by the precompensator R. In all cases R. is either 
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known or, for inpul-output decoupled systems, may be computed [4, r >] 
readily. Hence, flic problem reduces to finding an optimal control k 
which minimizes the performance index 

oo 

J - / (x 1 £ X + u l P u)dt . (1-1) 

0 

This is simply the optimal regulator problem except that now £ and £ 
are not specified but must be determined for a given closed-loop condition. 
Matrix £ is symmetric and positive and semidefinite while £ is symme- 
tric and positive definite. 

2.2 Development of the Specified Optimal Control Technique 

The conventional approach to optimal control problems is to specify 
£ and in the performance index and to solve for the optimal control 
policy. It is well known that the optimal control policy for infinite 
time problem is 

u = - £ _1 B/ $ x (1-8) 

where S_ , the Riccati matrix, is the steady state solution of the Riccati 
equation 

S * - S A - A' S-a+llff 1 !; | . (1-9) 

£ is symmetric and positive definite. Comparison of Equations (1-3) and 
(1-8) yields 

K = P" 1 B' S . (1-10) 

Thus, given a set of weighting matrices £ and £ » a corresponding op- 
timal solution is readily obtained. However, there is no guarantee that 
the resulting closed-loop system has the^esS'red dynamic characteristics. 
Instead, one has to evaluate the closed-loop system for several sets of 
values of £ and £ before arriving at a desired solution. 
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To develop the method of optimal specified response, consider a 
system matrix A_ and a desired closed-loop matrix /[ defined by 

A = A - B K (1-11) 

If the inverse of the system input matrix B_ exists, then the control ma- 
trix K may be readily found. For the case when the rank of is less 
than the order of the system, several courses of action exist. One is to 
specify a set of eigenvalues and to compute j( by any technique which 
will provide the desired set of eigenvalues. Finally the determined value 
of K, is used to compute from Equation (1-11). This, of course, re- 

moves some degree of arbitrariness in the specification of A_ , but still 
provides for a compromise between the choice of desired set of eigenvalues 
and the amount of effort necessary to obtain the set. 

A second course of action is to specify A , consistentwith' permis- 
sable values of j( , and to determine £ and £ for the closed-loop sys- 
tem. This is especially attractive when closed-loop behavior like decou- 
pling is desired. 

The mathematical equations necessary to solve the problem are obtained 
by first considering the optimal controller K. , given by Equation (1-10), 
substituting it into Equation (1-11), and solving for A - A. : 

A = A - B P' 1 B‘ S 

A - A = BP' 1 B‘ S. (1-12) 

The last equation, when substituted into the steady state Riccati equation, 
yields 

SA + A 1 S+£-S(A-A)=0 , (1-13) 

which reduces to 

A' S + S A + 0 = 0 


(1-14) 
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When the transpose of Equation ( I ~1 4) is added to it the final result, 
recalling that S_ and £ are symmetric, is 

(A + A) ' £ + S (A + A) = - 2 £ . (1-15) 

Equations (1-10), ( I - 1 1 ) and (1-15) are fundamental to the problem of spe- 
cified optimal control. The method is to determine from (1-11), cal- 
culate £ from (1-10), determining S, simultaneously, and determine £ 
from (1-15). Once P. , which is symmetric and positive definite, and £ , 
which is symmetric and positive semi-definite, are determined, they may be 
modified to maintain a degree of resemblance to the desired closed-loop 
matrix while minimizing state variable and input variables magnitudes. 

The nomdeal systems are analyzed m the next section. 


2.3 Analytical Development for Nonideal Systems 

Consider now the case where the number of manipulated variables is 
less than the number of state variables (p < n), and hence, the matrix 
is nonsquare. Assume that the p inputs are linearly independent from 
each other; consequently, the rank of is p. 

Suppose a change of variables is made with 

x - Gz (1-16) 

where is an n x n nonsingular matrix to be defined later. When this 
is substituted into the system equation 


x = Ax + Bu 
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the result can be written as 

z = G _1 A Gz + G" 1 Bu 


(1-17) 


Def i ne 


A 4 g" 1 AG 


B = G 1 B 


(1-18) 

(1-19) 


Further, define G such that 


B = 


B 


where is an n x p matrix, 0. represents an (n-p) x p matrix and is 
a nonsingular p x p matrix. This may be obtained as follows. 


If G = [G-j G^], then 


B_ - GB ~ [G-] G^] 


0 

B 


or 


and 


( 1 - 20 ) 




Gg = B(B) 


-1 


( 1 - 21 ) 


,-l 


The matrix G-^ can be any n x (n-p) matrix such that ^ exists. This 
demonstrates that a new system representation can be found of the form 


z = Az + Bu 


( 1 - 22 ) 


such that the first (n-p) rows of the matrix B^ are identically zero 
with the remaining p rows forming a matrix of rank p. Therefore, for 
the derivations and analyses which follow, the first (n-p) rows of IB 
can be taken as zero without loss in generality. The lower p rows of 
are assumed to be of rank p. 
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For the development of the control policy for the nonideal system 
where p < n, consider the closed-loop representation 

A = A - BK (I- 


4ll ^12 


-21 ^22 J 


-11 hi 


^21 h .2 ) 


The dimensions of the various matrices are as follows: 

A(n x n), A(n x n), B(n x p), K(p x n), A^ ( j x j), 

A-j ^ ( J x p)» Ag-] ( p x j) 5 Aqq( P x p) > A.'j'j^j ^ J ) j 

A-, 2 (j x p), A 21 (p x j), A^fp x p), 0(j x p) and 

B 2 (p x p), where j = n-p. 

The Riccati equation for optimal control is given by Equation (1-13) 
with the control policy given by 

K = P’ 1 B S . 


[K-, 

i — i 

-11 

hi 

-12 

I 22 

"an 

42 

_ q 12 

S 22 . 
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with K-jtp x j) , i< 2 (p x p), S n (j x j), S 12 (j x p) , S^p x p) , 

in (j x j), 3! 2 (j X p), £22 (p x p) and P (p x p) . 


Working with partitioned matrices. Equation (1-23) can be expressed 


'in 

l 3 * 1 

PO 


in 

A 1 
-12 1 

0 0 

| in 

a 

-22 


fei 

A 1 " 
~ZZ j 

5* hk 


Breaking down this matrix equation into components parts and solving for 
and JO, yields 


£1 = 3 2 


(^21 ~ Bsi ^ 


’21 


(1-24) 


-1 


^2 | 2 (Ag 2 - Ag 2 ) 


Provided JO, ^ exists. This was shown to be true since 10, is of rank p. 
Clearly A.^ - and A^ 2 = A)j 2 are restrictions imposed at this time be- 
cause of the limit in the number of inputs. After the transformation, the 
optimal weighting matrices £ and £ can be found. Repeating Equation 
(1-12) , we have 


BP 1 B S = A - A 


In the partitioned form this becomes 


0 



0 

,12 h t 1 I2 %2 
0 

-22 " -22 
l 12 = A 2 i - A 21 
=22 -22 " -22 


(1-25) 


(1-26) 
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2.4 Application to an Ideal System 

To illustrate the design procedure, consider an ideal second order 
system given by Equation (1-2) where 


A = 

-3 2" 

B = 

"l 

1 

, c = 

1 

-1 

0 

! 

4 - 5 


_0 

1 


_0 

1 


Let the desired closed-loop response be governed by 

-8 0 


A = A - BK = 


0 -8 


To determine and Q for the given A, the first step is to find 


K from 


1 


K = B 1 (A - A) = 


1 -1 

- - - j 4 3 

The next step is to determine a symmetric, positive definite from 

Equation (1-10) keeping in mind the condition that must also be sym- 
metric and positive definite. Thus 


where 


P = 


_1 p. K = 1 

* 


P P 

'll P 12 


c c ! 

a ll 5 1 2 j 


, s = 


[ P 12 P 22 j 

S 12 S 22 _ 

f the last equation 

the result 

1 + 4 P 12 

- 

■ p il + 3 P 12 

1 “ 3 P 12 + 4 P 22 

P 11 " 4 P 1 2 


22 


j S 11 S 12 


S 12 S 22j 


Thus 


S 11 " P 1I + 4 P 12 


(i) 


S 10 - Pn 3 P no 4 Poo “ Pit - 3 P -j ^ (li) 


12 


11 


12 


22 r ll 


S 22 P 11 ' 4 P 12 + 3 P 22 


(ill) 
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Because of symmetry of £, Equation (ii) yields the relationship 

P = ^ P 
*12 1 K 22 ' 

Equations (i) - (ii) become 

5 11 = P 11 + I P 22 

5 12 = “ P ll + 2 P 22 
S 22 = P 11 + 1 P 22 


(iv) 

(v) 

(vi) 

(vii) 


Analysis of the equations suggests that P 22 may be chosen arbitrarily. 

4 

Since £ must be positive definite, P-j-j > P 22 • Thus an acceptable 
choice for P is 


P = 


1 2/3 ' 

2/3 1 


which yields 


S = 


" 11/3 
1 


1 

4/3 _ 


From Equation (1-15), £ is found directly to be 



109 

17 

1 



£ “ 3 




17 

46 


It may be readily verified, using Equation (1-13), that the above values 
of £ and £ will yield the desired value of £ and hence the feedback 
matrix K. 

The time responses for the closed-loop system K ~ A - £ £ are shown 
in Table 1-1 for three sets of initial conditions. While the transient 
behavior is as desired, let us assume that the maximum magnitude of the 
control variable u 2 is excessive. To reduce it, P 22 , the weighting asso- 
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ciated most closely with u 2 is increased from P 22 = 1 to P 22 = 2 while 
all the other elements of P^ and ^ are held constant. The resulting 
feedback matrix, obtained from Equations (I -10) and (1-13), is 


i 2.90 
i 

0.32 

[ 1.56 

1.39 


and the closed-loop system matrix is 


-7.76 


A = 


2.44 


-0.30 ' 
-6.39 


The time responses for this modified closed-loop system are shown in Table 
1-2 for the same three sets of initial conditions used in the original 
system. Note that the desired reduction in the magnitude of u 2 has been 
obtained, but at a loss in the degree of decoupling between the states. 

A system designer would have to choose a compromise between excessive in- 
put and excessive coupling. With the foregoing information, he is in a 
position to make an objective choice of a control policy that most nearly 


meets his needs. 
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TABLE 1-1 

TIME RESPONSE FOR THE SECOND ORDER 
EXAMPLE PROBLEM, P(2,2) = 1.0 


X - j ( O ) 

= o , x 2 ( o ) = 

10 




TIME 

X 1 

X 2 

u i 

U 2 

(V + (u 2 ) 

0.0 

0.00 

10.00 

10.00 

- 30.00 

40.00 

0.3 

0.00 

7.86 

7.86 

- 23.59 

31 .46 

0.6 

0.00 

6.18 

6.18 

- 18.56 

24.75 

0.9 

0.00 

4.86 

4.86 

- 14.60 

19.47 

1.2 

0.00 

3.82 

3.82 

- 11.48 

15.31 

1.5 

0.00 

3.01 

3.01 

- 9.03 

12.04 

1.8 

0.00 

2.36 

2.36 

- 7.10 

9.47 

2.1 

0.00 

1.86 

1.86 

- 5.59 

7.45 

2.4 

0.00 

1 .46 

1 .46 

- 4.39 

5.86 

2.7 

0.00 

1.15 

1.15 

- 3.46 

4.61 

3.0 

0.00 

0.90 

0.90 

- 2.72 

3.62 

X ^ O ) 

= 10 , x 2 ( 0 ) 

= 0 




TIME 

X 1 

X 2 

U 1 

U 2 

jdV + <u 2 ; 

0.0 

10.00 

0.00 

- 10.00 

- 40.00 

50.00 

0.3 

7.86 

0.00 

- 7.86 

- 31 .46 

39.33 

0.6 

6.18 

0.00 

- 6.18 

- 24.75 

30.93 

0.9 

4.86 

0.00 

- 4.86 

- 19.47 

24.33 

1.2 

3.82 

0.00 

- 3.82 

- 15.31 

19.14 

1.5 

3.01 

0.00 

- 3.01 

- 12.04 

15.06 

1.8 

2.36 

0.00 

- 2.36 

- 9.47 

11 .84 

2.1 

1.86 

0.00 

- 1.86 

- 7.45 

9.31 

2.4 

1.46 

0.00 

- 1.46 

- 5.86 

7.33 

2.7 

1 .15 

0.00 

- 1.15 

- 4.61 

5.76 

3.0 

0.90 

0.00 

- 0.90 

- 3.62 

4.53 

x ^ o ) 

= 10 , x 2 ( o ) 

= 10 




TIME 

X 1 

X 2 

U 1 

U 2 

(Ut ) + (U 2 

0.0 

10.00 

10.00 

0.00 

- 70.00 

70.00 

0.3 

7.86 

7.86 

0.00 

- 55.06 

55.06 

0.6 

6.18 

6.18 

0.00 

- 43.31 

43.31 

0.9 

4.86 

4.36 

0.00 

- 34.07 

34.07 

1 .2 

3.82 

3.82 

0.00 

- 26.80 

26.80 

1.5 

3.01 

3.01 

0.00 

- 21.08 

21 .08 

1.8 

2.36 ' 

2.36 

0.00 

- 16.58 

16.58 

2.1 

1 .86 

1 .86 

0.00 

- 13.04 

13.04 

2.4 

1.46 

1.46 

0.00 

- 10.26 

10.26 

2.7 

1.15 

1.15 

0.00 

- 8.07 

8.07 

3.0 

0.90 

0.90 

0.00 

- 6.35 

6.35 
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TABLE 1-2 

TIME RESPONSE FOR THE SECOND ORDER 
EXAMPLE PROBLEM, P(2,2) =2.0 


1. X ] (0) = 0, x 2 (o) = 10 


TIME 

X 1 

X 2 

u i 

U 2 

(u.,) + (u 2 ) 

0.0 

0.00 

10.00 

-4.88 

-13.35 

18.23 

0.3 

0.04 

8.27 

-4.16 

-11 .11 

15.28 

0.6 

0.06 

6.84 

-3.55 

- 9.26 

12.81 

0.9 

0.08 

5.66 

-3.01 

- 7.71 

10.72 

1.2 

0.09 

4.68 

-2.56 

- 6.42 

8.98 

1.5 

0.09 

3.88 

-2.17 

- 5.34 

7.51 

1.8 

0.09 

3.21 

-1 .83 

- 4.45 

6.29 

2.1 

0.08 

2.66 

-1 .55 

- 3.70 

5.26 

2.4 

0.07 

2.20 

-1.31 

- 3.08 

4.40 

2.7 

0.07 

1.83 

-1 .10 

- 2.57 

3.68 

3.0 

0.06 - 

1.51 

-0.93 

- 2.14 

3.07 

X^O) 

it 

2 

o 

X 

ro 

o 

= 0 




TIME 

X 1 

X 2 

U 1 

U 2 

UV + (u 2 ) 

0.0 

10.00 

0.00 

-29.84 

-17.80 

47.65 

0.3 

7.92 

0.53 

-23.91 

-14.82 

38.73 

0.6 

6.28 

0.87 

-19.17 

-12.34 

31 .51 

0.9 

4.98 

1.06 

-15.38 

-10.28 

25.66 

1.2 

3.95 

1.14 

-12.35 

- 8.56 

20.91 

1.5 

3.13 

1 .16 

- 9.92 

- 7.13 

17.05 

1.8 

2.48 

1 .12 

- 7.98 

- 5.93 

13.91 

2.1 

1.97 

1 .06 

- 6.42 

- 4.94 

11.36 

2.4 

1.57 

0.99 

- 5.17 

- 4.11 

9.29 

2.7 

1.24 

0.90 

- 4.16 

- 3.42 

7.59 

3.0 

0.99 

0.81 

- 3.36 

- 2.85 

6.21 

X^O) 

- 10 , x 2 (o) 

= 10 




TIME 

X 2 

X 2 

U 1 

U 2 

(u 7 ) + (u 2 : 

0.0 

10.00 

10.00 

-34.73 

-31.15 

65.89 

0.3 

7.96 

8.81 

-28.08 

-25.94 

54:02 

0.6 

6.34 

7.71 

-22.72 

-21 .60 

44.32 

0.9 

5.06 

6,72 

-18.40 

-17.99 

36.39 

1.2 

4.04 

5.83 

-14.91 

-14.98 

29.89 

1 .5 

3.22 

5.04 

-12.09 

-12.47 

24.57 

1.8 

2.57 

4.34 

- 9.81 

-10.39 

20.20 

2.1 

2.06 

3.73 

- 7.97 

- 8.65 

16.63 

2.4 

1.64 

3.19 

- 6.48 

- 7.20 

13.69 

2.7 

1 .32 

2.73 

- 5.27 

- 6.00 

11.27 

3.0 

1.05 

2.33 

- 4.29 

- 4.99 

9.29 
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CHAPTER 3 

LIMITATIONS IN APPLICATION OF SPECIFIED CLOSED-LOOP 
OPTIMAL CONTROL TO X-14B AIRCRAFT 

3.1 Control Design Problems in VTOL Aircraft 

The major problems in the control of VTOL aircraft, occur at the hover 
and near-hover conditions. Because of the low speed of the aircraft at 
these conditions, insufficient aerodynamic forces exist to produce an 
effective control system relying on aerodynamic control surfaces. Hence, 
supplementary controls are added to provide lift, roll moment, and pitch 
moment. With the added controls the pilot is required to perform a large 
number of tasks to maintain the aircraft at a desired attitude. 

To simplify the task of controlling the airplane, a decoupling con- 
troller scheme was sought. The specified closed-loop optimal control meth- 
od, developed in the previous chapter, was selected. The primary reason 
was that while linear decoupling techniques gave adequate results and these 
were used to obtain a completely decoupled system, it was necessary to 
minimize the magnitudes of some of the control and state variables, and 
it was necessary to maximize the envelope control system operation in 
the linear domain. The latter implied avoiding saturation of the control 
variables as much as possible. It is obvious that the handling character- 
istics of a controlled system are drastically changed when saturation is 
reached. To maximize the envelope of operation, it is necessary to dis- 
tribute the control effort in an equitable manner. 

3.2 Description of the X-14 VTOL Aircraft 

The Bell X-14B VTOL aircraft is a fixed wing, jet-propelled, deflected- 
jet vehicle. It is a small (4200 pounds) two-place, side-by-side airplane 
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with two jet engines mounted in the fuselage. The vertical thrust is 
derived from cascade-type diverters on the exhausts of each engine which 
enables the pilot to select any desired thrust angle from vertical to hori- 
zontal. In addition to this thrust diverter, it is also assumed that a 
side vane is present in the exhaust enabling direct lateral movement. 
Angular control about the three axes is maintained through a continuous 
bleed from the engine compressors and discharged at the tail and wing tips 
of the airplane. 

3.3 Equations of Motion for a VTOL Aircraft 

The complete development of the nonlinear equations of motion for a 
VTOL aircraft in general and the X-14B in particular is made in Appendix 
A-I. A general method to linearize an n tr order nonlinear equation is 
presented in Appendix A-I I - Appendix A- III is a listing of the linearized 
system equations for several flight conditions of the X-14B. 

A common practice in analyzing aircraft is to separate, or decouple, 
the longitudinal mode of motion from the lateral mode. The linear equa- 
tions of motion for the longitudinal mode of the X-14B in hover, after 
dropping the effects of the elevator, can be expressed as 

x = Ax + Bu 

where 

0 1.0 0 0 

0 -.150 -4.4 x 10~ 4 -1.5 x 10 

4 = -32.0 0 - .020 - .038 

000 - .021 
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where 


x = 


0 

.035 


0 

.020 


0 

0 


B_ 

= ] 

0 

.56 


0 



V. 

0 

0 

3.6 ^ 


C 


r 





X 1 


e 


~ 5 y ' 


~ u i 

X 2 

= 

q 

> u = 

a 

= 

u 2 

X 3 


w 











_ U 3 _ 


The physical variables corresponding to the state variables are 8 (radians), 
q (radians/second), u (feet/second), and w ( feet/ second ) , respectively, 
with the control vector corresponding to 6 (reaction nozzle about the y 

V 

axis), o (thrust angle measured positive from the vertical), and n (engine 
RPM as a percent of maximum) . 


3.4 Design of a Controller for Hover Flight 

In the open-loop system there is strong coupling between roll and the 
vertical velocity. Let the feedback control policy be given by 


r 





k n 

k 12 

k 13 

k 14 

i 


u_ = - 

Kx = 

I 

K 21 

k 22 

k 23 

k 24 

X 



j 

k 31 

k 32 

k 33 

k ^ 
K 34 J 


To obtain 

the desired decoupling a 

closed-loop sys'tem matrix may be 

specified 

as 

~ 0 

1 

0 

0 " 





-4 

-3 

0 

0 




A » 

0 

0 

-.6 

0 


(1-27) 



L o 

0 

0 

-.6 . 
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It must be determined that the system is completely controllable and com- 
pletely observable. Since state variable output is assumed, observability 
is complete. It is controllable since the rank of the controllability 
matri x 

[B AB A 2 B A 3 B] 

is four. 


Next, it should be determined that a feedback control matrix K, 
exists that will produce the closed-loop response, A. The general form 
of this nonideal case was treated in the previous chapter. The closed- 


loop equation is 

fo 1 


0.-035k n 

-.15-.035k-|2 

- 02k 2l 

-.02k 22 

32 . — . 56k 2 1 

0- . 56k22 

0-3.6k 31 

0-3.6k 32 



(1-28) 

0 

0 

-1 .5xl0~ 4 -.035k 13 

-1 .5xl0~ 4 -.035k 

■* 02k 23 

-.02k 

- .02- . 56k 

- .038- . 56 k 

°: 3 - 6k 33 

- . 021-3. 6k^j 


Comparing Equation (1-27) with Equation (1-28), the control matrix is found 
to be r - 


146.94 

81.43 

-0.60 

0.03 

-57.14 

0 . 

1.03 

-0.06 

0 . 

0 . 

0 . 

0.16 


(1-29) 


The next step is to determine P and Q which, when used in conjunction 
with optimal control theory, would yield the same matrix K. Following 
the same procedure as given in the example of Chapter 2 a set of positive 
semi- definite matrices P and Q , which yields Equation (I- 10) as the 
optimal solution was found. This set one solution to inverse quadratic 
optimization problem for the X-14B airplane is 
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50.000 

28.416 

-.32301 



28.416 

17.150 

- 244.47 

5 

(1-30) 

-.32301 

-244.47 

60676. 

"I 


649794. 

0 

-20.918 

71.434 


0 

39392. 

116.83 

-3872.8 

(1-31) 

-20.918 

116.83 

1.1395 

-43.576 


_ 71.434 

-3872.8 

-43.576 

1683.6 . 



The corresponding solution to the steady state 


Riccati equation is 


given by 


489834. 

163518. 

-133.74 

3867.3 

163518. 

116327. 

-22.548 

-7.3061 

-133.74 

-22.548 

1.8539 

-70.279 

2867.3 

-7.3061 

-70.279 

2715.3 


(1-32) 


As a check, the above values of and Q were used in the per- 
formance index (1-1) and the resultant optimal solution, obtained by 
conventional optimal control theory techniques, was found to be 


A = 


0.00 

1.00 

0.00 0.00 

-4.00 

-3.00 

4.77xl0“ 4 -2.17xl0~ 9 

-2.35xl0" 6 

-1.61X10" 6 

-0.60 -l.llxlO" 5 

-6.08xl0' 8 

-4.17xl0” 8 

7.31xl0“ 9 -0.60 


Note that the above matrix is nearly identical to the desired matrix 
given by (1-29) and, hence, the first part of the controller design 
problem is complete. 
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02 
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OA 
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a 0 
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28.416367 
-0. 37300667 

Oft/M = 50.000000 


28,416362 

»* • #« 

-244.46929 
P62.73 « 12.1497930 


•0 <32300667 
-744.46929 


P t 3 * 3 ) » 60675.5120 


>46,939 
-57,]4’9 
1 ,5d«4Ur*Ofl 


41 ,6?A6 "0 , 604 40 ft 

7,8 754*72-06 1 *03571 

1 . 1575662-00 “2 . 0 31 6 86E-09 


7(1,13 * "50.000000 


3. 4478542*07 
-6.783733E-02 
0 .160033 


,0 1.00000 >0 .0 

-6.00000 "3.00000 4.7756482-11 •? . 1 7 1 76 82-09 

-2.3507882*06 -l , 6 102502 *06 -0.600000 -l . 10V674E-05 

-6.0833912*08 -4 . 1 6 7 1 6 6 2-08 7.31 3349E*09 -0.600000 
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P ( 3 . 3 1 * 61232.9670 
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1 ,(16<1632 r - o » 


41,4766 -0.604400 

1,1831542-06 1.03571 

7,514=412-09 -2.1036792-1? 


PC1.1 ) * 50.000000 
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-fi. 65068 
0.125699 


-4.00000 
-1 .500635E-06 
-3.8470752*03 


1.00000 

-3.00000 

-1,0545662-06 

-2.7053432-08 


P ( 2,2) * 17,1497930 


P ( 3, 33 = 60063.7570 


• 0 

5.1189572-13 

-0.600000 

7.5737332-12 


.0 

9.395772E-04 
4 ,80638 
*0.473516 


146.939 

•57.1*78 

4*056542^-06 


*1*4286 -0.604408 

6 , 44 1 J?0E-O6 1.03571 

7, 6174002-08 -5.B99339P-17 


P(i,( 3 » 50.000000 
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,0 1.00000 .0 *0 
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-1.4603552-07 -9.4226402-08 2,1737622*11 *0*956442 
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•2.4867372*02 
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.0 
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28 b. 
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To determine whether the above controller was acceptable, a linear 
model and the controller were simulated on the digital computer. Further- 
more, the diagonal elements of were varied by 1 percent sequentially 
and the corresponding control matrices were determined. These are shown 
in Table 1-3. The resulting closed-loop systems were observed for trans- 
ients due to initial conditions and the results are shown in Table 1-4. 

The parameters recorded were the maximum values of the state and the in- 
put variables. Examination of the data for the exact decoupling control- 
ler reveals that the control effort u-j (reaction bleed control, Sy) is 
excessive and should be reduced. The penalty on this variable, , 
when increased by 1 percent causes a near 90 percent decrease in the 
maximum value of u-j . This was the evaluated to be the most desirable 
of the controllers presented in Table 1-3. With this choice of p^-j , 
the final closed-loop configuration selected is 


0 

1.00 

0 

0 

-4.023 

-3.013 

9.701xl0“ 5 

-5.536X10" 6 

-126.6 

-67.40 

1 

— i 

CO 

, — J 

-2.408xl0~ 2 

-3.277 

-1.744 

1.092xl0~ 2 

-0.6005 


The disadvantage of this design is that there is some unilateral 
coupling. That is, pitch and pitch rate (0 and q) disturb the for- 
ward and vertical velocities (u and wj , but the last two affect 
0 and q only minimally. Thus, the "best" controller selected is 

r i 



18.416 

13.040 

-0.176 

0.010 

K = 

168.949 

120.348 

0.282 

-0.025 


0.911 

0.484 

-0.003 

0.161 
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3.5 Conclusions 

One of the goals of this study was the development of a suboptima 1 
decoupling technique for general systems. Even though the theory of op- 
timal control is well defined, it suffers from the fact that it is dif- 
ficult to prescribe a meaningful performance index. When a degree of 
decoupling is desired in a controlled system, there is no means of 
specifying a performance index which would provide a decoupling control- 
ler. Thus the problem was investigated in a reverse order. A decoupl- 
ing control policy was determined and then an attempt was made to find 
a performance index which, when used in a standard fashion, yielded 
the desired closed- loop system. The objective was to manipulate the 
known performance index to strike a balance between the degree of tol- 
erable coupling and some criteria on magnitudes of states and inputs 
of the system. 

The elements of the performance index were perturbed a finite 
amount, one element at a time. The resulting control policies were 
evaluated and the most desirable policy was finally selected. There 
was an increase in coupling with the new policy, but that was off set 
by the large decrease in the amount of control needed to maintain the 
airplane in the proper attitude. 

The technique developed in the foregoing appears to be quite pro- 
mising. Additional work, however, is necessary to define the most de- 
sirable type of control action. 
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APPENDIX A 

DEVELOPMENT OF THE EQUATIONS OF MOTION FOR THE X-14 AIRCRAFT 

A. General Equations of Motion for VTOL Aircraft 

This development of the equations of motion for a VTOL aircraft as- 
sumes the aircraft is rigid and the origin of the coordinate system is at 
the mass center [1, 2, 3J. Additional assumptions used in this develop- 
ment of the equations of motion are: 

1. The VTOL aircraft has diverters and vanes affecting the 
engine exhaust. 

2. Any wind disturbance is irrotational . 

3. Engines rotate in the same direction at the same speed. 

A standard aircraft body axis system is used where the aircraft is 
symmetric about the x-z plane with the positive y-axis pointing out the 
right wing, the z-axis down, and the x-axis. in the forward direction of 
flight (stability axes). This is a right-handed coordinate system. 

For the system, the velocities, angles, angular velocities, forces 
and moments associated with each axis are shown in Table A-I-l and Figure 
A-I-l [4]. 




Table A- 

'1-1 



Axis 

Linear Velocity 

Force 

Angle 

Angular Velocity 

Moment 

X 

u 

F x 

<f> 

P 

L 

Y 

V 

F y 

e 

q 

M 

Z 

w 

F z 

y 

r 

N 
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From Etkin [1], the time derivative of the angular quantities can be 
expressed as 

6 = q cos ij> - r sin <j) 

<p = p + q sm <j> tan 0 + r cos <j> tan 0 
^ = (q sin (p + r cos cj>) / cos 0 

Since the analysis is for a VTOL aircraft, it becomes necessary to in- 
troduce some notation not normally found in airplane equations. A listing 
of this additional notation is in Figure A-I-2. 

Inertial Forces 

In writing the equations of motion, all inertial components are ex- 
pressed as D'Alembert forces, i.e., 

EF . . + EF . . . , = 0 

external inertial 

where 

^inertial ~ ~ ma 

Using this notation the inertial forces along the three axes are 
F . = -m (u + wq - vr) 

A I 

Fyi = -m (v - wp + ur) 

F . = im (w - uq + vp) 

Likewise, the moments are 

M xi = - P T x + ' I z '> + + Pfl) l xz 

My, = - q I y + Pr (1 Z ‘ r x> + (r2 - P 2) ‘xz 

«zi = - r I z + pq ll x - I y ) + (p - qr> I xz 
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1 . 


2 . 


3. 


4. 


5. 


6 . 


7. 


8 . 


9. 


ru. 


11 . 


12 . 


13. 


14. 


15. 


16. 


17. 


18. 


19 . 


Figure A-i-2 

Definitions of Non-standard Notation 
w subscript is for wind velocity 
V_j = jet exhaust velocity from the diverters 

0 = thrust diverter angle measured in the verticle plane from the 
positive z-axis 

1 «= exhaust side vane angle measured positive about the x-axis by 
the right hand rule 

U 0 “ [ < u ~ V 2 + (W ~ W w )2]1/2 

U Q = [(u - u w ) 2 + <w - w w ) 2 + <v - v w ) 2 ] 1/2 
a ** angle of attack; sin a = (w - w^)/Uq 
8 «= side slip angle; sin 8 = (v - v^J/Uq 

x^ ** distance engine intake i3 in front of the center of gravity 
z ^ => distance engine intake is below the center of gravity 
x^ - distance engine exhaust pivot point is in front of the center of 
gravity 

Z£ «= distance engine exhaust pivot point is below the center of gravity 
2^ is the effective length where the exhaust will impinge on the diverter 
vane from the end of the jet engine 

I “ moment of inertia of one engine about the axi3 of rotation 

*» angular velocity of an engine measured positive in the direction 
of the x-axis 
N.E. =* number of engines 
T ** total thrust from the diverters 

L , M , N apparent aerodynamic damping near hover not directly 
p q r 

attributable to the effects caused by C. , C , C 

<■ m n 

p q r 

Lg , , N Moment effects due to reaction control nozzles 

x y ^ z 
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20 . 


21 . 


22 . 


23 . 


V *=jet velocity at the engine exhaust 
-*o 

T » total installed jet thrust 

*= longitudinal dLverter efficiency factor 
= lateral diverter vane efficiency factor 
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where 


I “ /(y 2 + z 2 )dm, 

X 


xz 


/xzdm 


I = /(x 2 + z 2 )dm, 

I = /xydra 

xy ' 


T A 2 , 2 X 
I « /(x + y )dm 

z 

I » fyzdm 
yz 


with 


I “0 and I 
xy yz 


0 


since the aircraft is symmetric about the xz plane. 

In addition to aircraft inertial terms, there is a contribution from 
the rotation of the engines to the moment equations. If the engine speed 
is almost constant, the contribution to the moment equations can be expressed 
as 

K * 0 
e 

x 

M * (N.E.) I « r 
e e e 

y 

M « -(N.E.) I n q 
e e e 

z 

Gravitational Forces 

Since the origin of the coordinate system is at the center of gravity, 
the gravitational forces do not enter into the moment equations. Their 
contribution to the force equations is 


-W sin 0 


g 


F =* W sin <{> cos 0 

y 

g 

F = W cos tj> cos 6 
z 

% 

Mass Flow Effects 

There are two contributions to mass flow effects, one at the engine 

i » 

intakes and a second effect at the engine exhausts. If the mass of the fuel' 
burned is neglected, the forces due to the mass flow can be expressed as 


t 

X. 

m 

— ta 

(u - u - 

w 

v 3 

sin <f cos 



* ~m 

(v - V - 
w 

V J 

sin a) 


» 

z. 

m 

a — m 

(w - w + 
w 

V J 

COS 0 COS 

A) 
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where ra is the entering mass flow rate which is approx tmatcJy equal to the 
mass flow rate out, is the relative exhaust 1 velocity of the jet, ° is 
the diverter angle, and X is the side vane angle. Recalling that the gross 
thrust equals mass flow rate times relative velocity, 

T « m 
or 

T 

m “ v 

j 

Upon making this substitution, the mass flow contributions to the equations 
of motion can be written as 

F *= T [sin a cos X - (u - u )/V ] 
x. w j 

m 

F - T [sin X - (v - v )/ V. ] 
y » w 1 

•'m J 

F =-T [cos o cos X + (w - w )/V.] 
z. w 1 

m J 

The mass flow contributions to the moment, equations enter in because the 
intakes and exhausts are not at the center of gravity. The moments are 


M 


x. 

m 


in 


M 


m. 


T 2 i (y ' v „ )/v j 

T (w - w w ) - Z 1 (u - U b )]/V 3 


xn 


M 


z. = ' T X 1 <v - v w >/? j 
"in 


where and are the distances from the center of gravity to the center 
of the intakes. The moments caused by the exhaust are 


M 


m 


out 


M 


m 


out 


M 


z. 

m 


^-T sin ^( Z 2 cos °) 


T cos A (z^ sin o + cos «*) 


T sin A (x^ — sin a) 


out 


where x ^ and are the distances from the center of gravity to the center 
of the exhausts, and is the effective length where the exhaust will 
impinge on the diverter vane. 
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The total moments caused by engine mas3 flow terms can now be 
expressed as 

M =■ T [-sin X (z 0 + Z cos a) + .z (v - v )/V, ] 

z 1 1 w J 

M a T { cos X ( z 2 sin a +■ x^ cos a) + [x^ (w - w^) - z^ (u-u w )J/V.] 
a 3 

M ** T [sin X(x,. - sin a) -■ x, (v - v )/V_] 
z » z i 1 w i 

m J 

Additional moments are caused .by discharge nozzles in the tail and on the 
wing tips. These reaction control moments can be written as 
M 
M 


*= Lr 5 
*R.C. 5 x x 

=•• M. <5 
y R.C. °y y 

M « N,. 6 

Z R.C. 6 z Z 

where <5 , 6 and 6 are the reaction nozzle openings expressed in degrees, 
a y z 

Aerodynamic Expressions 

me complete aerodynamic expressions for a convene onal airplane can 
be written as 

x » S-u 2 s c 

aero z o x 

Y * U 2 S C 
aero 2 o y 


Z =|D^SC 
aero 2 o z 


2 

£u _ 2 S b C, 


aero 2 o 

M = ^ U 2 S c C 

aero 2 o m 

N »|? 2 SbC 

aero 2 o n 

where C , C , C , Co , C and C are the non-dimensional aerodynamic 
x y z * m n 

2 * 2 1/2 

coefficients, o is the air density, U q ** [(u - u ) + (w - w^) ] , 

U » [(u - u ) 2 + (w - w ) 2 + (v - v ) 2 ] li ^ 2 , S is the wing area, b is the 
o w w w 


span and c is the wing chord. 
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Now use a Taylor series expansion of the non-dimensional aeroaynnmic 

forces retaining only the first terms In the series [1, 2] and let the 

A 

vehicle velocities relative to the air mass be denoted by u ^ = u - u^, 

A A 

v tr, v - v and w , => w - w . Then, 
rel w rel w 

Sc 9 C 9C 

C = C + • 5 -£cz+ -wAi , + -r~ £ 
x x„ da rel So e 

0 e 

SC SC SC SC 

C = C + 8 + P + "5“^" r + 5 

y y^ 38 Sp Sr So^ r 

3C SC SC SC SC 

C = C + - a + -r— ^ u . + -r-^- q + a + ~ <5 

z 3a Su rel 3q 1 Sa So e 

0 e 

sc sc sc sc sc A 

C I ■ C £_ + UT 6 + w P + 3T r + TS- S 0 + IT 6 r 

0 r a r 

3C SC SC SC SC 

~ . iu. in , m ni ♦. , ni ^ 

C = C + -T — a d -5 — u . + -5 — q + -57— ■ a + r-r— ° 

m iDq Sa Su rel Sq 4 Sa S6 e e 

SC SC SC SC SC 

C n - C „„ + IS* 6 + 3T » + ?? 1 + IT 5 a + Tp 6 r 
0 a r 

where a is the angle of attack (sin a = W rel/ ^0^ * ^ C ^ e S ^P 

angle (cos 8 = U / U~ ) , and 6 , 6 and 6 are the elevator, aileron and 
U ' U 6 3. r 

rudder deflections, respectively. 

Following Etkin [2], if the partial derivatives are written as 
non-dimensional stability derivatives, the aerodynamic portions of the 
equations of motion become 


X =|U Z S[C +C a + 7 : C u + C 5 3 

aero 2 o x x U x rel x. e 

o a o u 6 


Y = £ U 2 S[C + C 8+-C p-h-C r + C 6] 
aero 2 o 7 0 2 u Y p 2 u y r ^ T 

0^0 r 


: •= fu 2 S[C + C a + C u„ . +„£ C q +._£ C a+C 5 , 

aero 2 o z z U z rel 2U z 2U z. z.ej 
oaou oq o a o 

e 
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L 

aero 


t"o 2 

c. 


Sb[c 


+ c 


0 + ~ C 
2U i 


P + ~ c 
2U 


r + C 


fi r ] 


6 + 

V 


M 

aero 


a 

2 


IT Sc[C + C 


m 


m 


. I 

77 


a 


m 


u 


u 


rel 


+ 2F 


m 


q + 


— c 


2U 


m* 


a+ C 


a 


mp 


5 e> 


N 

aero 


O —? 

•r U Sb[C + C 
2 o n Q n 3 

+ C 5 ] 

n 6 r 
r 


3+ ~ C p+-C r + C 

2U n p 2U n r n <5 

o K o a 


In addition to these stability derivatives normally found in airplane 
equations, for VTOL aircraft it is desirable to add the rotational 
damping ooservea on a novenng airpiane. This can be expressed as. 

L aero VTOL L p P 

M aero VTOL * M q q 

N aero VTOL " N r r 

After summing the component parts, the complete non "linear equations of 
motion for a VTOL aircraft can be written as in Fig. A-I-3. 
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Figure A- 1 -3 

Summary of Complete Non-linear Equations of Motion for a VTOL Aircraft 


T 1 

u = vr - wq - g sin 9 + - [sin a cos X - ~ (u - u ) ] 

m v w 

j 

+ ~ pU 2 S [C +C a+^-C (u - u ) + C 5 J 
2ra 0 x rt x x w' x. e 

0 a 0 u 6 

e — 2 

T 1 ‘ pU 0 S 

v = wp - ur + g sin cos 9 + —[sin X - — (v -v ) ] + —z 

m v . w zm 

j 

[cl + c ^=- c Y p + V c v r + c v M 

0 x e 2U 0 X P 2U 0 r x « r 

T 1 

w " uq - vp + g cos <j> cos 6 [cos a cos X + (w - w ) ] 

m v_^ w 


(A-i-1) 


(A-I-2) 


PU - . 

+ — r [C -f C a + — - C (u - u ) + -rr r— C q 

2m 1 z n z z w 2U~ z 

0 a 0 u 0 q 


(A-I -3) 


+ — r c G 5 3 

21L. z - z. e 

0 a 6 

e 


(I V - > . I „ 2 . . . T 


qr + — ^ (r + pq) + y- [-sin A (z 2 + 2y cos a) 
x x 


U 2 Sb 


x 


+ V- c £ 

2U q P 


b 

2U 


- 

0 3 


■ — ^6 + p 

I X 

21 

“V 

X 

X 

r + C *. 5 a 


S r> 


(A- 1-4) 


(I - 1 ) I 7 ' 

q „ — A_ pr + (r _ p ) - A- (N.E.) I e H e r + — - 

y y 1 1 y y 


{cos A(z 2 sin a + x 2 cos a ) + — ■ - [x^ (w - w^) - z^ (u - c ; ) ]} 
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M, 


M 5 

q + 1 
y y 


+ y* q + <Sy + 


pu o 2sc 1 

[C + C a + ~-C (u-u) 

m IK m w 

y 0 a 0 u 


(A-I-5) 


+ ^r- C q + ^TT- C a + C 6 ] 

2U 0 q 2U 0 a 6 e 

n e 


(I - X ) I , 

r ° pq + y^ (p - qr) + — (N.E.) 1^ g q 

z z z 


x N 

+ y~ [sin * (x 2 - ^ sin o) - ~ (v - v w ) ] + ~ r 
z j z 


N - 2 

V p u/Sb , . 

+ -~6 + ■ ° [C + C 8 + — C p + — C r 
1 z 21 z n 0 n 8 2U n p 2U n r 


(A-I-6) 


+ c { :c 5 } 

n 5 a n 5 r 
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Figure A-I-4 


Summary of the Reduced Non-Linear Equations of Motion for the X-14 Aircraft 


T 1 

u - vr - wq - g sin 6 + — [sin a cos X - — (u - u ) 1 

m V , w J 

3 




p. off 


+ A C ) sin a - C cos a] 


v = wp - ur + g sin <J> cos Q + ~ [sin A - ~ (v - v^)] + ■— PU Q 2 SC 

j W ™ 


(A-I-7) 


(A- 1-8) 


• T 1 

w = uq - vp + g cos 4> cos 0 [cos 0 cos X + (w - w ) 1 

m V . w 

3 

- 2 ~ PU q 2 S [(C L + ACL) cos a + C sin a] 

p. off 


(A-I-9) 


a - ij i z 

P - )L — q r + -2L?. + p q) -f __ jL s i n X ( z ^ + o cos a) + — ^ ( v - v )] 

V v ^ * ’ i W 


X 


X 


j 


+ I 2 P + ~ 6 X + 2f- pU ‘Sb (C + c 5 > 
x x x * *0 *5 

a 


1 2 , 
L 

x 


<i - 1 ) i . 

z x xz , 2 2. 1 T 

q _ — ~ pr + — ( r - p ) _ — 2 I ft r + -A— cos X 

J- i e e I 


(A-I-1C) 


(z 2 sin a + x 2 cosa ) + ~* [x^ (w - w w ) - z (u - u tf ) ] 

M J 

M 5 

+ t * + S y + 2T pU 0 SC »V +1C B > 

y 7 7 0 „ m 0 

p« Oil 


(A-I-ll) 


+ (C_ + A C ) a + -rrr~ C q + C 6 } 

m 2LL m m , e 

a F<r a 0 q 6 

p. off n e 




(I - I ) I 

$• » — — ^ pq + (p - qr) + -j— 2 Ifl q 

z z z 

T x.. N 

+ ~ [sin A (x - £ sin a)--^(v-v)]+~r + 
£ _ 2 1 V. w I 


N 


. I z 

j z z 


+ ~-p U- 2 Sb (C + C A + C 6 ) 

21 z ° n 0 n 6 r n 6 3 

r a 


(A-i-12) 
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B. Application of the Equations of Motion to the X-1.4 Aircraft 

Several of the above aerodynamic stability derivatives have a very 
small effect in the representation of the dynamics of the X-14. The 
determination of which stability derivatives have this insignificant 
effect was based upon examination of wind tunnel data and a knowledge 
of the physical airplane. The final justification far retaining or 
neglecting different stability derivatives comes from the fact that an 
analytical model was obtained which when used in simulator studies closely 
approximated the physical airplane {5j« Using these findings the above 
aerodynamic contribution to the equations of motion can be expressed as 


X = ^ U _ 2 S (C T sin cx - C Q cos a) 


aero 2 o 


Y » u 2 s c v 

aero 2 o Y 

2 

Z = -pJJ S (C T cos a + C sin a) 

aero 9 o 1 D 

L - “ U 2 Sb (C 0 + C„ 6 ) + I, p 

aero 2 o £ , a' P 

a 

M «=^-U 2 Sc(C + C ot+ ~ C q + C <5 ) + M q 

aero 2o m in 2U m ^ m,- e q 

o a o q 0 


N U 2 Sb (C + C 6 +C 5 ) + N r 

aero 2 o n n„ r n„ a' t 


where C is the non-dimensional lift coefficient and C is the drag 

i) 

coefficient [ 2 ] . 

Continuing to follow the development in reference [5]* some of the 
above "power on" stability derivatives are expressed as "power off” plus 

V _ 

" 4” terms where the " A" represents the difference between power on and 
power off effects. The aerodynamic portion of the equations of motion 
now becomes 


X 


aero 


P_ 

2 


U/ S[«L 


+ A C ) sin ct ~ C cos ct] 


"power off 
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Y 

aero 



S C 


Y 


2 = - — U S [(C + AC ) ,cos a + C sin a] 

aero 2 o L . c L D 

power off 


L ® ~ U I 2 Sb (CL + C. 5 ) + L p 
aero 2 o & H, a p v 

0 <5 r 

a 


M = %■ U 2 Sc [ (C 
aero * o 


+ + (c + a c > « 

O „ * 1 „ % 

power off power off 


+ C q + C 6 ] + M q 

2 U m n % e q n 
o q o H 

* A 


ry 

N » 0 Sb (C + C <S+C 5 ) + N r 

aA v a n v* n r 


Upon incorporating these expressions for the aerodynamic forces in the 
complete equations of motion, the reduced non-linear equations of 
motion for the X-14 can be written as in Fig. A-I-4. 

These equations can be put in state variable form by eliminating 
the derivative terms from the right hand side of equations (A-I-10) and 
A- 1-12) . Rewrite the equations as 

I p = I r + L* 

X xz 

I r « I p + N* 

z xz r 


(A- 1-13) 



who rc 


L* = (I - I ) qr + I pq + T [-sin \ (z^ +£- cos a) 
y z xz i. 


+ h (v - v ' v + s p + x s x + 

+ 1/2 e D 0 2 Sb < C f 0 + % s a> 


N* = (I - I ) pq - I q r + 2 1 !! q 
x y xz n e e 

+ T [sin \ (x. - £ 1 sin o) - x, (v - v )/V.) + N 
/ x x w 3 r 


+ N f 6 + 1/2 p U 2 Sb (C + C 6 + C 5 ) 

6 z o n n, r n* a 

2 O O O 


Upon solving Eqns. (A- 1-13) 

p s- ( 1 , L* + 1 ^ K*)/(I 
r - (I L* + I N*) / (I 

XZ X 


r a 

for p and r it is found that 
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C. Numerical Evaluation of the X-14 Equations of Motion 


Numerical’ values for the constants in equations (A-I-7) to (A-I -12) 


are given in Table A- 1-2. In addition, the following stability 


derivatives can be expressed as functional relationships [5j: 

C. = 0.0030 cos 2 g/ deg 
* S 

a 

C - -0.00014 cos 2 g/ deg 
n 6 

a 

C ** 0.00118 cos^g/ deg 


C = -0.0178/ deg 
e 


Furthermore, Figures (A- 1-5) through A-I-12) show graphs of C. , 

*0 

Cyt C » C ? /l > 

x n 0 m L r r D m,. c 

0 a rc power off 0 ,, 0 

power off r power off 

AC. /T and T_vs. RPM. T is the non-dimensionalized thrust with T = 
L c 0 c c 

t 

T/(l/2 S) . T is the thrust at the engine exhaust and T is the 

U o 

thrust at the end of the diverters. Hence, 

1 ' T „ 

where £ is the longitudinal diverter efficiency factor and is the 
O A 

lateral side vane efficiency factor. Also, 

V V J 
J J 0 

which demonstrate that the jet velocity is somewhat less at the exit 
of the diverters than at the exit of the engines. In order to reduce 
these graphs to a useable form, functional relationships must be ob- 
tained. It was attempted to fit all of the graphs with polynomial 
curves. This, however, was impractical for the more irregular func- 


tions and straight line approximations were used. The functions used 


are shown In Table A-I-3. 
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The reaction control nozzle effects were evaluated by two independent 

methods. The X-14's response to sinusoidal inputs was analyzed by NASA 

at Moffet Field. The second method is to note that 

moment = M. 5 - Force . Length 

6 y 

y 

or 

<5 = F . length 

o y max 

y max 

where the length is the distance to the discharge nozzle. F^^ was 


measured directly using transducers, and <5 was also measured di- 

^raax 


rectly. Numerical values associated with these calculations and a com- 
parison with the other method are shown in Table A-I-4. This agreement 
is quite good, especially when considering the errors and uncertainties 
involved in measuring quantities like the moments of inertia and the 
maximum force at the nozzle outlets. 



Table A- 1-2 


11.4 
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C => -0.0178 deg 

m 

6 

e 


A 

m 

a 

= 0.576 rad _1 

M „ /I 

Sr y 

= 0.0350 rad/sec^/deg 

N /I 

4 z 

2 

= 0.0166 rad/sec /deg 

L /I 

& x 
0 X 

2 

= 0.0681 rad/sec /deg 


V, = 2000 ft/sec 

= °* 91 

= 1.00 



Table A-I -3 


T * T (£2) with £2 in percent of maximum RPM and T in pounds (for one 
engine only) 


T - 44.15011 * £2 - 1869.205 
o 

= 68.96552 * £2 -4184.483 
« 80 * £2 - 5278. 

■ 18.18182 * £2 + 910. 


£2 < 93.3 
93.3 < £2 < 99.1 
99.1 < £2 < 100.1 
100.1 < £2 < 101.2 


£2 - £ 2(1 ) 
o 


£2 = 0.02265 * T 4- 42.3375 
o 

« .0145 * T + 60.675 
o 

= .0125 * T + 65.975 
o 


T < 2250 
o 

2250 < T < 2650 
o 

2650 < T < 2730 
o 


= .055 * T - 50.05 
o 


ic L /? c = f <U 0 > 


2730 < T < 2750 
o 


AC t /T = -2.248398 x 10 7 IT 3 + 2.483009 x 10~ 5 U_‘ 
L c 0 0 


1.332149 


x 10 -3 U 


0 


A V T c = f <V 


AC^/T = 4 * 900354 x 10 7 U Q 3 - 3.164620 x 10 -5 U Q 2 + 1.529505 
x 10~ 3 u 


C D " W 


D 


.11 + l./(2. 152195 U Q + 12.5) 


V 


= C, 


M, 


0 


<V 


power off power off 

= -0.6 + 0.00237 l) 


power off 


« -0.52 + 0.0148 (U Q - 34) 


U Q < 34 (ft/sec) 
34 < U Q < 51 


« -0.27 + 0.001579 (U 0 - 51) 51 < U Q < 101 


0 
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: M " C M <0 0 J 

a power off °power off 


M 

a 


-4.1 + 0.04145 U 


0 


power off 


» -2.7 + 0.1007 (U Q - 34) 


- - 1.0 


Uq < 34 (ft/sec) 


34 < U r 




< 51 
> 101 


C Y = C Y <B) 


= -1.8518 



0 < g < 

0.189 (rad) 

« -0.35 + 0.393 

(6 

- 0. ] 89) 

0.189 < B < 

0.418 

= -0.26 - 0.362 

(8 

- 0.418) 

0.418 < 8 < 

1.22 

= -0.55 



8 > 

1.22 


C = C (8) 
n 0 n U 

C * 0.113 0 

n 

U 

= 0.022 

j 

= 0.022 + 0.1348 (8 - 0.594) 
- 0.13 


0 < 8 < 0.1945 (rad) 
0:1945 < 8 < 0.594 
0.594 < B < 1*395 

B > 1.395 



C « -0.083 8 
*0 

=-0.014 + 0.0066 (8 - 0.169) 
=-0.012 - 0.00865 (0 - 0.471) 
*=- 0.02 

C L = C L 

power off power off (U Q ) 

C = 0.61 + 0.002763 U 

power off 


0 < 8 <0.169 (rad) 
0.169 < 8 < 0.471 
0.471 < 0 < 1.395 
8 > 1.395 

U Q < 51 (ft/sec) 


= 0.75 



Table A-I-4 


Axis 

Discharge 

Point 

Length 

<5 

max 

F 

max 

Inertia 

Calculated 

Sensitivity 

NASA 

Sensitivity 

X 

Wing Tip 

16.9 ft. 

20° 

200 lbs. 

2 

2340 slug ft. 

2 

4.13 deg/sec /deg 

2 

3.9 deg/sec /deg 

Y 

Tail 

18.75 ft. 

ro 

o 

0 

130 lbs. 

2 

3400 slug ft. 

2 

2.05 deg/sec /deg 

2 

2.0 deg/sec /deg 

Z 

Tail 

18.75 ft. 

20° 

90 lbs- 

5400 slug ft. 2 

2 

.90 deg/sec /deg 

2 

.95 deg/sec /deg 


cn 

cn 
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D. Solution of the X-14 Equations of Motion for the Reference Control 
Variables 


The reference values for the elements of the control vector u_ are 
desired for any given state vector jc. Since a steady state reference 
condition is sought, set k = £. Upon examination of the non-linear 
equations of motion for the X-14 (see Fig. A-I-4), it is seen that 
u, v, and w are only functions of T, o and X . Therefore, an explicit 
solution for these variables should exist since there are three equations 
and three unknowns. To find this solution, first rewrite the equations 
as 


0 = A + T (sin o cos X - u ,/V. + f sin a) 

rel j 

0 = B f T (sin X - v _/V.) 

rel y 

0 • C - T (cos c cos X 4 - v , /V. f cos ~ ) 

rel j 


(A-I-14) 


where 


A s» ra (vr - wq - g sin 8) + 1/2 p S (C^ 


sin a - 


power off 


Cp cos a ) 


- 2 

B = m (wp - ur + g sin <$> cos 6) + 1/2 p S 


C = m (uq - vp + g cos <f> cos 6) - 1/2 pU^ S (C^ 


power off 


cos a + C D sin a) 


£ - 4 C i/ x c 

After squaring, rearranging and adding equations, the result can be 


written as 


DT + 2 ET - F - 0 


where 

D = 1 - 0 2 /V 2 - f 2 + 2f (u n sin a - w . cos a)/V. 

0 j rel rdl j 



57 . 


E ° (A u , + ® v + C w . ) /V - f (A sin a - C cos ft ) 
rel rel rel j 

2 2 2 
F * A + B + C 

Hence, 

f 2 

T » (-E + / E + DF )/D. 

Now since T is known, 

A = sin 1 ( v rel /Vj - B/T) 
and now from Eqn. (A- 1-14), 

0 - sin 1 [ C u re i/V^ - f sin a - A/X)/cos X ] . 


(A-l-15) 


(A-I-16) 


(A-I-17) 


(A-I—18) 


Some additional information is needed to solve the moment equations 
since there are six unknowns and only three equations. Since both the 
three reaction controls and the three aerodynamic controls have approximately 
the same range of travel 4 let us asyimp for determining the mot-mi <5PtHn« 
for reference flight conditions that the aerodynamic and the reaction con- 

J 

trols are of the same magnitude, or that 

6 =6 
x a 

6 =6 
y e 

6 = 6 
z r 

When the above numerically determined values of T, X and o are 

substituted into the summary equations for the X-14, the equations can 

be solved for <S , 6 , 8 , 6 , <5 , and 6 . The results are 
z e r x y z 

<5 =6 = - { (I - 1 ) qr + I pq + T [-sin X (z„ + l. cos a) 

x a y z xz L 2 1 

+ z i VkV C o Z Sb V } /( h 

J * Ox 

+ 1/2 p U Q 2 Sb ) 

6 

a 


(A-i-19) 
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6 

y 


5 — — { (I ~ I } pr + I (rr 2 - p 2 ) - 2 1 0 r 

e z x r xz v ee 


+ T [cos X (z 0 sin o + x,, cos a ) + (x, w , - z, u , ) /V . ] 

2 2 1 rel 1 rel j 


+ M q + 1/2 p U_ Sc [(C 
q U ra r 


+ A C ) 




power off 


+ (c 


m 


+ A C ) a + C q ] } / (M 


o 


m 


power off 


2U n m 
0 q 


(A-] 


+ 1/2 p U n Sc C ) 
0 m 6 


5 — 6 — — { (I - I)pq-I qr + 2 I « q + N r 
z r x y r xz ^ e e M r 

+ X [sin X (x„ -JL sin o) - s!, v ,/V.] + 1/2 P U 2 Sb (C 

1 rel j U n ( 

+ C <5 ) } / (N + 1/2 p U 2 Sb C ) 

n K a 5 * 0 n 6 

a r 


- 20 ) 


(A- 1-21) 
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APPENDIX B-I 

LINEARIZATION OF NON-LINEAR EQUATIONS OF MOTION 
FOR A TIME INVARIANT DYNAMIC SYSTEM 


The non-linear equations of motion for a time invariant dynamic 
system can be expressed as: 

A “ i ( “) (B— 1-1 ) 

where x is the state vector, u is the control vector, and f is a vector 
function of x. and u [ 1 ]. 

Suppose that for the system represented by Eqn. 05-1-1 ) there is a 
given _x (tg ) = Xg and u_ (t) for t > tg, then the unique solution 
is determined and can be found. Now consider small perturbations <$x (tg) 
in Xg and cfu (t) in u (t). Hence, small perturbations in the solution 
£ (t, Xg) can be expected and _x + - _f (_x + 6x, _u + 6 u); <$x (bg) ~ ^Xg 

< B-T-2) 

Expanding the right-hand side of Eqn. (B-I-2) in a Taylor series and 
retaining terms only of the first order, and after substituting equation 
A- 2-1 for x, we have the variational equation 


5x “ A 62 £ + B_ 5 u; 6x. (tg) = 6 xq 


where 


3f* 


ij 


9x 


j 


„ 9fi 

B . « i 

ij 3u, 

j 


<B-I-3> 

(B-I-4) 


with the partial derivatives being evaluated along the known solution, 
x (t) = $ (t, Xr.) fl nd control u (t). Hence, the motion of the dynamic 
system about a known path for a small perturbation is seen to be governed 
by the linear ordinary differential equations (B-I-3) and (B-I-4). 
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The numerical evaluation of the partial derivatives in equation 
(b-I~ 4) can be accomplished by using the central difference approxi- 
mation. [2] A review of this method follows. The Taylor series for 
the function f^ expanded about (x^ + A.x, u^) with 
(Ate) T = (0 , . . 0 x_. 0 ... 0 ) 


is 


3f . 


3 2 f. 


f i ( io + 4 £’ V ■ h <2o« V + + 1/2 + ••• 

3f . 3 

with t — evaluated at (x., u n ). Similarly, the function expanded about 
dX . — U — U 

j 


(Xq - Ax, Uq) is given by 


6f . 


3 2 f. 2 

f i ( 2o ' 2-' V ' f x ( V V " aT x j + 1/2 ~~2 4x j + ••• 

4 3x j 

3 

Since Ax is small, (Ax_. ) and higher order terms are neglected. If 
equation (B-l’-S) is subtracted £iom t£-I -5) , and the resulting equation 
solved for the partial derivative evaluated at the point (x^, Uq) with 
Ax^ - 0 for it ^ j , the result is 


( B-T-6) 


A. . = 

3k* 


3£. + Ax, uj - C (io - Ax, Ujj) 


2Ax. 


(b-I~7~) 


In a similar manner the elements of B are found to be 


B., = 
ij 


3f i f i ( V ^0 + A - } ' f i ( ^0’ ^0 " 


ou , 

J 


2(Au.) 

J 


(B-I-8 ) 


where (Au) =(000 


Au. 

J 


. 0) 



A digital computer program that performs the above operations has 
been developed and used in the linearization of the equations of motion 
of the X-14B airplane. 
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II. Dynamic Estimation of an Unmeasurable 
State for the X-14B VTOl Aircraft 



Dynamic Estimation of an Unmeasurable State 
for the X-14B VTOI. Aircraft 
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NOMENCLATURE 


MATRICES 

A = continuous time system matrix 

j3 - continuous time system input matrix 

C - continuous time system output matrix 

£ = modal continuous time output matrix 

G - controller gain matrix 

G = modal controller gain matrix 
*=m 

^ = modal continuous time input matrix 

= identity matrix 

K = full observer gain matrix 

L, » partial observer gain matrix 

7 = eigenvector matrix 

A = modal continuous time system matrix 

$ = discrete time system matrix 

$> , = desired discrete time modal matrix 
*=a 

= modal discrete time system matrix 

* c tn 

2 = discrete time input matrix 
v = modal discrete time input matrix 


VECTORS 

u_ - input vector 
= eigenvector 
x ~ state vector 
~ output vector 
z = modal state vector 
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NOMENCLATURE 


SCALARS 

q = pitch rate 

r = yaw rate 

t = time ' 

u * forward velocity 

v = side velocity 

w = vertical velocity 

6 = pitch angle 

X ~ continuous time eigenvalues 

A^ - desired eigenvalues 

(A)pis - discrete time eigenvalues 

t = time interval for difference equations 

<j> = roll angle 
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CHAPTER 1 
INTRODUCTION 

Accurate measurement of the complete state vector is an integral, 
and sometimes difficult, part of implementing the powerful techniques of 
modern control theory. This research is concerned with the development 
of a method to estimate the side velocity of NASA's X-14B vertical takeoff 
and landing vehicle during hover. When the aircraft is flying at low side 
velocities, particularly during hover, measurement devices are unable to 
give reliable readings, and thus are usually not used. Hover is defined 
as low speed flight where aerodynamic and inertial effects are not impor- 
tant. This corresponds to flight with forward velocities of less than 
30 ft/sec. Determination of the side velocity is important for 
both completing the state vector for control and for general knowledge 
of flight. 

The approach is to investigate the application of linear dynamic ob- 
server theory around flight operating conditions to estimate the side 
velocity. Both full and partial dynamic observers are developed for 
two linearized hover conditions. 

The linear differential equations are discretized to hold computer 
time to a minimum. The observers are tested on the linearized systems. 

A linear modal control policy is also developed primarily to test the ob- 
servers’ ability to track the real system. Finally the developed linear 
observers are applied to the full nonlinear equations to determine their 
accuracy and range around the given operating conditions. 
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CHAPTER 2 
THE MODEL 

The model for the X-14B VTOL was developed in previous research by 
Hoffman, Loscutoff and Seevers [1} and has recently been improved by 
Roesener [2], The resulting eight nonlinear differential equations are 
found in Figure 1. The state variables for the model are 
6 = pitch angle, 
q = pitch rate, 
u = forward velocity, 

<f> = roll angle, 
p = roll rate, 
r - yaw rate, 
v = side velocity, 
w = vertical velocity. 

This set of nonlinear differential equations was linearized around several 
operating conditions characterizing hover, transition and full -flight 
in [!]. 

In this research, the system is linearized at the following operating 
conditions : 

Case 1 : u = .01 ft/sec, 

6=q = w= 4> = p = r = v = 0 

Cas e ?. : u = 20 ft/sec, 

e-q=w=<£=p = r= v = 0 . 

T'-iese conditions are used for development of linear observers and control- 
ler, They will be referred to as Case 1 ar.d Case 2 respectively m sub- 
sequent sections. The linearization at u = .01 ft/sec was chosen oecause , 



6P. 


it was believed to be the most difficult flight condition for estimating 
the side velocity, due to the limited coupling between the measurable 
states and the side velocity. The u = 20 ft/sec flight condition was 
chosen because it is a typical hover condition. These linearized flight 
conditions are found in Tables 1 and 2. The output matrix, in Table 1A, 
is the same for both Case 1 and Case 2. The output matrix shows that all 
states are directly measurable, except for the side velocity. 
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FIGURE 1 

.Summary of the Reduced Non-Linear Equations of Motion for the X-14 Aircraft 
.u ** vr - wq - g sin 9 + - [sin a cos \ - (u - u )] 

m Vj W J 

1 2 

* "2m p S ^ + 4 C ) sin a - C cos a ] 

^ ° V off L D 

* - wp - ur + g sin * cos 6 + I [sin \ - ( v - v^j + PU^SC^ 
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T :J " 

*" 2 m pU o^ S - ^ C r + icL ) cos a + c_ sin a] 

p. off u 

- (I V ~ V \z . T Z, 
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.f ^ q ^ 6 + ■— pU 2 Sc [(CL +AC ) 

y y y 2I y 0 X m o 

p. oxr 


+ (C 


m 


p# off 


+ A C ) a + •—* C q + C { ] 
tn 2U n m - ~ 1 


J r\ 

0 q 


m 6 e 


(I ** I ) I 
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N 
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TABLE 1 


THE LINEARIZED SYSTEM - CASE 1 


REFERENCE CONDITION 

THETA =0.0 RADIANS 
Q = 0.0 RAD/SEC 

U = 0.01 FT/SEC THRUST 

W =0.00 FT/ SEC 
PHI = 0.0 RADIANS 
P = 0.0 RAD/SEC 
V = 0.00 FT/SEC 
R = 0.0 RAD/SEC 


TRIM PARAMETERS 

ELEVATOR*= 0.47 DEGREES 
REACTION NOZZLE DY = 0.47 DEGREES 
DEFLECTION ANGLE (SIGMA) = 0.00 DEGREES 

ENGINE RPM =98.57 PERCENT OF MAXIMUM 
AILERONS*= 0.00 DEGREES 
RUDDER*= 0.00 DEGREES 
SIDE VANE ANGLE (LAMBDA) = 0.00 DEGREES 
REACTION NOZZLE DX * 0.00 DEGREES 
REACTION NOZZLE DZ = 0.00 DEGREES 


A MATRIX 


0.000000 

1.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

-0.150113 

-0.000436 

-0.000147 

0.000000 

0.000000 

0.000000 

-0.500992 

32.035897 

0.000000 

-0.020130 

-0.038363 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.010000 

0.000050 

-0.020319 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

1 .000000 

0.000000 

0.000000 

0.000000 

0.024327 

0.000000 

0.000000 

0.000000 

-0.451157 

0.000306 

-0.015424 

C. 00 1000 

0.000000 

0.000000 

C. 000000 

31 .832373 

0.000000 

-0.022425 

-0.010000 

Q.OCOOCO 

0.316251 

0.000900 

0.000000 

0.000000 

-0.015039 

-0.002452 

-0.200514 
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MATRIX 





0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 



0.034907 

0.019679 

0.001822 

0.000000 

0.000000 

0.000000 



0.000000 

0.562003 

0.000000 

C. 000000 

0.000000 

0.000000 



0.000000 

0.000004 

3.600197 

0.000000 

0.000000 

0.000000 



0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 



0.000000 

0.000000 

0.000000 

-0.033893 

0.068243 

0.001279 



0*000000 

0.000000 

0.000000 

0.562003 

0.000000 

0.000000 



0.000000 

0.000000 

0.000000 

-0.001310 

0.002275 

0.016623 



♦NEGLIGIBLE AT LOW SPEEDS, THEREFORE NOT INCLUDED AS INPUTS IN THE B MATRIX. 



TABLE 1 A 


THE OUTPUT MATRIX - CASE 1 AND CASE 2 


C MATRIX 


1.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

1.000000 

0.000000 

0.000000 

0.000000 

0,000000 

0.000000 


0.000000 

0.000000 

1.000000 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

1.000000 

0.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 



TABLE 2 


THE LINEARIZED SYSTEM - CASE 2 


REFERENCE CONDITION 


TRIM PARAMETERS 


THETA =0.0 RADIANS 
Q = 0.0 RAD/SEC 
U =20.00 FT/SEC 
W = 0.00 FT/SEC 
PHI =0.0 RADIANS 
P = 0.0 RAD/SEC 
V = 0.00 FT/SEC 
R = 0.0 RAD/SEC 


ELEVATOR 
REACTION NOZZLE DY 
THRUST DEFLECTION ANGLE (SIGMA) 

ENGINE RPM 
AILERONS 
RUDDER 

SIDE VANE ANGLE (LAMBDA) 
REACTION NOZZLE DX 
REACTION NOZZLE DZ 


-1.96 DEGREES 
-1.96 DEGREES 
0.87 DEGREES 

98.76 PERCENT OF MAXIMUM 
0.00 DEGREES 
0.00 DEGREES 
0.00 DEGREES 
0.00 DEGREES 
0.00 DEGREES 


A MATRIX 


0.000000 

0.000000 

-32.035897 

0.000000 

0.000000 

0.000000 

0.000000 


I. 000000 
-0.375150 
0.000000 

20.000000 
0.000000 
0.024373 
0.000000 


0.000000 

-0.001513 

-0.028319 

-0,026590 

0.000000 

0.000000 

0.000000 


0,000000 

-0.014664 

-0.007426 

-0.024535 

0.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

31.832873 


0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

-0.451157 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

-0.004610 

-0.082783' 


0.000000 

-0.501932 

0.000000 

0.000000 

0.000000 

-0.015424 

- 20.000000 


B MATRIX* 


0.000000 

0.037416 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

0.019777 

0.564681 

0.008554 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

-0.015868 

-0.009580 

3.552470 

0.000000 

0.000000 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

-0.034060 

0.564746 

-0.001351 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.072015 

0.000000 

0.002325 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.001329 

0.000000 

0.017267 


*(DY-ELEVAT0R) .(DX+AILERONS) ,AND (DZ+RUDDER) ARE COMBINED AS INPUTS IN THE B MATRIX. 
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CHAPTER 3 

CONTROLLABILITY AND OBSERVABILITY 

Controllability and observability are dual criteria that must be 

satisfied to allow feedback control and state estimation. Controllability 

is the ability to transform any state at time to any other state at 

time t. by applying the appropriate control input jj(t). Observability 
0 

is the ability to identify any previous state by observing the output 
jy(t) for a finite time interval. 

Several methods for determining whether a system meets these criteria 
exist. The method used in this research determines controllability and ob- 
servability of individual modes. This is a superior method since modal 
input-output signal flows can be determined and the degree of controlla- 
bility and observability of modes can be determined. For example, if a 
system is found to be marginally controllable or observable it is important 
to know if the offending mode or modes are stable or unstable and to what 
degree. 

The linearized system equations 

x = Ax + B u ( II- 1 ) 

y_ = Cx { 1 1— 2 ) 

are transformed into the modal domain through the similarity transformation 
x_ = T 2 

where T is a matrix of the system eigenvectors _v such that 

J ~ f.¥] Ip . . . v^] . ( 1 1 - 3 ) 

The system equations become 
Z = Az i Hu 

y. " 


( n-4) 
(H-5) 
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where 

A = T^AT 

H = T-'b 01-6) 

F = CT . 

A is the modal system matrix, 



where X's are eigenvalues. 

In the event that complex roots are present, the T matrix takes the 

form 

t-IR -IC ’ * * J 

where 

. v_ IR = the real part of the eigenvector 

v^£ - the complex part of the eigenvector. 

an ^ are C0 ' n P^ ex conjugates. The resulting modal matrix. A, 
becomes 
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where 

Xj-p = real part of eigenvalue 
= complex part of eigenvalue. 

The treatment and interpretation of systems with complex eigenvalues 
is the same as for any other system except for the variations noted. The 
case of repeated eigenvalues is neglected throughout this discussion. 

H is the modal input matrix. If any complete row of H is zero, the 
mode is unaffected by any control policy implemented, and is uncontrollable. 
Therefore, the controllability criterion is satisfied if every row of H 
has at least one non-zero element. 

Observability can be determined by examining the modal output matrix 
JF. The requirement is that every mode must appear in the output vector 
To accomplish this, every column in £ must contain at least one non-zero 
element. 

For Case 1, the A, H and £ matrices are as shown in Table 3, while 
for Case 2, the A, £ and £ matrices are shown in Table 4. From both A 
matrices, it can be seen that the linearized systems are inherently un- 
stable. However, both linearizations are completely controllable and com- 
pletely observable, and therefore can be stabilized by feedback control. 

The most important result is that since the side velocity is observable 
through the measurable state variables, dynamic observers may be used to 
estimate it. 



TABLE 3 


THE LINEARIZED SYSTEM IN MODAL FORM - CASE 1 

A MATRIX 


-0.13660 

0.42355 

0.00000 

0.00000 

-0.42555 

-0.18650 

0.00000 

0.00000 

0.00000 

0.00000 

-0.30996 

0.020 1 9 

0.00000 

O.COOOO 

-0.02019 

-0.30995 

0 . 00000 

o.oocco 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

o.ocooo 

0.00000 

0.00000 

O.OOOGO 

0.00000 

0.00000 

0.00000 

0.00000 


3.33837 

1.97663 

0.31828 

3.31723 

1.34020 

0.06916 

-2.18415 

-1.30327 

-0.34636 

-11.50003 

-6.75670 

-1.41944 

4.30468 

2.72745 

1.61925 

2.01065 

1.05477 

0.04308 

-3.50121 

-1.78538 

-0.13110 

0.00431 

0.00251 

3.86946 


0.00520 

-0.01328 

0.00905 

-0.00063 

0.00468 

0,00469 

-0.00279 

0.00038 

1 .00000 

0.00000 

1.00000 

0.00000 

0.00002 

-0.00023 

-0.00008 

-0.00002 

0.00009 

-0.00078 

0.00383 

-0.00092 

0.00035 

0.00011 

-0.00117 

0.00036 

0.00345 

-0.00364 

-0.00175 

0.00023 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

-0.12060 

0.00000 

0.00000 

0.00000 

0.00000 

0.15912 

0.00000 

0.00000 

0.00000 

0.00000 

0.11029 

0.00000 

0.00000 

0.00000 

0.00000 

-0.02084 

I MATRIX 

-2.88713 

5.24580 

1.94411 


-3.45986 

6.14046 

-1.88379 


-65.27823 

127.52590 

1.04069 


-289.15382 

569.48136 

5.42458 


100.33871 

-197.84024 

-7.27932 


-44.71371 

94.88475 

-0.54612 


-76.76226 

159.69359 

3.74697 


0.05583 

-0.11182 

-0.00548 


' MATRIX 

0.00314 

0.00558 

-0.00407 

-0.00112 

-0.00038 

0.00089 

-0.00045 

0.00002 

0.99977 

-0.99772 

0.99973 

-0.36601 

-0.00046 

-0.00023 

0.00035 

0.93061 

0.00006 

0.00038 

0.00009 

0.00000 

- 0.00001 

0.00006 

0.00001 

0.00000 

-0.00085 

0.00032 

-0.00064 

0.00005 
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THE LINEARIZED SYSTEM IN MODAL FORM - CASE 2 

A MATRIX 


-0.28951 

0.66661 

0.00000 

0.00000 

-0.66661 

-0.28951 

0.00000 

0.00000 

0.00000 

0.00000 

-0.74830 

0.00000 

0.00000 

0.00000 

0.00000 

-0.17565 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


1.82635 

1.03373 

3.08609 

1.41645 

0.67930 

-2.88325 

0.12454 

0.07176 

0.18701 

0.12258 

0.11941 

1. 09 256 

-1.94683 

-0.58549 

-4.16373 

-1 .91553 

-0.79425 

22.85470 

-0.24056 

-0.12457 

-0.22805 

-0.18510 

-0.08084 

0.S5420 


0.00811 

-0.02072 

0.00020 

0.00436 

0.01146 

0.01140 

-0.00015 

-0.00077 

1 .00000 

0.00000 

0.00911 

0.96188 

0.19109 

-0.38004 

0.00455 

0.27063 

0.00103 

-0.C0059 

0.02073 

-0.00649 

0.00009 

0.00086 

-0.01551 

0.00114 

0.00459 

-0.00607 

-0.00027 

-0.01050 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.06400 

0.14200 

0.00000 

0.00000 

-0.14200 

0.06400 

0.00000 

0.00000 

0.00000 

0.00000 

0.10626 

0.43161 

0.00000 

0.00000 

-0.43161 

0.10626 


H MATRIX 


18.60148 

-0.03779 

-0.00078 


-18.16568 

-0.05271 

-0.00090 


-2.63280 

-2.52081 

-0.05069 


-40.15281 

-0.02941 

-0.00043 


20.09309 

-0.01219 

-0.00002 


65.47933 

0.08503 

0.00160 


-6.41398 

-2.51054 

-0.03326 


-11.89252 

4.99281 

0.09281 


MATRIX 




-0.00285 

-0.00444 

0.00023 

0.00016 

0.00045 

-0.00069 

-0.00004 

0.00012 

1.00000 

O.COOOO 

-0.01555 

0.01246 

-0.12565 

0.04577 

0.00397 

0.00224 

0.00000 

-0.00057 

0,00598 

0.01345 

0.00008 

-0.00004 

-0.00517 

0.00401 

' 0.00007 

-0.00086 

0.00007 

-0.00018 
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DISCRETIZATION OF LINEAR DIFFERENTIAL EQUATIONS 


The linear differential equations are transformed to discrete time 
equations because the aircraft has an on-board digital computer for con- 
trol. Thus discrete time control is more applicable. Digital computers 
are also more suited to adding and subtracting difference equations than to 
numerically integrating differential equations [3]. 

Difference equations are the time domain solutions to differential 
equations. Therefore, the set of differential equations (1) and (2) 

x = Ax + Bu (II-l) 

£=Cx ( 1 1— 2 ) 

can be transformed into the difference equations 

*k-H ' &k * ffifc > (n ' 7) 

y.*£Sk in-8) 

Assuming iu(t} is constant over a time interval (0 < t < x, where i is 
tne time iureivtO), it can be shown that for linear, constant, coefficient 
systems, $ = e^ T (1 1-9) 

(11-10) 


T = [e- T - I1A - 1 _S . 


It will now be shown that $ can be simply formed from A when in the 
modal d;r«.«-it>. first, it must be shown that the similarity transformation 
matrix, 1, regains invariant and continues to diagonalize the system 
whether in difference or differential equation form, 
is described by the series 


0 - e-- T - l + At + I,A ? t 2 + ^,A 3 t 3 + 


( 11 - 11 ) 

To show that 1 diagonalizes c, premultiply and postmultiply by T'^ and T 
respect i i ely . 


-1 




1 t - 1^,2 


i $T = i -JI + X '^It ’A lx + .. 


I H-12) 
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Now, simplify (12) using the identity 
A 2 = MI“ ] A , 

ear i 

T"V * I + At + 1 A 2 t 2 + ... 
«= = = = 2 != 


(H-13) 


This technique can be readily extended to all higher order terms; there- 
fore 


T -1 $T = * * e^ T 

“ *=“ =m 


e X l r 0 . . . . 0 

0 e X 2 T . . 0 


(11-14) 


The discrete time eigenvalues are defined as 


X DIS 


1 2 2 

+ At + ^-,A t + 



(11-15) 


in equation (13) . $ can then be obtained from h by the following steps. 
First, transform the continuous time system into the modal domain, a step 
done previously to determine controllability and observability. Then 
create by transforming the individual eigenvalues of A into by 
using the scalar equation (15). Now, simply transform out of the modal 
domain, 

I^TST 1 . (11-16) 

— -nr 

The input matrix, ¥, can be approximated easily by noting that 
I = [(I + Ar ■» -|,A 2 t 2 + ...) - 

I - T[A + i (11-17) 

I « xLl + |-,At]B . 


For this system, the time steo was .05 seconds. For this second order 
approximation the error is cn the order of 0.25 oercent when the higher 


order terms are left out. The system equations for Case 1 are found in 
Table 5 and those for Case 2 are in Table 6. 
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TABLE 5 


THE LINEARIZED DISCRETE TIME SYSTEM - CASE 1 


l.CCCOO 

0.04980 

0.00000 

0.00000 

0. 00002 

0.99200 

-C.00CQ2 

-0.00001 

v * OUv/L. 

-C.CX1Q 

O .99 SO 0 

-0.00192 

c.orco 

0.0 X 50 

** S\ 

V/ * w 

0.S99C0 

0.0 X 00 

0.C5C03 

0,00000 

o.oocoo 

o.ocoqo 

0.G0119 

0.00000 

0.00000 

0.00050 

0.01570 

O.OOGOO 

0.00000 

0.00000 

0.C0005 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 

0.00175 

0.00098 

0.00009 

0.00000 

0.02810 

0.00000 

0.00000 

o.coooo 

0.18000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


MATRIX 


o.coooo 

0.00002 

-0.00062 

0.00000 

0.00000 

o.oooco 

-0.02*80 

0.00002 

0.00912 

0.00202 

0.00040 

-0.00001 

w , vO Job 

0.00000 

-0.0 X 0 1 

0,00000 

1.0 CC 00 

0.04940 

-0.00002 

0.00000 

0. C 0001 

0.97300 

-0.00077 

0.00 QC 2 

-0. C 0010 

-0.09075 

0,99000 

-0.0 C 012 

1.59090 

0.03870 

-0,00043 

0.99900 

MATRIX 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


-0.00169 

0.00341 

0.00006 


-0.00007 

0. 00011 

0.00083 


0.02810 

0.00000 • 

0.00000 




TABLE 6 


THE LINEARIZED DISCRETE TIME SYSTEM - CASE 2 


4> MATRIX 


1.00000 

0.04950 

0.00000 

-0.00002 

0.00006 

0.98100 

-0.00007 

-0.00073 

•1 .60000 

-0.04000 

0.99900 

-0.00036 

0. 00110 

0.99000 

-0.00137 

0.99800 

0.00000 

0.00003 

0.00000 

0.00000 

0.00000 

0.00119 

0.00000 

0.00000 

0.00000 

0.01560 

0.00000 

-0.00001 

0.00000 

-0.00782 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 

0.00187 

0.00099 

-0.00079 

0.00000 

0.02820 

-0.00048 

0.00094 

0.00092 

0.17700 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

-0.00062 

0.00000 

0.00000 

0.00001 

-0.02470 

0.00000 

0.00000 

0.00000 

0.00041 

0.00000 

0.00000 

0.00000 

-0.01240 

0.00000 

1.00000 

0.04940 

-0.00002 

-0.00001 

-0.00018 

0.97800 

-0.00066 

-0.00023 

0.00000 

-0.00074 

0.99000 

0.00000 

1.59000 

0.03980 

-0.99300 

0.99600 

MATRIX 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


-0.00170 

0.00360 

0.00007 


0.02820 

0.00000 

0.00000 


0.00003 

-0.00006 

-0.00043 
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CHAPTER 5 


LIMITED INPUT MODAL CONTROL VIA STATE VARIABLE OUTPUT 


The linearized set of difference equations, given by equations . (?) 
and (8) and repeated here, 

(II-7) 

( I 1-8) 


J< k+ , * % + Xdk 


ik ' .&k ’ 

specify a system of eighth order, with six inputs and seven outputs. For 
the purposes of controller design, though, state variable output will be as 
sumed. The missing states will be provided later by dynamic observers. 

The control matrix G will be designed to place six of eight eigenvalues 
wherever desired.- The remaining two eigenvalues will be considered suf- 
ficiently stable, although a readily available extension of the following 
technique allows for placement of all eight poles. 

To condition the first six eigenvalues, the system is transformed into 
the discrete time modal domain, using 


Thus 




= 4 ^ + 




= CTz, 


(11-18) 

(11-19) 


where 


w = 

^ ™ “ 

s = T -1 $T . 

=m *= = 

Now, let us partition the system equations into two segments. The modal 
state vector becomes 

_z = the first six modes, all to be conditioned 


: 78 


= the last two modes, not to be conditioned. 
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The $ matrix is arranged so that 
“in 


0 


$ = 
=m 


*Vdis 

0 ^Vdis 
! (xJ 


0 . 


r- 8 


3 y DIS 


^ VdIS 


^ni 78 


M 


5 / DIS 




6 y DIS 


^Vdis 

(Xq) 


8'DIS 


J^ 7 g consists of the last two poles, which are the most stable and there- 
fore are ignored in the conditioning process. 

Let us define a matrix of desired eigenvalues, 




^ X 1d^ DIS 0 * 
0 ( A 2d^DIS 


0 . . 


^ X 3d^DIS 

^Wdis 


^ X 5d^ DIS 


^ X 6d^DIS 


0 


(X 7 ) 


V DIS 


^8^ DIS 


=^dL 


% 8 i 



These eigenvalues will be the closed loop system poles after implementa- 


tion of the control matrix G. Partition T such that 

■=* “Hi 


¥ = 
“in 


¥ 


■±m78 


where is a nonsingu1ar^6x6 matrix. Then the control matrix in modal 
form is 

Sm ' [<4 l ^ - &dL £>] (II - Z0 > 

to condition the first six modes. To illustrate that this will give the 
desired results, let us examine the system in modal form under the control 
policy 


A “ -Vk 


We have 


^k+1 ^m^k 


- T f”, 1 f [j> , 0] - Q 0] z, 

^mmiiL I ^miu =" - di. ~= j — ? 

Simplifying and rewriting in a partitioned form 

® ~ ^dL ^ 


(11-21) 


(11-22) 


» * 

—16 


r 

^mL 

0 1 


■ x * 

r 

—78 

k+1 


_0 

4i78- 


If/ w 

LJ W£HmL -J 

- 


h* 


■7 

-16 

. 


—78 

* 


j. * 


Expanding, we obtain 


% 


f— 16 

’ 

i 

In 

L® 

k+1 . 


£ 


W^[[1<IL ^ - ^nL £ ] ] % 


1 

r» «• 


—16 

A 

/78_ 


k 

(11-231 
( 11-24) 


Note that despite the f ] term ’ the modal c1osed loop system 


matrix 
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* 

f^ld^DIS 0000 

0 

0 

o 

t 

0 ^WdIS 000 

0 

0 

0 

0 0 ^ A 3d ) DlS 0 0 

0 

0 

0 

0 0 o (x 4d ) Dls 0 

0 

0 

0 

-0 0 0 0 (A gd ) 

DIS 0 

! o 

0 

0 0 o 0 - 0 


DIS 0 

0 

4ya4L^dL ' 


* { Vdi$ 0 
« 

^Vdis 

• J 


is in triangular form and its eigenvalues are the six elements of and 
the original eigenvalues and (Ag)^. Thus, the first six modes 

nave been moved to their desired values, and the last two modes have re- 
mained unchanged. To transform jG into the control matrix, reverse the 
modalizing process: post multiply by T~^ , 

-I 


— “ml.'^nL -dL -- 


( II-25) 

The resulting closed loop equations are 

x k+1 r (1 ~ TG)x^ . ( II -26) 

Using this technique the control matrices derived for Case 1 and Case 2 are 
as shown in Table 7. 

The modal technique oVlo/s the designer to place eigenvalues arbitrar- 
ily. The onl-y restriction is that eigenvalues must be chosen to give the 
required response while remaining »'itliin control power restraints. 

It should be noted here from previous development that 
Xx 


X DIS = c 


where 


X.. T( . ~ disci cui tin., eigenvalue 
Dio 

X = conti nucr. iir-ie eigenvalue 


x = time int' i ol . 
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For the stability limit on the complex plane is the unit circle. 

Positi ve, real parts of magnitude less than one correspond to continuous 
time negati ve, real parts. Zero corresponds to infinity as the fastest 


response. Any eigenvalue within the real, negative half of the unit circle 
will oscillate with each time step, and thus, is not very useful -for con- 
trol. Time constants for discrete time eigenvalues can be approximated as 



( 11 - 27 ) 


by recalling the infinite series for an exponential function. 
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TABLE 7 

THE CONTROL MATRICES FOR THE LINEARIZED DISCRETE TIME SYSTEMS 

CASE 1 
G MATRIX 


139.00000 

72.00000 

-0.70000 

1'. 74000 

-4.24000 

9.70000 

24.30000 

0.66200 

-239.00000 

-83.10000 

1.24000 

-3.44000 

25.20000 

3.38000 

-53.00000 

-1.20000 

0.27000 

0.30000 

0.00017 

2.77000 

-1.70000 

-3.9300 

-0.48400 

0.00114 

27.10000 

1I.70C00 

0.00479 

-0.40100 

127.00000 

76.3000 

12.30000 

1 . 38000 

12.80000 

5.50000 

0.00217 ■ 

-0.18900 

67.30000 

48.10000 

4.91000 

0.68800 

6.22000 

12.70000 

0.00877 

-0.05930 

13.7C000 

5.04000 

49.30000 

-0.02310 





CASE 2 








G MATRIX 




-2.16000 

0.00000 

-0.37500 

0.04670 

25.40000 

0.00000 

-0.03490 

-0.00006 

25.00000 

0.00000 

0.70800 

0.00156 

-38.50000 

0.00000 

0.10400 

0.00019 

-2.02000 

-0.05610 

-0.00171 

0.11300 

-0.27000 

-0.00117 

0.00028 

0.00000 

-0.152C0 

0.00C04 

0.00000 

0.00000 

-0.20000 

-0.00085 

0.70800 

0.00000 

7.24000 

0.00260 

-0.00003 

-0.00018 

-13.7000 

5.48000 

0.34700 

0.85700 

-379.00000 

-0.13100 

0.00182 

0.00968 

1050.00000 

3.67000 

-0.61400 

-46.50000 
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CHAPTER 6 


STATE VARIABLE OUTPUT VIA DYNAMIC OBSERVERS 


To implement the control scheme discussed previously, a full state 
vector is assumed for feedback. On the X-148 VTOl, seven of eight state 
variables are available for measurement during hover. A dynamic observer, 
which is a special form of filter, is used to estimate the missing state, 
side velocity. 

The full dynamic observer is a time invariant linear model of the 
system to be observed. The model is driven by the same inputs as the 
system. Therefore, if the model is perfectly accurate, i.e., no dis- 
turbances or noise are present and the initial conditions of the system 
known, the system state vector is known at any time t. As can be readily 
seen, these arc rather demanding assumptions and a method for updating the 
model states to the measurable system states is needed. This is the key 
to the usefulness of the dynamic observer. 

To design a dynamic observer, assume a system is governed by equa- 


tions (7) and (8), 


Let us define ^ as an estimate of 


The model 


defined by £ is 


(H- 28 ) 

£=£4 • (11-29) 


This is tne open loop model of the system. To force tiie model states to tne 

A 

system states an error feedback term, j((y. - ^ ) , is implemented. K is 
the observer control matrix. Thus, the closed loop model becomes 

" .M}- + ,2^ + ' 11-30) 




(11-31) 

To illustrate that a control matrix, K, can be chosen to drive 
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observer error in measurable states to zero, let us define an error state 
vector 

X = X - X . 

Subtracting equations (7 ) and (8 ) from (30) and (31), 

4+i ' 4+i ’ *4 + *4 + &4 ' 4> ' *4 - XHk 

4 -4 “4 

which becomes 

4+i = *4 + W4 ' 4> 

4=£4 

or the single set of equations 

4 +1 = (i - * <H-33) 

The eigenvalues of this homogeneous set of differential equations are 
determined by the arbitrarily chosen j<. By choosing K to make this set of 
equations stable, it can be seen that 
x(t) ->0 as t ■+ «° 



The rapidity of convergence depends solely upon the error system's 
closed loop eigenvalues determined by K. 

The dynamic observer, then, models the system, continually forcing 
the measurable states' error to zero at a convergence rate determined by 
J<, and gives an estimate of unmeasurable states based upon the measurable 
states. Thus, a dynamic observer, supplying estimates of unmeasurable 
states, can be used in conjunction with control schemes based on state 
variable feedback to create control. A block diagram o f the complete 
closed loop system can be found in Figure 2. 
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The matrix K can be determined by any method desired by the designer. 
This system, by the simplicity of its output matrix £, lends itself to a 
very simple procedure. First notice that £ is of the form [£ j 0]. The 
homogeneous error system matrix becomes 

- KC) = {$ - [K j 0]) . 

The effect of K upon the closed loop system, in this case, is clear. 

Choose 

K ij = $ ij 1 t j ; i = 1, 2, ... , 8 

K • • - $ • • *■ A j i = j ! j - 1 , 2, . . . , 8 
J ij 

where A^ is the desired eigenvalue, and the convergence rate is deter- 
i j 

mined. All but one eigenvalue can be placed using this method. To com- 
plete the arbitrary placement of all eigenvalues, determine the state whicn 
is most strongly coupled to the unaffected vector in J>. Then, devise a 
second order system with the two states. For example, for Case 1, it is 



(11-33) 


Obviously, K 77 


and Kgy can be chosen to arbitrarily set the eigenvalues 


of this subsystem. Since all of the other states are already decoupled. 


the eigenvalues of t.ie total system will be the main diagonal plus the 
eigenvalues determined by i.he second order subsystem. 

This design method v<ss used to calculate observer gain matrices for 
both linearized cases, i or Case 1, the only coupling between the side 
velocity vector and the rest ef ~.he system was through yaw rate. The re- 
sulting K matrix is shown m Table 8. For Case 2, the major coupling was 
between side velcuty r.ns 1 roll iate. The resulting observer gain matrix 


is also shown m i able 8, 
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The results of applying the dynamic observers in conjunction with the 
control systems described earlier are illustrated in Figures 3 and 4. 

These simulations are done by modeling the linearized systems, setting 
initial conditions at u = 10 ft/sec, v = 4 ft/sec, and all other state 
variables at zero. Recall that these velocities are relative to the 
linearized operating conditions. The observer started with all initial 
conditions at zero, so that for Case 2, u is actually 30 ft/sec. Since 
the system itself, in the computer simulation, was a model, a comparison 
between the system's siae velocity and the observer's estimated side velo- 
city is possible. Time constants for both observers are around one-fifth 
second and by one-half second have completely converged. This is a very 
rapid rate of convergence. The observer converges to the other state vari- 
ables, which are measurable and directly controllable as fast as or faster 
than the unmeasurable side velocity. It is interesting to note that side 
velocity is much easier to control for Case 1 (T.C.. - 1.75 sec) than for 
Case 2 (T.C. * 6 sec). The inertia terms in Case 2 couple the system more 
strongly, but also make the control more difficult. 

The full observer for Case 1 was also run with noise introduced into 
measurement (Figures 5 and 6). The measurement error was randomly intro- 
duced ±1.0 ft/sec for u and +10 percent for all other state varialbles. The 
estimated forward velocity follows the actual forward velocity much better 
than the measurement of u. This is due to the filtering effect of a full 
observer. Since an observer converges to a measurable state at a rate 
determined by its eigenvalues, the fluctuations with measurement error 
will be slower and smaller in magnitude than the actual measured values. 

For instance, a measurement e''ror of one ft/sec will result in a .4 ft/sec 
error for a system with a time constant of .05 sec over a .05 sec time in- 
terval. This property of observers should be considered when choosing ob- 
server eigenvalues. The faster tne eigenvalues the faster the estimated 
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states converge on measured states, and the less filtering is accomplish- 
ed; slower eigenvalues give more filtering and slower convergence. The un- 
measured side velocity estimate is also shown to follow the real state; 
even at its worst it is off by only one ft/sec. 

If filtering is of no importance, a second type of dynamic observer 
can be more useful. The full observer discussed here is redundant in that 
it recreates states which are already measured. The partial observer, 
discussed in the next section, eliminates tin's redundancy by creating esti- 
mates of only the unmeasurable states. The partial observer also uses much 
less computer time. 



TABLE 8 


THE OBSERVER CONTROL MATRICES FOR THE LINEARIZED DISCRETE TIME SYSTEMS 

CASE 1 


K MATRIX 


0.650 

0.050 

0.000 

0.000 

0.000 

0.000 

-0.001 

0.000 

0.600 

0.000 

0.000 

0.000 

0.000 

-0.025 

-1.600 

-0.040 

0.550 

-0.002 

0.000 

0.002 

0.000 

0.000 

0.000 

0.000 

0.500 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.450 

0.049 

0.000 

0.000 

0.001 

0.000 

0.000 

0.000 

0.400 

-0.001 

0.000 

0.016 

0.000 

0.000 

0.000 

-0.001 

1.000 

0.000 

0.000 

0.000 

0.000 

1.590 

0.039 

-2000.000 


CASE 2 


K MATRIX 


0.650 

0.050 

0.000 

0.000 

0.000 

0.000 

-0.001 

0.000 

0.600 

0.000 

-0.001 

0.000 

0.000 

-0.025 

-1.600 

-0.040 

0.550 

0.000 

0.000 

0.000 

0.000 

0.001 

0.990 

-0.001 

0.500 

0.000 

0.000 

-0.012 

0.000 

0.000 

0.000 

0.000 

0.450 

0.049 

0.000 

0.000 

0.001 

0.000 

0.000 

0.000 

0.977 

-0.001 

0.000 

0.016 

0.000 

0.000 

0.000 

-0*001 

0.620 

0.000 

-0.008 

0.000 

0.000 

1 .590 

-1000.000 

-0.993 



DELAY 





















FIGURE 5. A CLOSED LOOP FULL OBSERVER SIMULATION OF FORWARD VELOCITY WITH MEASUREMENT ERROR 
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CHAPTER 7 

THE PARTIAL OBSERVER 


This method for the development of the partial observer is explained 
by Luenberger [4], 

Again, assume the linear system given by equations (7) and (8) which 
is completely controllable and observable. The system equations, possibly, 
by a coordinate transformation, must be transformed so -that 
C - m 0] 

where ^ is the identity matrix. In this case the system is especially 
amenable to this method because it is essentially already in this form. 

The transformation is done to provide a direct correspondence between 
the output and the transformed state variables. The new state variables 
are partitioned 


x = 



where x is the vector of measurable state variables (by the definition of 
£, x is the same as in (8 )) and w is the vector of unmeasurable state vari- 
ables. 

The state equations can then be written as 

JI k Cl 1-34) 

The partial observer is defined from the above as 

v, • * V* + xa * ■ Vk - • V* 3 - {,I " 35) 

Notice that the term in brackets is a constant and is defined equal to 
zero for the correct values of w. Pearranging the partial observer equa- 




r<i> 

*=yy 

$ i 
=yw 


'Z 


*< 

i 


k+1 

<i> 

•p 


w 

k 



tion 
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^k+1 ^=^ww + =^wy^4; + + =^k+l " ^ww^k ~ =y^k^ (II-36) 

it becomes apparant that the eigenvalues of the observer are the eigen- 
values of 


($ - [_$ ) 
v =ww ^=yw ; 

and all other terms are inputs into the partial observer equations. To 
eliminate the need to calculate » assume the transformation 
2^ = w - Ly 

The final form of the observer is 

4+1 ' <iww- J-%vA + ( 4v - Wyw>i4 + %y ' Wyy>4 + ' My> ' 

(11-37) 

A block diagram of the system is found in Figure 7. 

The partial observer control matrix, is chosen by any desirable 
method to place the eigenvalues. In the X-14B VTOL only the side velocity 
is unmeasurable and thus the partial observer is a first order system. 

For Case 1 , 


L = [ 0 0 0 0 0 2000 -2000] . 


(\Pdis 


This placed the partial observer discrete time eigenvalue at 

0.725 

or a time constant =0.18 seconds. 

For Case 2, the partial observer gam vector was chosen as 
L = [ 0 0 0 0 0 -1000 0] . 


The partial observer eigenvalue was approximately .769. As seen 
through these examples the observer gains can be chosen for any desired 
convergence rate. Observer gains are not limited by control power con- 
straints. 

Figures 8 and 9 compare the estimated and actual side velocities, 
when the initial conditions of the system are u = 10 ft/sec, v = 4 ft/sec, 
and all others are set to zero. The states are again relative to the 
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linearized reference states. The system controllers are the same as used 
for the full observer simulations. Convergence is rapid and is 
essentially completed in both cases in less 'than one second. 




FIGURE 7. A PARTIAL OBSERVER SYSTEM WITH FEEDBACK CONTROL 


'SOL 
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CHAPTER 8 

APPLICATION TO FULL NONLINEAR SYSTEM 

The linear full observer is next applied to the nonlinear model de- 
veloped by Roesener [2]. 

The purpose is to determine the usefulness of linear observers in 
estimating the side velocity for the actual craft. Trim parameters for 
the craft were held to the nearest linearized reference values. 

To test the usefulness of the linearized equations for Case 1, the 
nonlinear equations were initialized at u = 5 ft/sec, v = 2 ft/sec, 

0=q = w = p= <j)=r=O. The essential results can 'be found in Figures 
10 through 12. Figure 10 shows the observer estimate of forward velocity 
converging on the actual forward velocity within one-half second, but m 
Figure 11, yaw rate is shown undergoing small oscillations at approximately . 
six cycles per second. The oscillations can be attributed to computer er- 
ror as the oscillations are small (=10 radians/sec) and difficulties were 
encountered with balancing accumulation error and truncation error. The 
linearized system's eigenvalues predicted no such oscillation. The ob- 
server is tracking the yaw rate very well, but the necessary .05 second 
time delay becomes very noticeable in an oscillation this fast. The yaw 
rate in the linearized versions showed no such oscillations. Remembering 
that yaw rate was the only significant coupling between the system and 
side velocity for this linearization, Figure 12 is no surprise when it 
shows no reasonable estimate of the actual side velocity. If computer 
error this small can degrade the estimate of side velocity this severely, 
it can safely be assumed that no accurate estimate will be developed for 
this operating condition on the actual craft. 
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For Case 2 the results were more useful. The yaw rate still oscil- 
lated as before, but for this linearization roll rate was the the major 
coupling between side velocity and the system. In Figure 13, for initial 
conditions of u = 20 ft/sec and v = 2.0 ft/sec, the observer gives a very 
good estimate of side velocity. At u = 25 ft/sec, v = 10 ft/sec (Figures 
14 and 15) the observer forward velocity converges rapidly to the measured 
forward velocity. The observer side velocity rises rapidly to 6 ft/sec, 
and then waits until the real system side velocity slows before continuing 
to follow. Given the initial conditions u = 25 ft/sec, v = 5 ft/sec, in 
Figure 16, the observer estimate rapidly overshot the system side velocity 
by about 20 percent and followed increasingly better as side velocity went 
through zero. The case of u = 15 ft/sec, v = 5 ft/sec is found' in Figure 
17. The observer again waited, this time at around 2 ft/sec, for the 
system side velocity to slow. 



SYSTEM 


OBSERVER 


INITIAL CONDITIONS 


u = 5ft/sec 
v = 2ft/sec 
OTHER STATES =0.0 
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0.00 0.10 0.20 0.30 0.40 0.50 0. 

TIME (fee) 

FIGURE 11. THE OPEN LOOP NONLINEAR SYSTEM WITH FULL OBSERVER 
ESTIMATE OF ROLL RATE TOR CASE I ,u-5,v=2 
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0.00 0.10 C. 20> 0.30 0.40 0.50 0.60 


TIME (sec) 

FIGURE- 12. [HE OPEN LOOP KONLINL'AR SYSTEM WITH FULL OBSERVER 
ESTIIIAIL OF SIDE VELOCITY FOR CASE l,u=5,v.=2 



3.0 


INITIAL CONDITIONS 


u = 20. Oft/sec 



FIGURE 13. THE OPEN LOOP NONLINEAR SYSTEM WITH FULL OBSERVER ESTIMATE 
OF SIDE VELOCITY FOR CASE 2 s u=20,v*2 
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TIME (sec) 


FIGURE 14. THE OPEN LOOP NONLINEAR SYSTEM WITH FULL OBSERVER 
ESTIMATE OP FORWARD VELOCITY FOR CASE 2 ,u=25,V=10 
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FIGURE 16. THE OPEN LOOP NONLINEAR SYSTEM WITH FULL OBSERVER ESTIMATE 
OF SIDE VELOCITY FOR CASE 2 ,u=2S>v=5 




OF SIDE VELOCITY FOR CASE 2 ,u=15,v=5 
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CHAPTER 9 
CONCLUSIONS 

Dynamic observers can be used to create an estimate of unmeasurable 
states needed for implementing state variable feedback control policies. 
Although based on linear control theory, dynamic observers can be applied 
to some nonlinear systems, depending upon the nature and degree of the non- 
linearities. 

The method applied in this research can place observer eigenvalues 
arbitrarily. If the output matrix is not in [J. j OJ form, general modal 
pole placement techniques could be used, as in the control section. A 
full dynamic observer should be used when measurement error is large and 
filtering of measured states desired. The only faults of a full observer 
are a slight time lag and an excessive use of computer time. The partial 
observer eliminates time lag and minimizes computer time, but offers 
limited filtering capabilities. The decision between full and partial ob- 
servers must be made by the designer to fit the individual situation. 

The results from the linearization at 20 ft/sec forward velocity in- 
dicate that the side velocity can be estimated for the full nonlinear 
system using linear observer theory. The eifective range of each linear- 
ization must be determined for criteria such as the desired accuracy in 
the estimate, balanced with the computer time and memory needed to imple- 
ment it. For the X-14B VTOL, these ranges seem to extend further above 
the operating points than below. For instance, the linearization at 
20 ft/sec forward velocity did a much belter job of estimating the side 
velocity at u = 25 ft/sec and v = 5 ft/sec chan at u = 15 ft/sec and 
v = 5 ft/sec. The two linearized flight conditions also seem to indicate 
that the 1 inearization ranges become increasingly smaller as flight becomes 
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slower, culminating in the condition at u = .01 ft/sec which seems to have 
no ability to estimate the nonlinear side velocity because of a lack of 
coupling between the state variables. 

It is recommended that a method for estimating side velocity over the 
entire operating range be obtained by generating a matrix of operating con- 
ditions similar to Figure 18. The aircraft, as it changed flight condi- 
tions, would switch from one set of parameters for a linear observer to 
another, continually using the most applicable. 

Figures 13 through 17 illustrate that to make the desired estimate 
more accurate, more operating conditions are needed. The limiting case 
of mapping operating conditions would be to investigate fully adaptive 
observers. These would include parameter varying observer gains to con- 
trol the model to the system. 

Limited input modal control is a useful and relatively simple tech- 
nique for designing control policies. The designer must approach pole 
placement in a logical and systematic manner. It is suggested that one 
eigenvalue be moved at a time while keeping all others stationary to de- 
termine the individual eigenvalues' effects on the gain matrix G. Using 
this method for this system it was determined that some poles could be ar- 
bitrarily placed with little effect on while Gwas highly sensitive to 
others. 

One effect of arbitrary pole placement of which the designer should 
be aWare is derivative action. This is an undesirable initial overshoot 
during control caused by a zero dominating the initial dynamics. To elim- 
inate this problem, the designer must place poles so that they dominate. 

A most informative method for determining complete controllability 
and complete observability is through the modal transformation. It allows 
the designer to determine the degree of the system's input, output coupling 



117 . 


for individual modes. For the X-14B VTOl, the modal input matrix H, 
for Case 1 and Case 2, was studied in detail before actual design. It can 
be seen readily that the magnitudes are larger for Case 1 than for Case 2, 
which suggests that Case 1 might be controlled more quickly. This con- 
clusion was reinforced when control was later implemented and found to be 
slower for Case 2. The inertia terms in Case 2, which create a more 
strongly coupled system, also require more control power. Examining the 
output matrix, JF, the designer is also able to determine that side velo- 
city will be much more difficult to estimate for Case 1 than for Case 2 
as the maximum coupling is an order of magnitude smaller. This observa- 
tion was proven true when the linear observer was applied to the nonlinear 
system. 

All simulations of systems to be performed on a digital computer 
should be transformed into discrete time difference equations. The trans- 
formation is simple and easily performed through many techniques, just one 
of which is explained here. The method used here is again most applicable 
if modal control is to be used. Discrete time equations have been shown to 
be approximately an order of magnitude faster to execute on a digital com- 
puter than integration routines [3], If a digital computer is to be used 
to implement the control at a specified time interval, it is an additional 
reason to use difference equations. The application of control techniques 
in this research proved to be no different than for continuous time con- 
trol except for the interpretation of eigenvalues. Eigenvectors remain 
invariant. 
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7.“r: (1 ,i>) = rr = ! 1,6) 

T“.T(l,a)=r > !\( i ,<f) 

Hi: illLLOi 


EE-tAmx m l 


mui i,E)=ima, 3 ) 

COMIME 
/VHlTE<b,32) 


22 ZCPr BT ( 1 HO , SX,. .THE P OPA L SIHILARITY TPPXSFOPK TIUTRIX TWftlM 

” 3 1 T £ ( a, 27 ) ( f THAT (I , J ) . J-l, 8 ) , I - 17¥1 

DC Ml I; 1 ,P 

DO h J J- 


1,3 

_T a 7 : (i , j ) = T 5 .LILL*jLL 


80 COMi:ii)£ 

_£ ?.L L_ K I H V ( T r R TI ,8,D,l.m 


M. 


HKI 7 E ( b , 3 j I 

IFg^rP BilMK i l 

AtiUt ( 6 , JiH( 1KMH 

^csrAttmtJS] 


ro 

GO 


CALL 'PHD 1 TSAU, BMAT, HP AT . B,8,0,0i6) 

*f:U(f J 2U j. 

24 F04 V AT( 1H0,5X, * IHPv/T KATrIX HUAI' > 

'-auiz M«i, m 

call «‘>?D(rnAri # pki,phii.8,6,u,o,6) 

CALL KFRCtPH II. THAT. PHIH . 8 .8, Q ,0, 8 ) 

i8ii£(S,2S> 


*■ vAi«nAi ± — J1 J j. — mj.wnr.1 imum_JUAn.. rn i n i 

( 6 ,27 ) ( ( °Hir f I,J),J*1,6),I*1,8) 

DO SO 1*1.6 

DO 90 J* 1 , £ , 

- PH) Mt.(T . J 1 =PrlIM(I, J I 

Pill I L( I.J)*FHID(I,J) 
IMJ-?) 5,90,00 

S 

hL'U ,J)=HMAT(I, J) 
HI..1T ( I. J > = _HL2_M. J.) 


90 

CCNTISJE 
. JRlT£.{5,3n 


34 

FOSSA! fliiU.SX.'THF P ARTIQNEO OUTPUT HRTRIX HLfi* ) 

-"“EXE ( Sj 29 ) ( ( MLfl ( I . J ! , J*1 .6 ) ,1*1,6) 
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'«FITEC£,3 S> 

-15 T.0 ".MAT! 1 ri 0 i:^2^:cLJJI_tt£_^PXLC>i:n.r!_ 


BF IT t <£ ,29) ( (hlKI (I,J).J*1,6),I*1,6) 
_CA LL- .12 u 5 < FJHXCi. xEH 




i*,<tTEtt,40) 

£2 - r IT LIJJ, 5X*JL^E-JMA££iU^i m TY I S-JSILbaJL-lflSfL. ^ 


*’17 E <6,27 H (AKPATf 1,3) .J*l ,8 1 , 1*1 , S ) 

C AIL 1 » ? D ( H LUT ,A<^ AT.EA«Ar.5,6,0.0.3) 

** IT EC c#4i) 

-12 FC O'B : ( IH i, 5V,_* jtr HL*1I TIME S RKfiAT 


132 


r:(b.27i ( a'r?7< i,o>,j*i,s>,;=i,6> 

S&I.L-2S? j ( fzxp- r. ,;."j i; i.c?iAr.6.A.o,o.9> 


KrHATUFlD.O) 

_'£AEf5il321_:r.Artl,UL 

RtAI (5,132) CMAT 1 1 ,2 ) 
-R r A0( 5,132) CHA:U,3L 
Hc/r(b,i:>2) cri,m,4) 
Jil 1 H 4,-U 2 LMAILUS1. 

LEAKS, 122) Cfl-TU,-:) 
READ( 5,132) T“AT< 1,7) 
PEAKS, 132) CMATU.tf 
REAL (5, 132) THAU L 
RFAK 5,132) Of A? ( 5 , £ ) 


’« h 1 7f ( £ , 42 ) 

-42 FOrKAIf UJ,5X,: THF_CPH.TPOL “ATRJ K Vi ORIGINAL FORM 

'■» P 17 £ f f ,27 > ( f CHAT ( I » J ) , J* 1 ,6 ) , 1*1 , 6 ) 

CALL IP HD ( 3X AT , - 5.-J I , OK At, 6 , 5,0 .0.61 


■ CHAD 


CALL r’SU3UHI,0>lAT,CO!»,8,6,0,0) 


FCfr AI(1HJ,5A, ’THE CLOSED LOOP SYSTEM MATRIX 

_XRII£,{ 5. 27J CHOtf- . U,J) . I«JL,A1. 


CON' ) 


CALL EIGE.4P< 0,8, CC*.. 39, ROR.ROC, VECR, VECI , IND1C) 


124 . 



Ill 

ITT 

In 

1U. 


mutG.ui) — 

FO^RTU IJ,5<,*THE PEAL HOOTS OF THE CLOSED LOOP SYSTEM 1 

i r. LI LQj 1 LU_ri25 — 


hC ’. V AT<6F14.S> 

VFIT c( C * 121 ), 


FC^ATl 1H J,5X, 'THE 

/PI ZZSS ,22LUi£!L_ 


COMPLEX ROOTS OF THE CLOSED LOOP SYSTEM* ) 


,/HITsVS , lal 1 
fC^nnvj.SX.’TPF INDICATOR*) 


• fi tTFfft .221 ) INDIC 


IIP CCNIIHl'F 
CALL SEX 


%?- 


SUipOMINE SEX 

DI M F'lSIOH XF[8>,XFl(3),K«At(8^G) f CPMAT{6.ei«,PXlLa J il 

1, f -AT( r , « » . vKUf 6 I .FCMA’i ( 6 / b ) , LGR( S ) , ROC ( (i ) 

2 . b“k T (a.0nOM,\r(0.i*I.VF.CKI8,e).YECIi3,8l 

3.IMICI8) 

fear?'- jiiLjMitm 

DOUBLE PRECISION XK , XKl , KMAT .C PNA I , FflAXt XKO , KCfiAT , GMAT , BOAT, PHI 

R^ADfS,! )i J.f 

10 (OSXATimO.O ) 

~ - Q . 0 

L = 9 v 

CC’.Os 

CALL NFREfBNRI.GXAT.KCMAT.'d.S.O.Q.S} 

C 1 LL_nS OBJ PHI , XCRST. FMA~,3 ,g,C,Ql 

, 30 20 1=1,500 

_ i . ~ . <■ I .... 

Cn_L .IrRDf F” AT, XK ,XX1 , 3 , 8 , 0 ,0 ,1 ) 

IKI .JIE. 10)_0AJ.Q_Z5 

nHIEf 8,23) T 

_? \ fO-RAT! F 1 0 . 5 ' . 

Lk;U<5,27! XX 

l o : 

21 FORMAT! 3£9. 3) 

18 00 n J =1 ,3 

19 XK<J) = XK1(J1 

I = T*TC 

20 CONTINUE 

STOP 

END 


125 . 



FULL OBSERVER SIMULATION 


.DIMENSION XK<9)*KHAr<8».Z).Hll<?,8.}iCPMATma>»P*! 

l*GMAT<6*8>,KCHAT<8,6),AxVa>*Pm<8,e>^ 


_2 

3*XR( 8} 

_RE*U KNAT*KCHAL*Ut_ 

RE*0< 5*9 ) XK 


10 


forma rear lo.oj ” — 

.D(J 997 1 JK*1 * 8 

DO 998 J«»l>7 


— KHATtl JK*JK)« PHICtJi<«.llO 


99)* 


CONTINUE 
_KKAT(!*I >..63 
KKATC2*2>».6 


_KH AT < 3 > 3 )•» , 55 
KMA 7 ( 4 * 0 )., 50 
-XHAU5*5)a.93_ 


KMA7(6*6)«,40 
. KHA T C7 * 7>«1 ,0. 


KNAT(8»/)**2000» 

ftEA0C5*n> CPHA1. 

11 F0R'(AT{7F10,0) 


■RC : AD( 5 >m((fSia<J>*.)M li.& 3 jlPl.n> 


11 FORMAT { 3E 1 4 , 7 } 

nEA0C5M0). GHAt 

T»0,0 

TC**. 05 


CALL HP«0(XHA7*C?FAT,KCHju»3*77o7oT8) 
XALL,KPROCPSI*CLMAIjrmO*e* 6 f 0 ,Q.A> 


X = 1 

XK1(3)«10, 

00 25 I « 1 . 200 ~ 

__ oo »i x-t.r 

9 J 4 K(K)«AK{K)«RAN 00 M(X)*.‘ 2 *XK(k)-, l.xktX) ' 

XK<3)»xmc»»iUNDIHm«2.«l. n 

CALL RPR0(PHI*XK,AXK.fl*«*0.0»l> 

— call HPRO(PSIG*XXO*1N*A*A*0*0*11 

CALL MSUR(AXK*IN*XXl*8»l*0*0) 

CALL HPRP(PHI*XK0*XK11* 9 *3*0*0* 1 > 

CALL HSUA<XK11*1N>XMP,8*1*0*0} 

£ALL-ti5UOtXilOjXK/X.1E*a> 1.0.0> 

/* A t I UOnR / UAu i . u..m ^ . 


CALL KPRD(XCMAT» XKE* XKEC *6«A*0*0* 1 ) 


ro 

<r> 


3 1 roau*T<no.5) 

MIIIEU* J3>. *R 

KR[T£C6,33) XK 

«RITC<»f 33 } JtXCI_ 

JJ fORHAKSF?, J> 

*a:TE(«,j}}_xKC_ 

T»T*TC 

00 .25_J«J,JS 

XR ( J 5 »XK 1 ( J) 

* K( J)*XKt (J? 
XK0{J>«XI<l0f J) 

Zi CONTINUE 1 

STOP 





PARTIAL OBSERVER SIMULATION 


-JHyEV&Xa -P S-LCA aS It U • C H A UT t B Xt 

1 rR<r)#*L*tt»n>8LflU>6)#CM*T<4/8),(.MAT(!,r)f*XK<8)<INfiCe) 
■2»ot)r(aj.*wrti»?3«AY?<r.»yjL^AyYt.ci«-rjjaj< ( wbtu i»-8. u . 

JAYHCr >»AYt(LC \ >»YU<l>#YAlAf t >/!NCn 

_3S;aL G4AT* LHAtAOUljulilfitiJt-AlSJ — - — — 

RLAOtS'O) XX 


9 rcemTC rr 10 , 0 ) 

P.CA3C51 12 ) .GHAT — 

8 FOR«AT( AT t 0, 0 5 

RE AD< 5> 9 3 C U AT 

BEAMSMOJUPSl (1* J)»J*5#63»1»1*8> 

-to ■ ..rORHAUJUA^JJ — — 

R£A0(5*t2><< PHKI»J)*Jol/n)#Iol»fl) 

A2 FDRvmaut.73 — 

RC*0(5«9)L W AT 

DO 6? J-J.7 — 

oa-62-i«i*i 

42 AlMfO»J)ePKI( I.J) 

CAUL- PPRQClMAr»AY.Y-*7O'YLfi.*?<0.».Qj.7J 

CALL HSURtAxY.AYYUi At.A»*/7»0»0) 

PP 6 3. - J 3 ! 4 — — - — - — i — ■ — — ~ 

BK< I > J ) «P S 1 ( <1 » J > 

.00 .4-3 I»i >1-. — — 

43 eru< j)»ps!C i*jj 

CALL HPRDUHAT»8Y<BYL#.i>J’.»0»0*Sl 

CALL H$!l9CP>f»8YL»BL8#l»6#0»0) 

’ 00 6A !»1»7 

6* ay*c n*pHia»o) 

A<iW*PHUi*8J — — 

CALL wPR0<L«AT»AYH»AY>‘La1»7»0aI>#1) 

AYVLS »AYWL< 1 3 = 

iLC*A*«*AY*Ll 

T»0,00 - 

TC-.05 

L » 0 ■ - — 

27 rDR*»7t7rl«.5> 

-27- — rcaMiKariA.5 3 

28 F0R M AT<*ri»,5) 

X«l- 

00 24 I«1*JOO 

OO-Sl-Xat.* 7 — — 

9 1 YRCK J.Ya(K)+RANOO«(X)*.2*Yft<K>*.I«YR(K3 

CALL P?RD<CHAT»XK»YR/7*8>C>*<i»t) 

CALL HP«0CPH:,XX*AXK»8>8»0»0f.iJ 

CALL MPR0<PSI»U#INR,4,4f0,0>l) 

QM ■_ MXlJ P<;AYY..tHfi .-YYl 
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YRlltO.O 

YALASsQ,0 

oo ?o J=!»r 

YA!.A\t>YALAU4LAC 

'•RL1*Y8U*LH»T(1# jmsco) 

-ao-cc^u we 

ZYRl = Z‘Y?!.1 

ZClOSOsALCiZYRi 

Zl*rAL-*l-JNUZCL03£ 

00 61 Jo 1 < 7 

61 0UT<J!nYR{J5 


— 0 O-U-6- 1 a Z ‘OU. . 

CALL M?RD(GhAT*CUTfUj&fO»0»O»i) 



00 101 J»1 # 6 

tU1nTWt»ftl ail. tHIIM.I 1 . . 

101 

CONTINUE 

rcanflE**/ fl.i-JY.fti 

in l .mc, i > go to oo 

ii 

rnRKitt rio.si 

KRITec6,>3>XK 

- u a 1 1 c > 7 

31 

fQPMATUEO.j) 

80 

'continue 

T*7iYfi 


z»n 


—oo as. j«i*o . 


XK { J ) »XK 1 { J) 

U CONTINUE - 

STOP 



SIMULATION OF NONLINEAR SYSTEM WITH FULL OBSERVER 


.REAL MA5S.tX.tY. lZt I XZ. LI1HE. IENOjLLDXIX.. LPIXaMCiYLIj 

l M03 Y.NOZU.NRtZ.N'JNO I M, !*'<*<>* 

REAL K u ATSta»n,K ,i ATCa.7).PHIS<a,ai*PJlICa.a)jCPMATC.7ja) 

l,PSIS(5.5)»P$rr<ft.8).XVJ<8).6M\TSC8*ft>*<>MATFC6*8).KCMAT$<8»8>* 
Z-KC’.AUa. 33*PSl-lSli'..fl.l-.P.SI.OE(-fl.A).TN(flI.X1*lprn>.Y Kf(fn.ymO(fll. 

jXXtt(0>.Xxec(6> 

4.YSH0RT(8). 1-IOf 8)»YRCrC35 

COVuCN/E$OfD>(/H,Xl,r I ( 99 ) * X * Y { 99 > . TR, H» KSTp , DYC 99) 

CC^PN/f 'JNCI/0»IHRUSt/RAS$»Si5HA>LA‘<fl0Ax 

!UO.XAREA,CL°O.DFLCL.ALPHA.rO.VKEL.lKmAH,CY*TX, IZ.m.IY, 

ZZIha-iLOSE. ZGNE.LP lX.LDXIX.nELXi.SPAH..NREL* CLQ. CU3ELA. 

3 r>E l A»CL B FM»IENO.n“E , >A«XTn(l*xO'(£»MR ! Z. MD Z 1 7. OEL It CNO. 

ACNnELR.&ELR.CNOELA.MOYlY.OELY.MJlIY.CHORD.CMOPO. 

50Elf “0*C«ALP0*OCLCMA, CMQ*C*")f LE . DELE* 3X Y A, J> UREL » V JET* RHO 

CO wu nN/CONT/I.N .. . 

CAT A G.wAS$*VjCT.RM(l,f(ARfA»$pAN.CHOfin»TX,JY»IZ.lXZ*ZONE*IENG*XONE, 

lUXlX/»3YIV>PgiY>0£LertA/CRCELEjH0Zl.Z.i.MRlZ.r-Wg/3?.17A. 13(1 .3537* 

2l6l 3.»o.no237a» 162.49. 3J.bY.5.58»2340.*3400«*5400.>180./0.583*O.S» 

36, 0*0. CAP 1 *0. 035o.-*O. t5* 0,5ZA/_*0. Q1ZB.*J..(HAA»-C.20 a* , -1.L..&X 

OAT* lPTX.A^5VEL/-0,45.t72fl./ 

DATA XTn(l>7-rWn.l-f;HC/.-»e, 01 33. 0 . 91 5 * 0 . 1 47/ 

UJ*1 C.O 

— la-tORyAtr.5nc.-a. 1,53 „ 

n Foppf t < ar i z. i ) 

h« .cl — 

RE AC (5*2001 >({PHI(I.K}.Knl.B)«I#</0) • 

2O0W0RMAT(4E}4,7) 

REACCi.Pooij-tjauu 

2003 f ORMATt 7F 10.0) 


N«8 

XIt.O.0 

XSTPtO 

xonf* t 

CC IS I?IY*1.3 

fFAn(5*ioz) r i c i > 

REio(5»t ozs-ma) 

YC3)»TI( 3) 

4-03 rORKATCFAO^CW 

M8)*YI(») 

00 15 J« 3 * 2000, 

URCLtYd) 

,‘Rfl«Y(f ) 

KRFl-YU) 

i:Cflia = 5CPTCi;2CLAA3*A!REL«A2yt«EXJLt2A- 

IIO»S9RT(>if(EL*»2+wnEL*«2) 

ALPHA*ARSTS(tlRfL/ua) 

BETA.ARSINCVKELZUORAR) 

iFaio.LT.siO-CL^osc.cita.aog^-iiiAi.’o 

If (U0.0t.51 . ) CL B 0*U. 75 

CX L CL «? .J AJip.JXo 95 ,llQ.<.;>-2.24 *39aX~07«UO«*3«! .3321 49r-03»U0 

C 0*0. 11*1./ <2. 1529 35* JO* 5 2,5) 

tnH$(P r I4;_,iT,q, !89Z^CY»-j ,35 l.pfiETA 

If (AR5(8;tA ),OE.0. 3 89.ANn.A85t BETAJ.LT. 0,4lSJCYn(”0.33A0,39Jtf 
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, It lT>5(fcCl& l-o, l B<sm( nnstlCT* VUST.il — — : 

!* (*«i 5 <nt t« » . nr.o« 4 ifl. t-vn.iusCrfSTU.Lr** » 2 ?> cYn{»o»ut-(j.36?* 

TC41ii.tKt.J4 J-.3,4I8 3 2*CA«StPET4 J-CaLTA.,3 . 

Tff *RSt«rT»J.RF.l .?«1 CT»«<i.55*< iS'S tlSET a JPUETjO 

in J?sf nrr* j . t. t. t>.. i «■?.'. eta : „ — _ 

;rr3?^(fif ta >, 5r,fl,!.i?.AVfl,M»sti!crA) .tT.o.sru csr>*i-n*oi4*-«>.oo<!«* 

U AR5T MtT I >-B,I.17Jl»(AR^(Fi£T4>/RETA) __ , _ 

IF< AuA(Ht.t*l,qr.a.4Tl.AHtl,aai5<3ETAJ*LT.lf J’Sl CT.O*t -5. 0 1 Z-fl . nOaaS* 
L4A3 Af PtTil-{l»4r! J.)*C.AR5tUE,l43./liCY.a.3. - 

Trfii>4(f!rr4 j.qr.s , tjsi Ci ABstafTA j/stta > 

CtOCl 4*0. S03*COS(«EIA l*4 2^~_„ .. . - — 

qL t tTj£*f!,o')A3*ci I ! , n*t',ci?? 

irtLc^Lf, }«, i -CROras'S.i5*Qtoq27lT(jn . . „ ._ 

If f 'JfliGr* o.«*,C.u*j.L‘i<S t » > CHnr»0<r“O t 5?*0»(H *S*TUf) , 0 ‘W 1 

ijr( un.at,si.L,,ivQ^ijr,,!. :-.uuo-J^3PS**4^2£».3^a/jiiriu , . 

Huf'-qi-i 

:ai,5,c’'3»o,*a,i5 ; 

DCs C'*P»«.Y'n35S'0?*(;n**3-3. H#'»?E-q5*uQ**J4‘l.52S505S-O3«U0 

tr(ufv.j.T.3»,j c^iLRo-A . ua.a»i4s*ua ~ 

rrfi!n*r.£. 10 . 1 . r .si , 3 ch/u_pow.. 7 *Oi. igi>r*uia«i 4. 3 

___ 

J"f ST < 5 *" T 47 .1 r. 0*17053 e!(floO*!n*flETA 

iTSAesifltTA.T.nr^Q.uss.ftNa.Auitijtza-uui-.a.isaj-csiiBaHULA.aiestsE.uu- 

J /P£TA ) 

, -ir< 6js5tf!Ef4J*r.£.a.57 4. ANf).AnStft£JAi.l r X.iv2.9iJ_*4*ia*LC0-Xlii3Aa*l-i*At*~ 

t C A»5< pfTJ ) -O.S'M 3 3 •{ AaSiUET JJ/SerA> 

~_Jr '. arstre M.i.ss, u_;7£3-.C- i uj.-;CuJ-UaX.i n s.;.f m ; /beta 3 __ _ . 

cvo>?i p*o,f"3iiB*cus(getA)*4? 

EMCt n-O.DrtOi 4*C0S(3E.TA4-AAJ 

*1 e*t c f S'liitTcxnoc.t j 

V DP E 1 3 — ■ 

PfUT*' Jftr<7}/7<>‘ 

. TAtUi 2 l>H J1/3B. -. ■ . ■ 

.pki 1 3 . r )»t t a ) /? a. . - „ 

P>Mt*.#/3t-»t33/20. 

P«HS. J34-Y(7 3/20.- „ 

l>K[<S>J3*T->'(C/{30,*y j£tl 

. s>xUl*4 3 5-T£I>JJll^. 

»rf| (<W0 )I £1Vf/(?0.«V)CT J 

PAlt fl*U3*’E>t»UiT/C2.9,M<isS«iUtJ7 

pitjf T'-’n-T'i'je/c?o < *Y 3 Erj 

iC«v»itatCilO >Tj-T . .. ...... 

00 Jo S ■ ! . * 

~i.» tr )»Y **>- . . — s „ ., - - - ■ 

oa to5 Loi>t> 

.03 lo5 a n \ , r ™ — — 

105 ■(Htru -■0»P>U<Lf>4> 

KV4TT 1 * 1 34 , £5 . , .. , ■ 

T«'»T <?*?>*. A 

~~*>1AT<J»13*.S5— 

**i rjA'i i«'.5i> 

^ slJtS.ilP.OI — ™». — ., — 

4HtrfA/«2«t.?T7 

™ XAT( 0 -f.)a-li ,;fv v t 


c«u M I, fiC{>!PiT»CP , <AT^i;‘<Ar»aAr>o^oAO> 

„C!iSL.'3Pptl£P4>I,S!10*ii',lP»S^»*a».0*.1.3 

C.V,S MS05CTK.0pT5U3'"Tf T r .E»Oo!*<fa0 3 


net 


CALL KPRCtKC“ATjXl'£*JCIt£Oaj&AlUOjaJ 

CALL KSlirc**lP*XKrC»* , <tD*«»l»0tO) 

ircxiticti ).L£.^.Ln..ir,atfl.-ijaa 

lOt DUEn-t.? 1 ' 

0£t Y«*l,90 — 

sinwA».«r/5r. j 

DUE GAn9 *.76*1728. <■ 

DCl. />•<>•') 

TWRl!SI*?.MAA*945W*«a^iv=Am^JjAJJJU-aJ5^ 

iruCP'JNT .LT, 5)00 TP 15 

GO in 4001 _____ 

too D£l£».«7 

OCLY-.U — - — 

SlO"**n,o 

GUEL*aJ«*5.HH«A- 

TH»UST*?.«(41.9455?«9S.17«*iaA.«8J)**a0A 

if <! count . .LY,_3)fto-ra-i3 — 

<500 1 Hf J !Tt{4*? f )05)X»('f(t>»I»t>8)*TR*XK0/XXE 


2005 FOR , '»T<tT|M£ !$*»F10.5./ f 'NYSTE* « **<>EU,3»/ 

t»-tOR4£?^L«*JUftE44-^W-*JXfiAlW — n J . , a £ 1 1 . 3 » f/i 

ICOLNT.O 

00 1003 XuljQ- 

XNO<* )°X<tO{X> 

AOQJ CONTINUE 

15 CONTINJE 

JX0Otl~lTDS — — — ■ — - ■■ — 

E»0 


.FUNCTION FC1 »X.Y) 

COMKON/fSO/UOOP 
.HEAL_YC991»IJIlfiJ*HASS.LAVat>Aj-IX..IZjaX2.*XiM.DNE.l 


1< IEf , G»NPt/fN07J2»NnYlY*M0tV 
_jCO ••.‘inH/JLO 1 T / IH _ 


C0u|.on/f-|i'Jrl/P.*ThRUST><'ASS>5T6'4A,LAM!|0A» 
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III. Wind Gust Analysis of the X-14B VTOL Aircraft 
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Wind Gust Analysis of the 
X-14B VTOL Aircraft 


NOMENCLATURE 


Variable or 


Constant 

Description 

Value-umts 

u 

forward velocity of the X-14 

fps 

V 

side velocity of the X-14 

fps 

w 

vertical velocity of the X-14 

fps 

P 

roll rate of the X-14 

rad/sec 

q 

pitch rate of the X-14 

rad/sec 

r 

yaw rate of the X-14 

rad/sec 

X 

longitudinal displacement of the X-14 

S 

feet 

y 

lateral displacement of the X-14 

feet 

z 

vertical displacement of the X-14 

feet 

0 

pitch attitude of the X-14 

rad 


roll attitude of the X-14 

rad 

* 

yaw attitude of the X-14 

rad 

^net 

net thrust 

pounds 

a 

diverter vane angle 

rad 

A 

side vane angle 

rad 

6 x 

roll reaction control nozzle angle 

degrees 

6 y 

pitch reaction control nozzle angle 

degrees 

6 z 

yaw reaction control nozzle angle 

degrees 


aileron control angle 

degrees 

*e 

elevator control angle 

degrees 


rudder control angle 

degrees 
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Variable or 


Constant 

Description 

Value-units 

si 

engine angular velocity 

rad/sec 

UW 

steady-state longitudinal wind 

fps 

VW 

steady-state lateral wind 

fps 

WW 

steady-state vertical wind 

fps 

UG 

longitudinal qust component 

fps 

VG 

lateral gust component 

fps 

WG 

vertical gust component 

fps 

g 

acceleration due to gravity 

32.174 fps 2 

in 

mass of the aircraft 

130.35 slugs 

! x 

moment of inertia of the aircraft with 
respect to the x-axis 

2340 slug-ft 2 

r y 

moment of inertia of the aircraft with 
respect to the y-axis 

3400 slug-ft 2 

i 

z 

moment of inertia of the aircraft with 
respect to the z-axis 

5400 slug-ft 2 

'xz 

product of inertia of the aircraft with 
respect to the x and z axes 

180 slug-ft 2 


moment of inertia of one engine with 
respect to the axis of rotation 

0.5 slug-ft 2 

s 

wing area 

182.69 ft 2 

b 

wing span 

33.83 ft 

c 

mean aerodynamic chord 

5.56 ft 

p 

density of the atmosphere at sea level 

0.002378 slug/ft 

X 1 

distance the center of the engine intake 
is forward of the center of gravity 

6.0 ft 

*2 

distance the center of the engine exhaust 
is behind the center of gravity 

-0.0133 ft 

Z 1 

distance the center of the engine intake 
is below the center of gravity 

0.583 ft 

Z 2 

distance the center of the engine exhaust 
is below the center of gravity 

0.916 ft 
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Variable or 

Constant Description Value-u ni ts 

effective length between engine exhausts 0.167 ft 
and diverter vanes 


C Lpo 


net exhaust velocity of the engines 1613 fps 

lift coefficient, power off 

difference in lift coefficient between 
power off and power on 


C Q drag coefficient 

C coefficient of side force due to side 

y velocity 



rolling moment coefficient 

rolling moment coefficient due to aileron deg -1 

deflection 


C^po pitching moment coefficient, power off 

AC Mq difference in pitching moment coefficient 

between power off and power on 



pitching moment coefficient due to angle 
of attack 


Da 

AC 


C Mq 

C MS 

e 


lift coefficient due to angle of attack 

drag coefficient due to angle of attack 

difference in pitching moment coefficient 0.576 rad 1 
due to angle of attack between power off 
and power on 

pitching moment coefficient due to pitch -11.4 
rate 

pitching moment coefficient due to eleva- -0.0178 deg 
tor deflection 



yawing moment coefficient 

yawing moment coefficient due to rudder deg ^ 

deflection 



yawing moment coefficient due to aileron deg ^ 

deflection 



rolling moment coefficient due to side 
vel oci ty 
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Variable or 


Constant 

Description 


Value-units 

VC 

rotational damping in roll 


-0.45 sec 1 

VC 

rotational damping in pitch 


-0.15 sec 1 

V‘z 

rotational damping in yaw 


-0.20 sec" 1 

L « V 

X 

reaction control rolling moment 


0.0681 rad/sec^/deg 

Vy 
y J 

reaction control pitching moment 


0.0350 rad/sec^/deg 

V 1 * 

reaction control yawing moment 


0.0166 rad/sec 2 /deg 

% 

longitudinal diverter efficiency 
based on engine speed 

factor 

0.9139 

% 

longitudinal diverter efficiency 
based on engine thrust 

factor 

0.807 

a 

angle of attack 


rad 

6 

sideslip angle 


rad 
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CHAPTER 1 
INTRODUCTION 


1.1 Air Turbulence. 

Nature confronts man with certain obstacles in his use of the air as 
a medium of transportation. One of those obstacles, turbulent motion of 
the atmosphere, was the subject of this research. An airplane is subject 
to random external forces which result in random variations of attitude 

i 

and trajectory. The time scale and intensity of these responses are gov- 
erned by the scale and intensity of the turbulence, as well as the charac- 
teristics of the aircraft. Their effect is to produce fatigue in both the 
pilot and the vehicle, to produce an uncomfortable ride, and to impair 

i 

precise control along the flight path. 

1.2 Vertical Take-off and Landing Aircraft. 

A vertical take-off and landing (VTOL) aircraft is affected to a 
greater extent by the turbulence of the atmosphere whenever the aircraft 
is in the hover or transition modes of flight than when in normal aero- 
dynamic flight. The aerodynamic forces used for control have negligible 
effect in hover. The only control available to a hovering VTOL is the 
stability augmentation system built into the aircraft. There exist dif- 
ferent control systems for different VTOL configurations. A particular 
control system for a specific aircraft is the subject of this research. 

The airplane studied in this project was NASA's X-14 research VTOL 
aircraft. The control of this particular aircraft is accomplished by 
vectoring the thrust of its two engines (to provide thrust for hover and 
propulsion) and by bleeding air from the engines and ducting it to control 
nozzles in the wing tips and the tail (to provide thrust for attitude 
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control). The control effort available is limited since too much air bled 
from the engines lowers the thrust capabilities. A trade-off between 
thrust capabilities and control effort available is a major factor in de- 
signing VTOL control systems. Therefore, the knowledge of control effort 
needed is vital in designing the VTOL. 

Control effort is defined as the angular acceleration produced by a 
control input. If an angular acceleration produced by an input disturb- 
ance, such as a wind gust, exceeds the control power available, then the 
aircraft is considered unstable. The line between stability and insta- 
bility depends not only on the criteria used to define stability, but also 
upon the characteristics of the aircraft. For the X— 1 4 the stability is 
based on the angular accelerations and the attitude of the’aircraft. 

These criteria were applied to the hover mode where they are far more 
important than for regular aerodynamic flight. 

It is important then to be able to calculate accurately the angular 
accelerations produced by the turbulence of the atmosphere. A model of 
atmospheric disturbances is necessary to perform this task. Modeling the 
atmosphere is difficult using explicit functions of time. Studies of 
VTOLs and other aircraft involving the wind most frequently use a statis- 
tical, probabilistic model for the wind model. 

1.3 Research Objectives. 

The objective of this research was to find the maximum allowable gusts 
which the X-14 could encounter while in hover and still maintain a stable 
flying condition. The general procedure for solving the problem was to 
generate a simulated velocity field from a standard wind model. This field 
was then imposed on the simulated aircraft as a disturbance. Due to the 
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nonlinearities of the aircraft, there were changes in' its aerodynamic char 
acteri sties. These changes resulted in a modified aircraft motion and 
were duly considered. The aircraft control' system sought to drive the air 
craft to its reference flight condition, hover. 

The basic solution was found by arbitrarily selecting a wind model, 
observing the output variables of consequence, and changing the amplitude 
of the wind model until the aircraft failed to return to its nominal 
state. 

The models for the aircraft and wind are presented in Chapters 2 
and 3 respectively. The solution procedure, in more detail, is presented 
in Chapter 4. The results of this study and a discussion and interpreta- 
tion of these results are presented in Chapter 5. 
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CHAPTER 2 
AIRPLANE MODEL 


2.1 X-14 VTOL Airplane. 

The aircraft used for this project Is the NASA X-14 VTOL research 
aircraft stationed at the NASA Ames Research Center, Moffett Field, Cali- 
fornia. The airplane was built by the Bell Aircraft Corporation as a re- 
search vehicle for the study of VTOL flight. 

The X-14 is a thrust- vectored VTOL airplane. Two General Electric 
085-19 jet engines are fixed on the aircraft. The thrust of the engines 
is vectored by diverter vanes at the jet exits. The thrust is diverted 
down for hover, whereas, for the transition of the aircraft from hover 
to normal aerodynamic flight, the thrust is rotated rearward. 

The X-14 is equipped with a stability augmentation system for control 
in the hover and transistion modes of flight. The control effort is ob- 
tained by bleeding air from the compressor stages of the engines. The air 
is ducted to reaction control nozzles located in each wing tip and in the 
tail. The aircraft attitudes, rates, and accelerations are monitored by 
gyros and accelerometers. Various types of controllers had been tested 
and used, but they were not of importance to this study. Some type of 
control system was needed to do the simulations of this project and an 
acceleration command (manual) mode was used. 

The airplane was assumed to be a rigid body and therefore had a body- 
fixed reference frame with the origin located at the mass center (Figure 
1). The aircraft was symmetrical about the x-z plane with the positive 
y-axis pointing out the right wing, the positive z-axis pointing down, and 
the positive x-axis pointing in the direction of forward flight. The 




z 


L = rolling moment 
M = pitching moment 
N = yawing moment 
p = roll rate 
q = pitch rate 
r = yaw rate 

[X, Y, Z] = components of resultant aerodynamic force 
[u, v, w] = components of velocity of 0 relative to the 
atmosphere 


Figure 1. Body-fixed Reference Frame. 
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translationa,l displacements, velocities, accelerations, and the components 
of the resultant aerodynamic forces were considered positive in the same 
sense as the coordinate axes (body axes). 

The airplane also had another reference frame attached to it. The 
vehicle-carried vertical frame (see Chapter 4) had its origin located at 
the vehicle mass center. The angular orientation of the vehicle was then 
the relation of the body axes to this vehicle-carried frame. The relative 
orientation was expressed by the Euler angles 0, <j>, and \p. 

The position angles were considered positive if an observer, located 
at the origin of the reference frame (the vehicle-carried frame in this 
case), looked along the axes and saw a clockwise rotation of the airplane 
axes (body axes). The yaw angle, ip, was positive for a clockwise rotation 
about the z-axis (i.e.» the right wing goes back). The pitch angle, 0, 
was positive for a clockwise rotation about the y-axis (i.e., the nose 
of the airplane goes up). The bank, or roll, angle, <J) , was positive for 
a clockwise rotation about the x-axis (i.e., the right wing goes down). 

The angular velocities and accelerations and the aerodynamic moments were 
considered positive in the same manner as the angular displacements. 

The stability augmentation system reacted to disturbances by creating 
reaction control moments opposite to the induced disturbance moments. The 
reaction control moments were dependent upon the angular displacements of 
the sleeves located in the ends of the air ducts (see Figure 2). For 
zero displacements, the sleeves were oriented so that the air was expelled 
through the top and bottom control nozzles to provide equal amounts of 
thrust. Zero displacement was the reference flight condition. 

The reaction control nozzles in the wing tips corrected any induced 
rolling or yawing moments. Rolling moments were generated by changing 
the difference between the left nozzle exit area and the right nozzle exit 




Figure 2. Definition of +6 
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area. The angular displacement, 6 X , of the sleeve in the control nozzle 
was considered positive if it created a positive rolling moment (Figure 
2). Yawing control moments were created by rotating the thrust vectors 
from the left and right nozzles in opposite directions about the lateral 
axis (y-axis) of the aircraft. Pitching control moments were created by 
a reaction control nozzle located in the tail. Changing the differential 
area between the top and bottom exit areas produced the pitching control 
moment. The reaction control angles for pitch, <5y, and yaw, <5 Z , were 
considered positive if they produced positive pitching and yawing moments 
respectively. 

The aerodynamic controls were directly coupled to the reaction con- 
trols through the pilot controls (the stick and the pedals) as follows: 

aileron control angle, $ a e <$ x , (111-la) 

elevator control angle, S e = and (Ill-lb) 

rudder control angle, 6 r = S z . (III-lc) 

Two angles were associated with the thrust vectoring. The diverter 
angle, o, had a range from zero degrees displacement (thrust vectored 
straight down for hover) to 70 degrees displacement for full aerodynamic 
flight. This angle was measured in the vertical plane counterclockwise 
from the positive z-axis (body axis). The exhaust side vane angle, 
was measured about the x-axis (body axis) by the right hand rule and had 
a range of ±25 degrees. 

2.2 Basic Model. 

The equations of motion for the X-14 were taken from [1]. The model 
in [1] was a ninth-order system with the state variables being u, v, w, 
p, q, r, 8, <J>, and i^. These variables are defined in Appendix 4. For 
this project three other state variables were added which allowed the 
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displacements to be calculated. The model was of the basic form 

x = f(t;x»u,t«) ftlI-2) 

where x was a twelve-dimensional vector of state variables, a was a nine- 
dimensional vector of control parameters, and w was a six-dimensional 
vector of input disturbances. 

The problem was to solve for the transient response of the aircraft 
subject to a disturbance(s) . The hover mode, k =0_» represented the ini- 
tial conditions. The disturbance vector, w, was known for each instant 
of time. The control vector was calculated by setting _x =0, the reference 
flight condition, and solving equation (2) for u at a specific instant of 
time. Equation (2) was then solved for x at the next instant of time. 

This was continued until the nature of the transient, whether a recovery 
or a crash, was determined. 

2.3 Changes in the Original Model. 

The equations of motion of the basic model were reviewed by the 
author, more for a better understanding of the model than to see if they 
were correct. The equations were found to be nearly correct. Only a few 
discrepancies were found in certain terms of the equations. 

The term accounting for the rolling moment due to side velocity was 

omitted from the original model. This term, the stability derivative, C 0 , 

6 

was derived following methods of [2]. The derivation of this term is pre- 
sented in Appendix 1. 

In the interim period between the time the basic equations were first 
written and this project started, the engines of the airplane were replaced 
with newer, more powerful engines. The newer engines, the GE J85-19 
engines, had greater thrust capabilities. With the new engine data incorpo- 
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rated into the model, a simulation for the hover mode was conducted and 
a discrepancy between simulation results and actual flying results was 
discovered. 

In hover the actual airplane was flying with the engines running at 
98 to 99 percent of maximum engine speed. The model calculated the engine 
speed to be only 90 percent. The conclusion drawn from this was the ef- 
ficiency of the diverter vane had been neglected in the model. An effi- 
ciency factor was added to the model to adjust the engine speed to that 
of tile actual aircraft. This factor also had a direct effect on the thrust 
and the exit velocity of the jet exhausts. 

Reference [1] also calculated this factor and reported it to be 0.91. 
This was based on a direct correlation between the engine speeds of the 
model and the actual aircraft. For this project, the engine speed, fl, 
could not be calcualted by any means other than from the thrust, = 
fl(thrust). The thrust was an integral part of the model and was included 
in the control vector u. Therefore, the efficiency factor was re-calcu- 
lated based on the ratio of the gross thrust for hover as calculated by 
the model and the gross thrust for hover of the actual aircraft. The 
factor was calculated to be 0.807. This factor was assumed to be constant 
over the range of operating conditions of the aircraft. The calculation 
of this factor, Cg, is presented in Appendix 2. 

The nature of the solution procedure (see Chapter 4) required the 
modification of two aerodynamic coefficients, AC L , lift interference, and 
AC Hq , pitching moment interference. The corresponding lift interference 
force appeared in the equations for forward and vertical velocity. The 
pitching interference moment appeared in the equation for the pitch rate. 
These terms required that the value of the net thrust be known before it 
had been calculated. Both coefficients were defined as functions of two 
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parameters, the relative velocity between the aircraft and the atmosphere, 

U r and a non-dimensional i zed thrust coefficient, T . The equations for 
o c 

AC^ and AC^ o were 

AC, 7 

— = f(U Q ) = -2.248398x10”' *U^ +2. 483009x1 0~ 5 *U^-1 . 3321 49 x10 -3 *U o (lll- 3 a) 

c 

Q 

—• — = f ] (U 0 ) - 4.900354 x10” 7 »U 3 -3.164620x10“ 5 -Uq+1.529505x 10' 3 *U {lll-3b} 
T c 

'where T = T ^/(T/2 . p . U 2 -S). 
c net w o' 

These aerodynamic terms were incorporated into the thrust terms of 
their respective equations by algebraic manipulation (see Appendix 3). 

This '‘rearrangement" eliminated the necessity of having to arbitrarily 
select an initial value for the thrust. 

The last change involved the aerodynamic coefficients C , C , and 

, f I _ X* 

o o 

(^ . These coefficients were "extended" to be defined for the full range 
of the sideslip angle, 8, +90 degrees to -90 degrees (see Figure 3). The 
extended portion of the curves were based on the known curves. 

2.4 Final Study Model. 

The model finally used for the study was a twelfth-order system. 

The state variables were the three translational velocities (u,v,w), the 
three angular displacements (6, d>» $ and associated angular velocities 
(p,q,r), and the displacements (x,y,z). The full set of equations of the 
system are presented in Appendix 4. 

These equations of the model that define the dynamics of the X-14 were 
nonlinear differential equations. More correctly, equation (2) was written 
as 

x « f{t;x(t),tt(t),«>(t)} . - (lU-4) 

With the control vector assumed to be a known function of state variables, 
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equation (4) then became 

x - f{t; (t),u(x(t)),i£(t)} . (1 1 1 - 5 ) 

and was the equation of a closed-loop system. 

The nonlinearities of the system were due in part to the inertia 
terms of the fundamental dynamical equations. The kinematic variables 
(angular positions and velocities) were nonlinear also. The external for- 
ces, especially the aerodynamic forces contained inherent nonlineari ties . 
Nonlinearities were also introduced into the system by the aerodynamic 
coefficients. However, in this project, the biggest nonlinear contribu- 
tion was made by the control system of the feedback loop, the pilot [6], 

To model the pilot control system is very difficult. A pilot model 
must be able to perceive rates, predict attitudes, and then provide the 
proper lead time in control inputs so that it can maintain some degree of 
precision. The control policy should be one that is capable of driving all 
the state variables to some reference condition, which, for hover, is the 
zero vector. For a nonlinear system, as was the X-14, the derivation of 
such a controller was a difficult task without first producing a linear 
model. For large disturbances as those considered in this study, many 
linear models would be required with each one linearized about a different 
operating condition. The analysis would then switch models as the condi- 
tions warranted. This approach would optimize the use of the control ef- 
fort available, but it would be extremely difficult to construct. There- 
fore, the control policy of [1] was used because of its simplicity. 

Reference [1] calculated the control vector by setting x = 0 and 
solving for u (see Appendix 5 for the equations). This control would not 
return the aircraft to its original position, but would control the air- 
craft to some steady-state reference. For example, if a disturbance caused 
the aircraft to roll, the present control would seek to drive the roll 
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acceleration, p, to zero. This meant that the roll rate would be a con- 
stant and that the roll angle, <f>, would continue to increase at a constant 
rate urvtil the critical roll angle was exceeded. Using this type of con- 
trol would tend to give results on the conservative side for the calcula- 
tion of maximum permissible gusts. 

The pilot introduces a time delay into the system due to his time 
response, which includes his reaction time and the time required for the 
pilot's response to be transmitted to the physical control system through 
servos [6], Reference [7j gave a nominal reaction time delay of 0.15 
seconds. A neuromuscular lag of 0.10 seconds was also given. 

The pilot represented the majority of the time delay. The response 
of the control to a command input was taken to be 0.05 seconds. Equation 
(o) was then written as 

x = f{t;x(t) ,u(x(t-T) t ) } . (1 1 1 -6) 

where t was the time delay. 

The control system just mentioned was affected by the time delay of 
the pilot. The control vector u was calculated every 0.30 seconds. During 
this response time interval, the aircraft was assumed to be in a state of 
flight in free conditions with locked controls. 

A second control was used in which the pilot was assumed to apply 
full control when a disturbance was sensed. Full control in this sense 
meant that full control effort was available to him through the stability 
augmentation system. The throttle for the thrust was assumed to be con- 
trolled just enough to keep the airplane flying. The thrust and the vane 
angle(s) were calculated for x. = (L 

This type of control, u. , . was more representative of a pilot's re- 
, action to disturbances. The design was very basic in nature and could not 
be regarded as sophisticated enough to be useful in other applications. 
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Figure 4 shows a schematic of the control policy using pitch as an example. 
A disturbance (in the form of a step) was introduced to the system and the 
pilot reacted (after the time delay) by applying full control effort to 
counteract the disturbance. This should have resulted in an immediate 
cnange in sign of the pitch acceleration. Full control effort was applied 
until the pi ten rate passed through zero. Then the control was returned 
to the "other" control. The pitch angle was then watched to see if it ex- 
ceeded its critical value. 

Tne preceeding u. control sequence was only applied to step dis- 

“UlclX 

turbances. A simpler sequence of applying u^ ax was used for disturbances 
other than steps. Whenever the aircraft exceeded a certain attitude, usu- 
ally one half the critical value, the control was changed from to 

a . When the attitude dropped below one half its critical value, then 
—max 

“max was returned to 

For the remainder of this study, the two control "schemes" will be 
designated as 


reference 

control 


i = 0. 

LL‘ - 

( 1 1 1-7a ) 

maximum control effort 

a 

-max . 

( 1 1 1— 7b) 


The nine parameters of the control vector were the net engine thrust 

(T ne t)» the two diverter angles (5, A), the three reaction control angles 

(6 ,5 ,5,), and the three aerodynamic control angles (6 ,6 ,6 ). For hover 
x y z a g r 

and low speed flight the aerodynamic controls have negligible effect on 
controlling the aircraft. 

The exhaust side vane angle, A, was set to zero because the present 
aircraft does not use it. Since the primary function of the side vane was 
to produce lateral forces for lateral translation, some simulations were 
conducted with it present. These simulations were restricted to the case 
of lateral disturbances only. Also, only the hover condition was considered. 
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Therefore, for the major part of this study, the nine-dimensional control 
vector was actually a vector of only five dimensions. 

The disturbance inputs represented the three orthogonal components 
of the wind and the three orthogonal components of the gusts. As discussed 
in more detail in the next chapter, the airplane was assumed to be facing 
into a headwind with no crosswinds or vertical winds. Therefore, the dis- 
turbance vector was reduced by two dimensions to a four-dimensional vector. 

The valid range of the model was dependent on four variables [16]. 

The variables with their limits were 

4 

forward velocity (u) - +101 fps (+60 knots), 

-17 fps (-10 knots) , 
side velocity (v) - z51 fps (±30 knots), 

angle of attack (a) - ±20 degrees, and 

sideslip angle (3) - ±90 degrees . 

The conditions for hover were defined as 

forward velocity (u) - +30 fps, 

-10 fps, 

side velocity (v) - ±30 fps, and 

vertical velocity (w) - ±10 fps . 

2.5 Stability Criteria. 

The objective of this project was to find an envelope of wind gusts 
in which the aircraft could operate safely. Therefore, certain criteria 
for instability were established. In hover, close to the ground, the at- 
titude of the aircraft is always very critical. It was assumed that the 
aircraft would be considered unstable, or unable to recover to its refer- 
ence flying condition, if either the roll or pitch angle exceeded 25 de- 


grees. 
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Since the control effort available was limited, the aircraft was also 
considered unstable whenever it exceeded its maximum control effort. The 
maximum available control efforts were 


$ (roll) 

- 

1.361 

_2 

rad-sec 

9 (pitch) 

- 

0.698 

rad-sec 2 

$ (yaw) 

- 

0.332 

rad-sec” 2 


The altitude for the aircraft was initially set at 100 feet. This 
gave the aircraft reasonable room to operate vertically and eliminated the 
altitude as a critical criterion in determining unsafe flying conditions. 

A change in altitude would be an important factor to watch since this would 
be critical at lower altitudes on the order of a few feet. Also, this al- 
titude was assumed to be great enough to eliminate any effect the ground 
might have on a hovering VTOL (e.g., re-ingestion of engine exhausts). 

2.6 Comment on the Aircraft Model. 

As stated in Section 2.2, the aircraft model was based on the work of 
[1] which in turn was based on the work of [16]. One aspect of the latter 
work raised certain questions regarding the aerodynamic forces and moments 
that appeared in some of the equations. In particular", the points in 
question concerned the angle of attack, a. Reference [16] included the 
pitching moment due to a in its model, but it neglected the lift and drag 
due to a. Figure 5 was reproduced from [16] and shows the moment, lift, 
and drag coefficients with respect to a. Three things were immediately 
obvious. First, unlike most aerodynamic procedures, these coefficients 
were given as functions of relative velocities and not as functions of a. 
Secondly, the drag coefficient had a negative value which indicated that 
the drag force acts as a propulsive force (i.e., a contribution to the 
thrust) and not as a resistive force. Thirdly, there existed some ambi- 
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guity in the heading of the graph. The statements "a=0" and "body axes" 
were left to be defined by the user. Digressing slightly in order to inter- 
pret this, [15] showed that the wings were at an incidence of 11 degrees 
at the root and 6 degrees at the tip with respect to the aircraft center 
line (body axes). This then meant the statement "a=0" implied that the 
aircraft was "pitched down" at an angle of 6 to 11 degrees (11 degrees 
will be used for the illustration). To the author, the statement "body 
axes" implied that the forces and moments represented by the coefficients 
were defined to be parallel to the body axes. However, the forces of lift 
and drag, L w and D^, were generated perpendicular and parallel to the rela- 
tive velocity vector (Figure 6). For a=0, this then meant that lift and 
drag were defined with respect to the mean aerodynamic chord. Therefore, 
lift and drag had to be rotated into the body axes system through an angle 
equal to the incidence of the wing with respect to the body axes (i.e., 11 

degrees). Rotation of these forces, and D ro , showed that drag with re- 
spect to the body axes was a positive quantity (Figure 6). This contra- 

dicted the information presented in Figure 5. At the present time, only 
a partial copy of [16] is available. The previous discrepancies may be 
resolved when a complete copy is referred to. 

For this study, the original equations as defined in [1] were used 

in which was included and and were neglected. However, a set 

of simulations, with C, and C,. included in the aircraft model, was con- 

La Da 

ducted. Also, since the atmosphere and aircraft reference frames were 
parallel (see Chapter 4), the angle of attack had an initial value of 11 
degrees. Therefore, another set of simulations was conducted where 11 de- 
grees was added to the angle of attack. The results of both sets of simu- 
lations are presented in Section 5.3.c. 
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Figure 6. Interpretation of "a=0" and "body axes". 
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CHAPTER 3 
WIND MODEL 


3.1 Background. 

A hovering aircraft, either a VTOL or a helicopter, is limited to 
hovering within the so-called Earth boundary layer. This boundary layer 
extends from the surface of the Earth through the lowest few hundred feet 
of the atmosphere. Due to the nature of many factors present in the 
planetary boundary layer, turbulence is nearly always present [9]. 

Turbulence is defined by [8] to be rotational, dissipative, three- 
dimensional, nonlinear, stochastic, diffusive, and a continuum phenomenon. 
Turbulence can then be defined as a random process that cannot be described 
by explicit functions of time; only a statistical, probabilistic approach 
can be taken [5]. Therein lies the problem of modeling the wind. 

Reference [5] stated that much of the extensive information avail- 
able on wind-induced turbulence near the ground is inconclusive and even 
contradictory. This made a low altitude turbulence model seemingly im- 
possible to derive. However, with certain assumptions, a reasonable re- 
presentation of presently available information was constructed. 

The first assumption was that turbulence can be broken down into a 
slowly-time-varying component and a rapidly- time-varying component. This 
was seen by observing a typical spectrum of wind speed near the ground 
{Figure 7). The spectrum showed that measured winds contain high-frequency 
and low-frequency modes with a wide gap of frequency where the wind con- 
tains little energy. Because of the long periods associated with the low- 
frequency mode (on the order of hours), this low-frequency mode was assumed 
to be constant over the relative short period of time (10 seconds) of the 



Cycles/hour 0.01 0.10 1.0 10 100 

Hours 100 10 1.0 0.10 0.01 


* tv=frequency 

S u (n) -power spectral density 


100 

0.00 


Figure 7. Schematic Spectrum of Horizontal Wind Speed Taken as a Function of Time at a Fixed Point. 
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computer simulation. This then defines a constant velocity mean wind. 

The high-frequency mode was assumed to represent the irregularly fluctua- 
ting wind gusts. 

Another assumption was that there was no dependence of the statis- 
tical properties of turbulence on time. Turbulence was therefore assumed 
to be a stationary process. Turbulence was also assumed to be homogeneous 
(i.e., the statistical properties were the same at each point in the gust 
field). The last assumption was that the gust velocity components were 
Gaussian. This was necessary in the development of the gust part of the 
model in Section 3.3. 

3.2 Wind. 

It was stated in Section 3.1 that the low-frequency wind was assumed 
to be a constant because of the short time period of the simulation. This 
was also carried over and used to assume that the direction of the wind 
would be constant over the time period of the simulation. 

The reference frame of the wind was an atmosphere- fixed reference 
frame (see Chapter 4), the axes of which point north, east, and vertically 
down. Hie wind model assumed the mean wind to be blowing from north to 
south (a negative direction) parallel to the north-south axis of the re- 
ference frame. Therefore, the crosswind and vertical wind were zero. 

The mean wind was assumed to be constant over the entire width and 
height of operation of the aircraft (see Figure 15). 

3.3 Gusts. 

A discrete model of gust velocity is described in [13], The discrete 
model has the "1 - cosine" shape and is shown in Figure 8. The equation 
of the discrete gust model is 
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v(x) = 


X 

{1 - cos{tf • — )} . 


( 8 ) 


This model is an arbitrary model that relates the gust magnitude in a 
rational manner to the expected intensity of continuous random turbulence. 
The model has a gust velocity v(x) defined spatially in terms of a magni- 
tude 'v , which occurs at a distance x=d , where d is a physical dimension 
m m m 

of the gust velocity field. This model can be applied to any of the three 
gust components [13]. This model describes an average of all conditions 
for clear air turbulence. The model neglects the effects on turbulence 
of terrain roughness, atmospheric stability (lapse rate), mean wind mag- 
nitude, and all other meteorlogical factors, except altitude. Too few 
data are available to incorporate these factors into the model. 

The validity of a discrete model for gust magnitudes was discussed 
in [2] and [13]. Reference [13] stated that a discrete gust provides 
spike- type inputs that may not be apparent in the simulated Gaussian ran- 
dom turbulence. These gusts affected a vehicle in a specific way, but 
one which was likely to be encountered. 

The discrete gust model jointly considers gust magnitude and gust 
gradient, both of which are important parameters. A study of these para- 
meters was the objective of this research. 

The discrete model was a function of three parameters, d^, already 
mentioned, L, the scale length of turbulence, and'a, the root-mean-square 
intensity of continuous random turbulence. The last two parameters were 
functions of altitude. A conditional probability density function of the 
random turbulence was the key to the discrete model. The derivation of the 
model is presented in [13]. However, a basic description and some of the 
assumptions used will be repeated here. 

The discrete model was developed by arbitrarily choosing reasonable 
values of the scales and then determining values for the intensities so 
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that the mathematical spectral form matched the measured spectra data. 

The root-mean-square intensity of vertical velocity, a s was the basis 

w 

for the model. The statistical properties of a were developed in [13]. 

w 

First, the discrete probability (P^ of encountering turbulence was 

obtained from Figure 9 (Figure 3, page 443 of [13]). As can be seen, it 

is a function of altitude. The value of a at a particular altitude was 

w 

taken to be that rms turbulence which was exceeded with a probability of 
exactly 0.01 (i.e., 99 percent of all time spent in flight at a given al- 
titude will be spent in either turbulence with less than the specified o w 
or in turbulent free air). If P(a v/ ) is the probability that a equals or 

exceeds a given value, then by definition of a , 

w 

P(cr w ) = P] • P(ff w ) = 0.01 (1 1 1 — 9 ) 

where P( a ) and P. are defined as above and P(a ) is the conditional prob- 
w 1 w 

ability of equalling or exceeding a given a once turbulence has been en- 

w 

countered. The function P(a ) was obtained from Figure 10 (Figure 2, page 
442 of [13]), which was based on fitting a Rayleigh distribution to known 
gust data. From the graph, the value of P(c ) necessary to make P(a ) = 

0.01 was found and the corresponding a , as a function of altitude, was 

W 

found . 

The vertical scale of turbulence, L , had been established to a 

w 

reasonably good approximation as a function of altitude. Reference [13] 
cites two different scales for clear air turbulence, the preferred von 
Karman scales and the Dryden scales. The vertical scale of turbulence was 
found from one of the following equations, 

L = 2500 feet altitude (h) > 2500 feet , (III-10a) 

or L - h feet altitude (h) < 2500 feet . (III-10b) 

The scales of turbulence for longitudinal, L u » and lateral, !_ v> turbulence 
were calculated from one of the following equations. 
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L = L = 2500 feet altitude (h ) > 2500 feet , (III-10c) 

u v 

or L = L * 184 . f ee t altitude (h) < 2500 feet . (lll-lOd) 

U v 

The root-mean-square intensities for the longitudinal and lateral clear 
air turbulence were derived from the von Karman relationship 


2 2 2 

!u _ v _ W 

, 2/3 " 75/3 75/3 

L u V L w 


4n-n) 


For any component of clear air turbulence, the scale of turbulence and the 

rms intensities were known. By arbitrarily selecting several values of 

d , corresponding values of v were obtained from Figure 11 (Figure 7, page 
m m 

429 of [13]). Figure 11 was derived by cross-plotting the values of v m /0 

and d /L from Figure 12 that corresponded to a probability of occurrence 
m 

of 0.01. Figure 12 is a Gaussian cumulative distribution, P{v /a), in 

a normalized form. The variable d /L was a parameter used in calculating 

m 

P(v /a). Therefore, for a given altitude and a given probability of being 
m 

in a specific turbulent field, a discrete model was found. 

This was a partial solution to the problem of determining the maxi- 
mum gust fields in which the aircraft can safely fly. The approach taken 

in this research was to find the envelope of maximum v 's and, then, 

m 

knowing the lengths, d^, from basic relationships regarding each v^, 

working back to find the probabilities of encountering such turbulence. 

A change was made in the discrete equation (8) in order to apply 
it to the problem. The procedure used in [13] was for an aircraft flying 
through the velocity field. The total time the aircraft was subjected 
to the turbulence was the length of the gust field, 2d m , divided by the 
airspeed of the aircraft. For this project, the length of time the X-14 
was subjected to the turbulence was the distance traveled by the gust 

divided by the speed at which the gust field was traveling. The length of 

time was specified based on previous knowledge about gusts [10]. The pro- 



Normalized discrete gust magnitude - 



Normalized discrete gust length - d x /L u » “V^v* ^z^w' 

[u,v,w] and [x,y,z] refer to longitudinal, lateral, and 
vertical gusts respectively. 

Figure 11. Obtaining v m Knowing L and cr and Selecting d^. 
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Figure 12. Probability of Equalling or Exceeding a Given Gust Magnitude v m for Various Values of d^/L. 
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duct of the time and the speed of the gust field then gave an appoximation 
to the length 2d n (see figure 30). 

Examination of Figure 7 shows that the period associated with the 
peak of the high-frequency component occurs around 0.05 hours (3 minutes). 
However, [10] stated that analysis suggests at least two- thirds of tur- 
bulent energy is associated with fluctuations lasting less than five 
seconds. This latter information was used as a basis for establishing 
testing frequencies of the gusts. Another conclusion was that, at mod- 
erate heights, the eddying energy was equally divided along all three 
axes. This then was the basis for modifying the discrete equation to be 
a function of time rather than length, and for it to be equally repre- 
sentative of any one of the three gust components. Therefore, equation 
(8) became 

v t 

v(t) - — - {1 - cos(ir * — )} , (til-12) 

2 t 

m 

where the distances were replaced by times. The time t was any instant 
of time during the simulation between t=0 and t=10 seconds. The time t 
was the "half-life" of the gust frequency and was when occurred. 

Three “base" times, t , were selected for the model. Times were 

m 

limited to the length of the simulation or less. The "base" times were 
ten, five, and three seconds. These times corresponded to frequencies of 
0.10, 0.20. and 0.33 cycles-sec - ^ (cps) respectively. The fourth repre- 
sentation of a gust, the step, was also used since a step represents the 
most severe type of turbulence an aircraft may encounter [2], [14]. 

3.4 Equations of the Wind Model. 

The final model of the wind consisted of a mean wind with gusts in 
the three component directions. The gust opposite in direction to the 
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mean wind was omitted for the reason that this condition would violate the 
validity of the model in that it would force 0 to be greater than 90 de- 
grees. The mean wind was considered to blow from the north at a constant 
velocity of UWI fps. The gusts were considered positive if they blew in 
the direction of the positive body axes. The mean wind was transformed 
from the atmosphere- fixed reference frame into the body-fixed reference frame 
by the Euler angles. The equations were 

U* -UWI .cos(*).cosfe> , ( III-13a) 

VW bf = UWI f .{cos(^).sin(e).sin(^)-cos(*)*sin(^)} » (lll-13b) 

WW bf * UWI af *{cos((j)).cos(ij))‘Sin(0)+sin(^)*sin((j))} , (lll-13c) 

where the subscripts bf and af referred to body reference and atmospheric 
reference frames respectively. 

The gust components were also transformed into the body reference 
frame of the aircraft by the Euler angles. The equations were 
UG^ = UG a ^»cos(^) •cos(e) + VG^*cos(0) *sin(ijj} - 
WG a ^.*sin(e) , 

= VG a ^.*{cos(^)*cos(ij))+sin(iJ))’Sin(tp)-sin(e)> - 
UG af *{cos((j))»sin(ijj)-cos{i|;)*sin((}>)*sin(0)} + 

WG a ^*cos(0)«sin(<}>) , 

WG, - = UG • {cos (ip) *cos($) *si n(0 )+si n( ip) *sin((J>) } - 
DT dt 

VG a ^.{cos(i|i)*sin(<f>)-cos(<j>)«sin(i{>) •sin(0) } + 

WG »cos((j)) »cos(0) , 

at 

where UG af * l/2*v mu *{l - cosU-t/t^)} » 

VG af = • 

UG af = • 

3.5 Input to the Aircraft Model. 


([ I I~T 4a ) 

(. III-T 4b) 

( III-l 4c) 

(III-l 5a), 
( III-l 5b) 

( III-l 5c) 


The approach to this study is illustrated in Figure 13. The wind 




Figure 13. Problem Approach. 
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model described the atmospheric turbulence, the output of which was a velo- 
city field. This velocity field acted as a disturbance on the aircraft in 
the form of changes in aerodynamic forces and moments. These forces and mo- 
ments were then fed into the aircraft model. The stability of the output, 
or the vehicle motion, of the aircraft model was used in determining the 
maximum gust field the aircraft could withstand. 

The velocity field of the atmosphere may be regarded, over the time 
and space intervals of interest, as composed of a steady mean value with 
turbulent fluctuations superposed. This leads to the assumption that the 
structure of turbulence takes the form of individual patches, in each of 
which the turbulence is approximately random, homogeneous, and isotropic. 
This is based on the fact that the statistical properties of the distur- 
bance input to an airplane flying through a turbulent field are not appre- 
ciably affected by the variation of that field with time. Essentially, 
turbulence may be treated as a frozen pattern in space. 

Another assumption must be made in order to deal with a hovering 
VTOL aircraft like the X-14. Reference [5] stated that an assumption of 
a frozen turbulence model was invalid for a stationary point. However, 

[8] and [12] stated that the hypothesis of a gust field frozen in time 
moving downwind with the mean wind speed (known as Taylor's hypothesis) 
was acceptable at low airspeeds. The lower limit of airspeed at which 
this assumption was valid was put at one- third the mean wind speed. 
Therefore, the model of the aircraft was assumed to have a forward air- 
speed equal to one-third the mean wind speed. This does not violate the 
definition of the hover conditions as a 30 fps forward velocity was de- 
fined as the limiting forward velocity for the hover mode. 

It was assumed that there was no variation of the gusts over the 
physical dimensions of the aircraft. The longitudinal gust was assumed 
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constant along the wing span, the lateral gust was assumed constant along 
the fuselage, and the vertical gust was assumed constant along the wing span 
and fuselage. The aircraft was in effect treated as a point [5]. 

Treating the airplane as a point simplified the model as far as defin- 
ing the input to the aircraft. However, it restricted the usefulness in 
calculating certain responses. The point approximation was valid provided 
the wavelengths of the gusts were much greater than the physical dimensions 
of the aircraft. This was valid for lower frequencies only. Reference [5] 
limited the upper frequency for validity to 2/(wing span) which, for the 
X-14, was 0.059 cps. This was the upper limit for considering the plane 
vanishingly small compared to the gust wavelength. The smallest frequency 
used in the model was 0.100 cps, or almost twice the limiting value. This 
would indicate that the point approximation assumed in this study would 
yield somewhat limited results. 

To consider the higher frequencies would then require treating the 
airplane as a finite plane in space where the gusts vary along the dimen- 
sions of the aircraft. Variations in longitudinal gusts along the wing span 
would result in induced rolling moments and pitching moments. Variations 
in vertical gusts along the wing span and fuselage would result in pitching 
and rolling moments also. Variations in lateral gusts along the fuselage 
would result in yawing moments. 

The gust model would need to be modified in order to produce these 
variations. The modifications would entail splitting the gust velocity 
field into smaller parts, each of which possessed different magnitudes and, 
each of which acted on a separate section of the aircraft. The aircraft 
model would then also need modification in order to calculate the response 
of the aircraft to each gust input. These modifications would increase 
the complexity of each model, which, at this time, is beyond the scope of 
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this project. 

Treating the aircraft as a point also neglects the effect of gust 
penetration (i.e., wing-to-tail delays). Gust penetration is important 
for high gust magnitude to airspeed ratios, which was the case for this 
project. However, this would require that the response of the tail to 
gusts be known. At the present time these responses are not available. 

For the altitude considered in this project, 100 feet, this model 
derived in the foregoing was assumed to be a reasonable representation 
of clear air random turbulence. 
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CHAPTER 4 

METHOD OF ANALYSIS 

4,1 Reference Frames. 

Four frames of reference were needed for the definition of this prob- 
lem. An inertial reference frame was needed in which Newton's Second Law 
is valid for motion of a particle or a rigid body. The inertial frame 

chosen was as Earth-fixed reference, F , with Earth surface axes 0 x y z 

e e e^e e 

(Figure 14). The origin of this frame was placed near the vehicle with 
the 0 g z e axis pointing vertically down, the 0 g x e axis pointing north, and 
the 0 y axis pointing east. The rotation of the Earth and its curvature 
were neglected, thus being called a Flat-Earth approximation (i.e., treat- 
ing the Earth as a stationary plane in inertial space). 

A reference frame defining the aircraft was also needed. In this 
case, two were defined. First, a vehicle- carried vertical frame, ‘F , with 

axes 0 xv 2 was attached to the aircraft at the mass center. The 0 z 
v rv v v v 

axis was directed vertically down and the remaining axes, and 0^, 
were chosen to point north and east respectively. For hover, the movement 
of the aircraft from its original position was small enough to consider 
the origin of F^ near enough to the origin of F g so that the axes of each 
could be considered parallel. 

The second reference frame associated with the airplane was a body- 
fixed reference frame, F^, with body axes Oxyz. The origin of this frame 
was located at the aircraft mass center (see Figure 1). The axes were the 
same as those defined in Chapter 2. As mentioned previously, the Euler 
angles gave the orientation of the body axes relative to F . However, 
by considering the origins of F & and F y to be coincidental, or nearly so, 
because of hover, the Euler angles also gave the orientation of the aircraft 
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relative to the inertial reference frame, F L relative to F . 

b e 

The fourth reference needed to define the problem was one to be used 

with the motion of the atmosphere. The atmosphere- fixed reference frame, 

F , with axes 0 x y z , was necessary to define the relative velocity be- 
a a a a a 

tween the aircraft and the atmosphere for the calculation of aerodynamic 
forces and moments. When the atmosphere was at rest relative to the Earth, 

F and F were the same. When the atmosphere was in non-uniform motion 

d 6 

relative to the Earth, F was chosen so that the average motion of the 

d 

atmosphere relative to F was zero [5]. Therefore, non-uniform and uniform 

d 

motion of the atmosphere relative to the Earth resulted in motion of F q 

relative to F . By selecting F parallel to F initially, F then moved 
g a e a 

parallel to F at constant velocity, 
e 

All reference frames were assumed to have the same origin at the start 

of the simulation, or at time zero, as shown in Figure 15. The changes in 

the references for an incremental time At are shown in Figure 16. Reference 

F has moved a distance of UWI»At feet in the negative x-direction. The 
d 

aircraft has flown along a trajectory path P and has undergone angular dis- 
placements in roll, pitch, and yaw. For purposes of derivations to follow, 
these displacements were assumed to be positive. The aircraft was no longer 
orthogonal to the other three reference frames. Disturbances parallel to 
F were transformed into F h by the Euler angles. The aircraft now "saw" 
orthogonal components of the wind and the gusts. 

The transformations were carried out in a particular sequence as shown 
in Figure 17 [5]. The sequence was 1) a rotation ^ about 0 y z v carrying the 
axes to a temporary position 0 x^ 2 z 2* 2 ) a Nation ^ about 0^ carrying 
the axes to another temporary position 0 x^y^Zg, and 3) a rotation 6 about 

0 x carrying the axes to their final position Oxyz. 
v 3 


4.2 Method of Solution. 






> 

Wind conditions are 
constant over the 
height and width of 
operation of the aircraft. 



Figure 15. Relative Positions of the Reference Frames for Initial Conditions. 
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Figure 17. F v ancJ F a Transformed into F^ Through 
the Euler angles 8, <f>, and ip. 
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The objective of this research was to find the largest gust field the 
aircraft could encounter and still maintain a stable flying condition. The 
stability of the aircraft was defined by certain constraints or limits on 
some of the state variables and control parameters. 

The procedure was to solve equation (6) 

x = f{tjx(t),u(x(t-x)),itf(t)} , 

for maximum w with constraints and limits imposed on x and u. Due to the 
nonlinearity and time variance of the equations, tv could not be solved for 
directly. 

A simple iteration technique to solve for w indirectly was used. An 
arbitrary wind model was selected and introduced into the aircraft model. 

If the airplane flew in a stable manner for ten seconds, the wind model was 
increased in magnitude. If the airplane could not recover, the magnitude 
was decreased. The iteration was carried out until the largest values of 
wind and gusts the aircraft could withstand were found. 

The problem was then an initial value problem with a known set of 
initial conditions and a set of 12 first-order differential equations that 
could be integrated by numerical methods. 

The initial conditions of the aircraft represented the hover mode. 

The forward velocity, u, was set equal to one-third the mean wind speed, a 
necessary condition for Taylor's hypothesis. The remaining state variables, 
except for the altitude, were set to zero. 

The length of time for the observation of the aircraft was arbitrarily 
set at ten seconds. The solution of the differential equations was obtained 
by integrating over the time interval in steps of 0,10 seconds. The inte- 
gration routine used a fourth-order Runga-Kutta iteration to calculate the 
initial values of the integration and a four-point Adams -Bashforth -Moulton 
predictor-corrector method to continue the integration. 
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The first part of the solution to the problem was to determine the 
maximum steady-state headwind the airplane could withstand. No other 
components of the wind were calculated because of the assumption that the 
aircraft was always facing initially into the steady-state headwind. The 
second part of the solution consisted of finding operational envelopes of 
maximum gust magnitude for various gust frequencies, headwinds, and con- 
trol policies. 

4.3 Computer Program. 

The flow chart for the program WINDY, used to calculate the opera- 
tional envelopes, is shown in Figure 18. The program was given the initial 
conditions for hover and the velocity of the headwind. For these condi- 
tions, the aerodynamic coefficients used in calculating the aerodynamic 
forces were calculated. These aerodynamic coefficients changed instantan- 
eously for each change in relative velocity of the atmosphere with respect 
to the aircraft. 

fit some time t the gust disturbances were introduced. These gusts 
were parallel to the atmosphere-fixed reference frame. The program trans- 
formed the gusts and the wind into the aircraft reference frame. 

The initial values of the control vector were then calculated by 
WINDY. These values were stored by the program. They represented the 
state of the aircraft before any disturbance was encountered. 

The solutions of the system of differential equations describing the 
aircraft motion were obtained and stored for a particular instant of time. 
Next, the angular accelerations were calculated and stored. These values 
were obtained by treating the equations involving angular motion (i.e., 
the p, q and r equations) as algebraic equations. For the particular 
time (of the simulation) in question, these equations take the form 
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Figure 18. Flow Chart for Windy. 
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Figure 18. (continued) 
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Pj = * (i I 'Ii- 16a) 

qj = f(Xj.«£j»«ij) , and (I II -16b) 

r. = . (I II -16c) 

The 'current values of the state variables, x., the control parameters, u., 

J J 

and the disturbances, w. , were substituted directly into the above equa- 

J 

tions. These equations also gave the control effort requirements since 
these requirements and the angular accelerations were equivalent. 

The stability criteria were checked as the program proceeded to the 
next part. If the aircraft was found to be unstable, the program printed 
the time history of the motion up to that point and then stopped. The time, 
the state variables, the control parameters, and the angular accelerations 
were printed for every 0.10 seconds. If the aircraft was stable, the time 
was checked to see if it was time to calculate new values for the control 
parameters. The control parameters were calculated every third step after 
the introduction of the disturbances. This simulated the time lag due to 
the pilot. 

Upon every calculation of a new control vector, the program checked 
the saturation limits of the various parameters. Any control parameter 
exceeding its saturation point was set equal to that saturation point. 

The time was then incremented by 0.10 seconds. If the time was less 
than ten seconds, the program introduced the gusts with their new velo- 
cities and then went through the same procedure as just outlined. If the 
simulation was complete (ten seconds) the program printed the complete 
time history. 

A copy of WINDY is presented in Appendix 6. 
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CHAPTER 5 

RESULTS, DISCUSSIONS, AND CONCLUSIONS 

5.1 Steady-state Wind (Headwind). 

The first part of the project was to find the maximum headwind (with 
no gusts) in which the aircraft could maintain a hovering flight condition. 
The computer simulations found the maximum headwind to be 61 fps (36 knots). 
The aircraft was found to remain virtually motionless with only the thrust 
and the pitch control nozzle angle, differing from hover conditions. 

The net thrust was lower by 449 pounds (obviously due to the lift created 
by the wind) and <5 was nearly saturated at -19.2 degrees. 

V 

The next higher value of headwind, 62 fps, resulted in the aircraft 
being unable to counteract the pitching moment created by the lift due to 
the wind. The pitch control angle was saturated (-20 degrees) for the en- 
tire simulation. The aircraft exceeded +25 degrees in pitch attitude after 
5.6 seconds. 

5.2 Maximum Gust Envelope. 

Table 1 presents the maximum gust envelope for the X-14. The values 
in this table represent the variable v m of equation (12) in Chapter 3. 

5. 2. a Analysis of Longitudinal Gusts. 

Table 1 shows the longitudinal gust magnitude (UG) to be independ- 
ent of frequency. All three frequencies have nearly the same magnitude, 

52 fps. The computer simulations showed that <3^ became saturated whenever 
the gust magnitudes approached their maximum values (52 fps). 

An analysis of the next higher magnitude for each frequency showed 



HEADWIND LONGITUDINAL GUST LATERAL GUST VERTICAL GUST FREQUENCY 


UWI 

-UG 

+VG 

-10 

-52 

+65 

-10 

-52 

+65 

-10 

■*53 

+72 

-10 

— 

+27 

-61 

maximum headwind 



-VG 

+WG 

-WG 


-67 

+23 

-19 

0.10 cps 

-66 

+23 

-17 

0.20 cps 

-74 

+23 

-21 

0.33 cps 

-27 

+16 

-15 

step 


All magnitudes are fps. 


Table 1. Maximum Permissible Gust Envelope. 
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the aircraft to exceed +25 degrees in pitch for the two lower frequencies 
(0.10 and 0.20 cps). The pitch control effort, equivalent to q, was ex- 
ceeded in the third case. 

5.2.b Analysis of Lateral Gusts. 

Examination of Table 1 shows that the permissible magnitudes for the 
lateral gusts to be the largest of the three component magnitudes. The 
limiting factors were the aircraft side velocity and the pitch angle, de- 
pending on the frequency content of the gust. At 0.10 and 0.20 cps, the 
side velocity exceeded the definition of hover for all gust magnitudes in 
excess of 65 fps. At 0.33 cps the pitch angle exceeded +25 degrees for 

j 

gusts in excess of 72 fps. It should be noted that if the side velocity 
was ignored, the maximum permissible gusts were 88 and 73 fps for 0.10 
and 0.20 cps respectively. 

Detailed analysis showed that a large yaw angle, \p, was indirectly 
responsible for the large pitch angle (0.33 cps case and for the two cases 
where side velocity was ignored). For a large i p t typically 50 to 60 de- 
grees, the lateral gust provides a component which resembles longitudinal 
velocity, U Q . Figure 19 shows the relationship between U Q and 9 for a 
longitudinal gust of 53 fps and 0.10 cps (minimum case of instability). 

Also shown are U Q and e for a lateral gust of 89 fps and 0.10 cps (minimum 
case of instability also--one fps above the limiting case). The figure 
shows a close resemblance between the two sets of curves as ip increases. 
This substantiates as being indirectly responsible for 9 exceeding its 
critical value. 

The lateral gusts in the form of a step function resulted in more 
limited responses because of excessive yaw control effort (ij> or r) demanded 
by the aircraft. The permissible magnitude was 27 fps. Ignoring aircraft 
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(fps) 



Time (seconds) 


Figure 19. U Q and e for Minimum Unstable 
Longitudinal and Lateral Gusts. 
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side velocity, the permissible magnitude was 53 fps. 

5.2.c Analysis of Vertical Gusts. 

Table 1 shows the aircraft to be the most sensitive to the vertical 
gusts. In all cases the critical parameter in determining instability was 
9 . For both positive (down) and negative (up) vertical gusts, the pitch 
angle exceeded +25 degrees before any of the other hover criteria were vio- 
lated. The maximum gust magnitudes were +23 and -19 fps for 0.10 cps, +23 
and -17 fps for 0.20 cps, +23 and -21 fps for 0.33 cps, and +16 and -15 fps 
for the steps. 

Figure 20 shows the time history for 0 for +23 and -19 fps (0.10 cps), 
the limiting hover case. The responses were normal for the first half of 
the simulation. Normal responses meant positive pitch from a down gust 
(gust striking a larger surface area behind the center of gravity of the 
aircraft and thus forcing the nose up) and negative pitch from an up gust 
(same reasoning, nose down). The initial -0 created by the up gust resulted 
in an increase in u, forward velocity. This increase results in a larger 
positi ve pitching moment which should tend to return the aircraft to an 
"even keel". Figure 20 shows this to be the case over the first six seconds 
of the simulation. However, at that time, the aircraft acquired a large 
positive pitch rate which rapidly drove the aircraft unstable. The pitch 
control effort did not respond as would be expected for this situation. 
Figure 21 shows 5 for both cases. As can be seen, 6 never approached 

V V 

saturation at any time for the up gust case. The reason for this can be 
found by examining the equation for <5^ (equation (A5-5), Appendix 5). In 
detailed analysis it was found that the angle of attack exceeded its valid 
range for the model, ±20 degrees. Typically a approached 50 to 60 de- 
grees. This was the cause for the unusual results. 
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(rad) 



Figure 20. 0 for Positive and Negative Vertical Gusts. 
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Figure 21. 6 y for Positive and Negative Vertical Gusts. 
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PRECEDING PAGE BLANK NOT FILMED 


(rad/sec 2 ) 



Time (seconds) 

term (i) (l^J^l/I^J.p-U^S-c^C^+AC^'a 

term (iii) T net (j)*(z 2 *sin(S{j))+x 2 *cos(o(j)))/I y 
term (11) T^U^AC^-c/I and term (iv) M . *6 /1 

y y y y 

Figure 22. Analysis of Selected Terms of q (up gust case). 
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q increased through term (iii) which resulted in an increasingly larger 
0 . As e increased (positively), d increased again and the cycle started 
again. This phenomenon can be seen by examining the equation for 5, 

a = sin"^(U re i/ v jet * A/T net (j) " M-sin(a)) . <in-17a; 

The second term, A/T nfit (j), was the dominating term of the equation. 

Examining it more closely revealed that e and 5 were closely related. 

The term, with the proper substitution for A, was 

•‘•-{m-(vr-w*q-g-sin{0))+l/2*p-Up‘S*(C Lo -sin(a)-C D *cos(a))}/T net (j)-“ 

The dominant part of this expression was m*g*sin(0 )/T ne ^.(j) . This term 
reduced essentially to sin(e) since m *g/T net (j) was nearly unity because 
m*g, the weight, was canceled by the thrust, which, for hover, was slightly 
less than the weight. Therefore, the expression for a was approximately 

a = sin _1 (.**sin(e)--0 (lll-17b 

which reduced to 

a*e (III-17c 

for large positive 9. 

This showed the model to be inherently unstable for any condition 
that resulted in a positive pitch angle. It must be remembered that 
this occurred because of the way the control vector, of which 5 was a part, 
was defined. The next two subsections attempted to alleviate this problem 
with a by restricting a through gust combinations and by utilizing more of 
the control effort available. 

5.2.C.1 Analysis of Longitudinal and Vertical Gust Combinations. 

The large angles of attack were responsible for the initial positive 
pitching motion. These large angles resulted when the vertical relative 
velocity component, ^ re i » was large compared to the relative longitudinal 
velocity, U Q . The equation for a. 



198 . 


a = sin" 1 (W rel /U 0 ) , (III-18) 

then forced the ratio w re -j / u 0 to be restricted to a value less than ±0.342 
in order for a to be within ±20 degrees. Placing a restriction on ^ re 7 /^ 0 
complicated the iterative procedure used to find the vertical gusts of the 
gust envelope. This was abandoned in favor of observing the response of 
the aircraft to "typical" gust fields. In these gust fields the longitudi- 
nal gust was also considered since it was an integral part of U Q (i.e., 

u o vu rel rel ' 

Reference [8] stated that the only successful non-zero cross-corre- 
lation between different velocity components recorded at the same points 
involved the longitudinal and vertical gust components. The procedure 
of Section 3.3 was used to find "typical" gusts for both components. These 
combinations were used to give an estimate of the maximum magnitude of a. 

Only the combinations that kept ct close to ±20 degrees or less were used. 

Table 2 gives the 10 combinations that were used. The gust combinations 
were calculated on probabilities of occurrence of 0.10, 0.01. and 0.001. 

Two different headwinds, 10 and 25 fps, were used in determining the 
gust combinations. The total effects of the two headwinds will be dis- 
cussed in more detail in a later subsection. 

The aircraft was unable to withstand the gust field for the last 
three cases of Table 2. In all three situations, the pitch angle exceeded 
+25 degrees. Table 3 gives the maximum values of the pitch angle for the 
10 cases. The data showed the aircraft to be more stable for the combina- 
tions with positive vertical gusts than those with negative vertical gusts. 
This was true only for the lower value of headwind (cases a to d. Table 2). 

The behavior of the aircraft "switches" for the higher headwind as the 
aircraft was more stable for those combinations with negative vertical gusts 
than for those with the positive vertical gusts. The degrees of stability 
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HEADWIND 

LONGITUDINAL 

VERTICAL 

PROBABILITY OE 

APPROXIMATION 

ACTUAL 



GUST 

GUST 

OCCURRENCE 

OF ALPHA 

ALPHA 

a. 

-10 

-8.4 

-7.2 

0.10 

21° 

16° 

b. 

-10 

-8.4 

+7.2 

0.10 

-21° 

-16.5° 

c. 

-10 

-19.1 

-16.3 

0.01 

28° 

21° 

d. 

-10 

-19.1 

+16.3 

0.01 

-28° 

-22.5° 

e. 

-10 

-29.8 

-25.4 

0.001 

30° 

* 

f. 

-10 

-29.8 

+25.4 

0.001 

-30° 

* 

g. 

-25 

-11.5 

-7.2 

0.10 

11° 

7.5° 

h. 

-25 

-11.5 

+7.2 

0.10 

-11° 

-T 

i. 

-25 

-25.7 

-16.3 

0.01 

18° 

16.5° 

J * 

-25 

-25.7 

+16.3 

0.01 

-18° 

-11.5° 

k. 

-25 

-40.0 

-25.4 

0.001 

21° 

37.5° 

1 . 

-25 

-40.0 

+25.4 

0.001 

-21° 

-13° 


* these two were not used. 

all gust frequencies were 0.10 cps; all magnitudes were fps. 


Table 2. "Typical" gusts, longitudinal and vertical combinations. 


CASE MAXIMUM PITCH ANGLE CHANGE IN ALTITUDE MAXIMUM 6 


a. 

-.041 to .117 rad 

+0.51 ft 

+3° to -4° 

b. 

.083 rad 

+1.26 ft 

-6.7° 

c. 

-.084 to .379 rad 

+4.22 ft 

+13° to -12° 

d. 

.202 rad 

+3.10 ft 

-15.2° 

g- 

-.040 rad 

-0.64 ft 

+2° to -1° 

h. 

.098 rad 

+4.81 ft 

-9.1° 

i . 

-.016 to .060 rad 

+4.53 ft 

+1.5° to -9° 

j* 

>25° at 6.4 secs 

+5.45 ft 

-20° 3. 9-6. 4* 

k. 

>25° at 6.2 secs 

+6.20 ft 

-20° 3. 9-6. 2* 

1 . 

>25° at 4.6 secs 

+3.41 ft 

-20° 3. 0-4. 6* 


*time during which 


S was saturated. 

*7 


Table 3. 0, Az, and 6 for “Typical" Gust Combinations. 
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here were based upon the relative absolute magnitudes of 8, which was the 
critical parameter. The three conditions at which the aircraft became un- 
stable were characterized by a rapid increase in a over the Tast few sec- 
onds of the simulations. Table 4. Again the aircraft model was shown to 
be inherently unstable due to the interdependence between a and 8. For 
these cases, the larger magnitude for the headwind created the positive 
pitch which led to the instability of the aircraft. 

5.2.C.2 Use of All Available Control Effort, u . 

Since the unstable cases for vertical gusts (limiting cases) involved 
0 , iimax was applied only to Whenever 8 exceeded 12.5 degrees (1/2 of 
its critical value), the control switched from to (i.e., 6 =±20 
degrees depending on the sign of 8). Table 5b compares the results of 
using u^ ax against not using it for the gusts of Table 5a. The effect of 
u^ ax was to reduce the maximum pitch rate and thus, allow the operational 
envelope to be expanded. The control was more effective for up gusts than 
down gusts. Figure 23 shows that a. did a better job in controlling 0 
by using more of the control effort available. The other control, 
never reached a saturated state and therefore, did a poorer job in control- 
ling 0. At this point it must be stated again that these limiting condi- 
tions of vertical gusts were beyond the valid range of the model in that a 
exceeded ±20 degrees. 

The results of using u in determining maximum vertical gust mag- 

“HliaX 

nitudes are presented in Table 6, which shows relative values between the 
vertical gusts of Table 1 and those using a . The results were affected 
by the fact a exceeded 20 degrees for a good portion of the time. However, 
Table 6 does demonstrate the effectiveness of using all the control effort 


available. 
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CASE 

TIME 

a 

j. 

3.3 secs 

.111 rad 


6.4 secs 

.474 rad 

k. 

4.8 secs 

.121 rad 


6.2 secs 

.411 rad 

1 . 

2.7 secs 

.117 rad 


4.6 secs 

.501 rad 


Table 4. Rapid Increase of o for Unstable Cases j, kj and 1. 


CASE 

HEADWIND 

VERTICAL GUST 

FREQUENCY 

a' . 

-10 fps 

+23 fps 

0.10 cps 

b‘ . 

-10 fps 

-19 fps 

0.10 cps 

c' . 

-10 fps 

+23 fps 

0.20 cps 

d\ 

-10 fps 

-17 fps 

0.20 cps 

e' . 

-10 fps 

+23 fps 

0.33 cps 

f'. 

-10 fps 

-21 fps 

0.33 cps 

g'- 

-10 fps 

+16 fps 

step 

h'. 

-10 fps 

-15 fps 

step 


Table 5a. Vertical Gust Magnitudes (repeated from Table 1). 


CASE 

MAXIMUM 

0 AND TIME, (r sQ 

MAXIMUM 

9 AND TIME : 

’ ^ax 

a'. 

.425 

rad 

at 

10 secs 

,335 

rad 

at 

7.3 

secs 

b'. 

.413 

rad 

at 

10 secs 

.266 

rad 

at 

9.1 

secs 

c' . 

.365 

rad 

at 

10 secs 

.328 

rad 

at 

8.9 

secs 

d'. 

.401 

rad 

at 

10 secs 

.223 

rad 

at 

7.3 

secs 

e'. 

.324 

rad 

at 

10 secs 

.290 

rad 

at 

8.5 

secs 

f. 

.411 

rad 

at 

10 secs 

.236 

rad 

at 

7.3 

secs* 

g*. 

.432 

rad 

at 

10 secs 

.240 

rad 

at 

4.9 

secs 

h\ 

.402 

rad 

at 

10 secs 

.238 

rad 

at 

7.5 

secs 


f o 

*q exceeded -0.698 rad/sec at 7.3 secs. 

Table 5b. Comparison of 0 max between and u^ ax - 
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The analysis of Table 6 showed that, with the use of u. . the air- 

nnaX 

craft became less sensitive to positive gusts and more sensitive to nega- 
tive gusts (step functions excluded) as the gust frequency increased. The 
results for the negative gusts showed the aircraft to call for more thrust 
than it was capable of producing (thrust saturation), to acquire a high for- 
ward velocity, and to call for radical changes in u^ ax » figure 24. As the 
figure shows, the forward velocity and the pitch angle diverge while the 
net thrust and the reaction control angle change radically. 

Step functions for vertical gusts were also investigated. The con- 
trol policy for u^ ax described in Section 2.4 was applied to this type of 
disturbance. Table 7 compares the maximum vertical gusts for (dif- 
ferent than the u max used with the step functions of Table 6) and u^_g. 
Tables 6 and 7 show that some form of control by a type of policy was 
definitely beneficial. 

5. 3. a Model Variation, the Effect of the Side Vane Angle, X. 

The present configuration for the X-14 does not include the side vane 
angle. The side vane angle was introduced in the study to see its effect 
in combatting the high side velocities induced by lateral gusts. The la- 
teral gusts were limited to positive-valued step functions and the control 

was u.« _n . 

“2C-U 

The resulting analysis is presented in Table 8. The maximum gust mag- 
nitude dropped slightly from 53 to 45 fps. The side vane angle reduced the 
side velocity of the aircraft considerably at a consequence of increased 
rolling and yawing moments. A slight loss of lift was noted since the 
aircraft required higher thrusts with X. 


5.3.b Model Variation, the Effect of Increasing the Magnitude of the Headwind. 
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FREQUENCY 

MAXIMUM VERTICAL GUST 

U‘ „ 

-x=0 

MAXIMUM VERTICAL GUST 
%ax 

0.10 cps 

+23 fps 

+23 fps 

0.10 cps 

-19 fps 

-35 fps 

0.20 cps 

+23 fps 

+28 fps 

0.20 cps 

-17 fps 

-26 fps 

0.33 cps 

+23 fps 

+30 fps 

0.33 cps 

-21 fps 

-26 fps 

step 

+16 fps 

+16 fps 

step 

-15 fps 

-15 fps 

Table 6. 

Maximum Vertical Gusts 

for- a^o and w 


MAXIMUM VERTICAL GUSTS (STEP FUNCTION) 

-i=0 +16 fps -15 fps 

^max +22 fps " 19 fps 

Table 7. Maximum Vertical Gusts (Step Function) for and . 


X INCLUDED WITHOUT X 


Maximum lateral gust velocity (step) +45 fps 
Maximum aircraft side velocity, v +3.23 fps * 

Maximum roll angle, <j> .402 rad 

Maximum yaw angle, \p , -.453 rad 

Change in lateral displacement. Ay 12.5 ft 

Change in altitude, Az „ 3.2 ft ? 

Maximum yaw control effort, ^ -.195 rad/setr 

and time of occurrence 0.2 secs 

Minimum net thrust 4194 lbs 


+53 fps 
+58.34 fps 
.201 rad 
-.250 rad 
265 ft 
24.2 ft 
-.322 rad/sec 
0.2 secs 
4122 lbs 


Table 8. Comparison of Lateral Mode Parameters for Control 
With and Without the Side Vane Angle. 
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Table 9 compares the maximum allowable lateral gusts (ignoring the 
limitations on aircraft side velocity) and vertical gusts for two headwinds, 
10 and 25 fps. Increasing the headwind had little effect on the verti- 
cal gusts, except for the step functions which showed slight decreases in 
magnitude. The pitch angle was the limiting factor in determining sta- 
bility. Again, as with the other simulations involving vertical gusts, 
it must be remembered that a plays an important role when interpreting 
results. 

The larger-valued headwind increased slightly the operational enve- 
lope of the aircraft when it was subjected to lateral gusts. The excep- 
tion was the step functions where an increase in headwind resulted in a 
big decrease in gust magnitude. For the step functions, the limiting fac- 
tor for stability was the yaw control effort. Large requirements for it 
occurred during the time delay before the control was initiated. 

The 25 fps headwind was felt to be a practical maximum which would 
be encountered in actual conditions. Hence, it was felt that there was 
no need for analyzing aircraft responses at still higher magnitudes. 

5.3.c Considerations of C La- W and (a+11 0 ). 

Section 2.6 pointed out the ambiguities and uncertainties of 
C^, Sla* anc * 01 as usec ^ ky D6]. The major concern was centered on the 
fact that [16] used in its model, but neglected and C^. The fol- 

lowing simulations were run to test the effects of adding these two coeffi- 
cients to the model. These coefficients appeared in the equations for 
forward and vertical velocities and were present in calculating the thrust 
requirements (see Appendix 7). 

The first simulation was concerned with the effect of and Cp a 
on the response of the aircraft to longitudinal gusts. The effect was 
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GUST HEADWIND 

(-lOfps) 

+VG (step) 

+53 fps 

+WG (step) 

+16 fps 

-WG (step) 

-15 fps 

+VG (0.10 cps)* 

+88 fps 

+VG (0.20 cps)* 

+73 fps 

+VG (0.33 cps) 

+72 fps 

+WG (0.10 cps) 

+23 fps 

+WG (0.20 cps) 

+23 fps 

-WG (0.10 cps) 

-19 fps 


(Table 1) HEADWIND (-25 fps) 

+10 fps 
+11 fps 
-13 fps 
+95 fps 
+78 fps 
+76 fps 
+22 fps 
+23 fps 
-19 fps 


♦ignores aircraft side velocity. 


Table 9. Comparison of the Limiting Lateral and Vertical 
Gusts for Two Different Headwinds. 
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negligible. The same result occurred when C, and C n were added to the 

La Da 

simulations involving gust combinations. Results of these two series of 
simulations are presented in Table 10. 

Section 2.6 mentioned that the aircraft possessed an initial angle of 
attack of 11 degrees when the aircraft and atmosphere reference frames were 
parallel. The last series of simulations were conducted with and C Da 
included in the model and with the angle of attack increased by a constant 
of 11 degrees [5], The simulation with the longitudinal gusts showed that 
the addition of 11 degrees to a resulted in slightly higher forward and 
vertical velocities and significantly lower thrust requirements (Table 11). 
The lower thrusts were a result of higher lift generated by the wing due 
to the higher angle of attack. As a consequence, the simulation also re- 
sulted in a slightly larger pitch angle which in turn required more pitch 
control effort (see Figure 25). 

The gust combinations that were considered were combinations (i) 
and (j) of Table 2, longitudinal gust of 25.7 fps and 0.10 cps and vertical 
gusts of ±16.3 fps and 0.10 cps The effect of (a+ll°) was to reverse the 
responses from the original responses (without a+ll°). The aircraft remain- 
ed stable when subjected to the down gust and it became unstable when sub- 
jected to the up gust. Table 12 compares the simulations without the addi- 
tions of Cj^, C^, and (a+11 0 ), the simulations with and only, and 
the simulations with C^ a , Cp^, and (a+ll°). Figures 26a and 26b show 0 
for the positive and negative vertical gusts. The figures show the re- 
versed responses due to (a+ll°). Figures 27a and 27b show the effects of 
(a+ll°) on the net thrusts for each case. 

5.4 Interpretation of Wind Gusts. 


The data presented in the first three subsections of this chapter 
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LONGITUDINAL GUST ONLY GUST COMBINATIONS 

LONGITUDINAL MAGNITUDE, 25.7 fps 
VERTICAL MAGNITUDE, ±16.3 fps 
FREQUENCY, 0.10 cps 
HEADWIND, 25 fps 

The most significant differences 
occurred in T t . 


25.7, +16.3 25.7, -16.3 


time 

without with 

without 

with 

(sec) 

(lb f ) Ob f ) 

(lb f ) 

(lb f ) 

0.0 

4130 4130 

4130 

4130 

4.2 

3942+ +4402 



5.1 


3720+ 


5.4 



+3022 

6.1 

+3861 



6.4 

4095+ 



10 


411 9+ 

+3985 


Table 10. Effects of C La and C Qa on Longitudinal 
Gusts and Combination Gusts. 


LONGITUDINAL GUST ONLY: MAGNITUDE, 19 fps; FREQUENCY, 0.20 cps 
HEADWIND, 20 fps 

without C La , C Da , (a+1 1 0 ) with C La> C Da , (a+ll°) 

6.67 fps 7.38 fps 

.285 fps at 2.8 secs -.496 fps at 2.7 secs (minimums) 

.000 fps at 5.2 secs .008 fps at 5.2 secs (maximums) 

.286 fps at 7.6 secs -.577 fps at 7.6 secs (minimums) 

.002 fps at 10 secs -.041 fps at 10 secs (maximums) 

4182 Ibf at 0.0 secs 4046 lbf at 0.0 secs (initially) 

3994 lbf at 2.7 secs 3716 lbf at 2.7 secs (minimums) 

4182 lbf at 5.1 secs 4056 lbf at 5.1 secs (maximumx) 

3993 lbf at 7.5 secs 3664 lbf at 7.5 secs (minimums) 

4180 lbf at 10 secs 4032 lbf at 10 secs (maximums) 

Table 11. Comparison of Relative Maximums and Minimums 
of u, w, and T ^ With Regard to (a+ll°). 



MAGNITUDE, 19 fps 
FREQUENCY, 0.20 cps 
HEADWIND, 20 fps 

No significant differences in 
any parameters between using 
or not using and C^ . 



2 


6 


8 


Time (seconds) 


note. Longitudinal gusts only. 



2 4 6 8 

Time (seconds) 


Figure 25. Effect of (a+ll°) on 0 and 6 y . 
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A: LONGITUDINAL GUST: 25.7 fps and 0.10 cps 

VERTICAL GUST: -16.3 fps and 0.10 cps 

B: LONGITUDINAL GUST: 25.7 fps and 0.10 cps 
VERTICAL GUST: +16.3 fps and 0.10 cps 


A. 

B. 


A. 


B. 


A. 


B. 

A. 

B. 


A. 

B 


without 

C La* C Da* 

wUh C La* C Da 01 ' 1y; 

w1th C La* W (o+in 


ok 

ok 

e>25° at 8.3 secs 

8>25° 

at 6.4 secs 

0>25° at 6.1 secs 

ok 

(“max 5 

8.33 fps 

10.2 fps 

17,8 fps 

( Vx> 

8.33 fps 

8.33 fps 

8.33 fps 

fyjiax 

-.62 fps 

-1.67 fps 

-3.00 fps 

( w max^ 

-.74 fps 

-1.07 fps 

-.595 fps 

< T net> 

4130 at 0.0 secs 
3720 at 5.1 secs 
4119 at 10 secs 

4130 at 0.0 secs 
3022 at 5.4 secs 
3985 at 10 secs 

3927 at 0.0 secs 
2246 at 6.0 secs 
3190 at 8.3 secs 

< T net> 

4130 at 0.0 secs 
3942 at 4.2 secs 
4095 at 6.4 secs 

4130 at 0.0 secs 
4402 at 4.2 secs 
3861 at 6.0 secs 

3927 at 0.0 secs 
3650 at 6.6 secs 
3788 at 10 secs 

% 

■'max 

) -8.9° at 5.7 secs 

9.2° at 5.1 secs 

+20% 3.6 to 6.5 

<«v 

) -20°, 3.9 to 6.4 

-20% 3.9 to 6.1 

-20% 4.8 to 5.6 


max 


Table 12. Effects of C, , C n , and (ct+ll°) on Combination Gusts. 

La Da 



[I] 


0 . . , 
original 


(a+ll°) 


Time (seconds) 


(a) Longitudinal gust. -25.7 fps, 0.10 cps 
Vertical gust: +16.3 fps, 0.10 cps 


original 


Figure 26. 


Time (seconds) 

(b) Longitudinal gust: same as (a) 

vertical gust: -16.3 fps, 0.10 cps 

Effect of (a+ll°) on 6 for Combination Gusts (+vertical gusts) 
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represented the maximum magnitudes of gusts for the discrete model. Real 
world wind gusts do not appear in such a well-behaved manner. Section 3.3 
defined a method for determining root-mean-square intensities from the 
discrete model. The only detail lacking was some measure of the variable - 
2d m . This variable was interpreted as the distance between the points 
where the gust velocities were zero (see Figure 8). The measure of this 
length was affected to a great extent by the speed at which the frozen gust 
field was moving. By considering the airplane as a point, only the move- 
ment of the gust field past that point was needed. A first approximation 
to 2d would.be the velocity of the gust field times the time duration of 
that field (Figure 28a). However, the point moves also and its displace- 
ment must be taken into account (Figure 28b). The last figure shows that 
if the displacement was in the same direction as the gust field was moving, 
the displacement was subtracted from the first approximation. For the oppo- 
site case, the displacement (of the point) was added. 

Now, at this point the speed at which the gust field is moving takes 
importance. For longitudinal gusts, this speed was assumed to be the mean 
wind speed. The longitudinal displacement of the aircraft was then needed 
to completely determine 2d . Once 2d was obtained, the rms intensities 
for longitudinal gusts could be determined. Table 13 presents the rms in- 
tensities for the longitudinal gusts of Table 1. The rms intensities were 
calculated for probabilities of occurrence (for 100 feet in altitude) of 
0.10, 0.01, and 0.001. This meant, for example, that for a rms intensity 
of 41,6 fps (24.7 knots), which was the lowest rms intensity for longitud- 
inal gusts, a gust with a magnitude of 52 fps (v from Table 1) would occur 
only 0.10 percent of the time. The conclusion drawn was that the aircraft 
could withstand a severe longitudinal gust environment. A sample calculation 
for rms intensity is presented in Appendix 8. 
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t=0 seconds 


t=lO seconds 


2d m = 10 f P s * 10 seconds 
- 100 feet 


(a) P stationary. 



If 26^-100 feet, point P will not be completely 
passed in 10 seconds. 

Therefore, 2d ffl must be altered in order to com- 
pletely pass P in 10 seconds. 


The displacement of P must be subtracted from 
2d m ^ ‘ e *’ 2d m * 2d m' " x ^* 


(b) P mobile. 

Figure 28. Interpretation of 2d . 

m 
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v ffi (of Table 1) rms Intensity 



.10 (10%) 

.01 (1%) 

.001 (0.1%) 

52 fps 0.1 cps 

84.0 

54.7 

41.6 

52 fps 0.2 cps 

115.5 

74.3 

57.7 

53 fps 0.33 cps 

151.5 

88.4 

75.7 


Table 13. rms Intensities for Longitudinal Gusts. 


PROBABILITY OF OCCURRENCE 

0.10 ( 10 %) 0.01 ( 1 %) 0.001 ( 0 . 1 %) 


0 u 

a v 

a w 

a u 

84.0 

60.0 

39.0 

54.7 

115.5 

82.5 

53.6 

74.3 

151.5 

108.1 

70.4 

88.4 


39.2 25.4 41.6 29.8 19.3 

53.0 34.5 57.7 41.2 26.8 

63.1 41.0 75.7 54.1 35.2 


Table 14. Lateral and Vertical rms Intensities, Based 

on Longitudinal Intensities, Using o u /a v /a w =2.8/2.0/l .3. 
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The length 2d did not exist for the lateral and vertical qusts be- 
m 

cause of the assumption that there did not exist any mean winds in these 
directions to transport the gust fields. Therefore, the rms intensities 
cannot be calculated for these gusts. Reference [12] stated that nominal 
values for component gust intensity ratios differ among the various liter- 
ature. However, [12] picked nominal values for the longitudinal /lateral/ 
vertical (o^/a^/a^) gust ratio of 2. 8/2. 0/1 .3 to agree with available data. 
Table 14 presents the lateral and vertical gust intensities based on this 
ratio and on the a^'s of Table 13. The smallest intensity for a vertical 

gust was 19.3 fps rms. From Figure 10 it was determined that the proba- 

-5 

bility of a equalling 19.3 fps rms was much smaller than 10 . There- 

fore it was concluded that the likelihood of the aircraft encountering a 
velocity field possessing enough energy to upset it was very remote. 

5.5 Closing. 

The results presented above must be compared against some standard 
in order to assess them. The pilots of the X-14 have arbitrarily set 
limits on wind conditions. They will not fly in winds over 12 knots (20.3 
fps— longitudinal direction) and they would prefer not to fly in winds over 
eight knots (13.5 fps). On research flights, data measurements will not be 
taken if the winds exceed five to eight knots [17]. 

Based on the assumptions made for this project (listed in Table 15), 
the following conclusion was made. The arbitrary limits set on the wind 
conditions by the pilots were within the limits determined by this study. 
The limits were 

52 fps, 0.10 cps 

52 fps, 0.20 cps 

53 fps, 0.33 cps for longitudinal gusts; 
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+65 and -67 fps, 0.10 cps 
+65 and -66 fps, 0.20 cps 
+72 and -74 fps, 0.33 cps 

±27 fps, step for lateral gusts; and 
+23 and -19 fps, 0.10 cps 
+23 and -17 fps, 0.20 cps 
+23 and -21 fps, 0.33 cps 
+16 and -15 fps, step for vertical gusts. 
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Assumptions, Aircraft Model 

1. Rigid body. 

2. Forward velocity equals one-third the mean wind velocity (initially). 

3. No variations in the gusts over the physical dimensions of the 
aircraft. 

4. No gust penetration effects. 

5. constant over the entire range of operating speeds of the en- 
gi nes . 

6. Considering the angular accelerations equivalent to the reaction 
control efforts required. 

7. Excluding C^ a and in the model (as does [1] and [16]) has 
negligible effect. 

8. A 0.30 second delay exists between the introduction of a distur- 
bance and the application of control. 

Assumptions, Wind Model 

1. Turbulence is stationary, homogeneous, and isotropic. 

2. Turbulence is a frozen pattern in space. 

3. Gust velocity components are Gaussian. 

4. The direction and velocity of the headwind are constant. 

5. Effects of terrain roughness, lapse rate, mean wind magnitude, 
etc. on turbulence are neglected. 

6. Periods of less than ten seconds represent most of the turbulent 
energy. 

7. Gust fields move at a speed equal to the mean wind speed . 

8. No crosswinds or vertical winds (steady-state winds) exist. 

9. Estimation of 2d is reasonable. 

m 

Table 15. Assumptions Concerning Aircraft and Wind Models. 



APPENDICES 
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APPENDIX 1 

Modification of the Rolling Moment Equation. 


The rolling moment, t, due to side velocity, v, was omitted from the 
set of original equations for the X-14. This characteristic of the airplane, 
known -as the dihedral effect, is represented by the stability derivative 
C„ in the roll equation (p equation). 


The derivation of C followed a procedure used in [2]. The stability 

3 

derivative is a function of five parameters: 1) influence of sWeepback, 

2) influence of fuselage, 3) influence of the tail fin, 4) contribution of 
the wing planform, and 5) influence of dihedral. 

1. Influence of sweepback of the wings. 

The X-14's wings are swept forward 4.6° at the quarter-chord line [15]. 
The interpretation of the effect of sweepforward (or sweepback, depending 
on the situation) is that a change in lift between the two wings occurs 
when the aircraft is yawed with respect to the wind. The derivative C 

*6 


is then a function of the lift coefficient C^. 

The difference in lift of two panels between the two wings is 
&L = C L ■ 1/2 * P • 1/2 * S • V 2 • {cos 2 ($ - A) - cos 2 (p + A)} (lII-Al-1) 

where AL is the change in lift between the two panels, 
is the local lift coefficient for the panels, 


p is the density of the atmosphere at sea level, 

S is the wing area, 

V is the local relative velocity between the panels and the atmosphere, 
3 is the angle of sideslip, and 
A is the sweepforward at the quarter-chord. 

Assuming p a small angle, equation (Al-1) reduces to 



223 . 


AL = C L • 1/2 • p • S • V 2 • b . sin(2A) . (III-A1-2) 

The rolling moment produced by the panels will be 1/2 the difference 
AL multiplied by the distance between the centers of pressure of the two 
panels 

L = 1/2 • d • C L • 1/? • p * S • V 2 • 3 • sin(2A) . (III-A1-3) 

The corresponding rolling moment coefficient is 

1/2 • d • C. • 1/2 • p • S • V 2 • B * sin(2A) 

C * , (III-A1-4) 

1/2 • p • • S • b 

which reduces to 

1/2 . d * C, • B • sin(2A) 

C. (lII-Al-5) 

L b 


where d is the distance between the centers of pressure of the panels and 

b is the wing span. 

The stability derivative C 0 is 

*8 

dC. 1/2 * d • C, • sin(2A) 

— =C • (lII-Al-6) 

88 8 b 


The total effect would be the integration of the changes ip lift over 
the entire wing span. For purposes of this derivation, it will be assumed 
that C„ is a representative average of the C. 's of all the panels. The 

A e e 


d for this average will be assumed to be the distance between the mean 
aerodynamic chords (approximately 2/5 of $he wing span). Therefore 

1/2 . 2/5 * b , C. • sin{2 • 4.6°) 

r* *— 


«* 0.031976 * C L . (lII-AI-7) 

2. Influence of the fuselage. 

The wing-body interference effect on C 9 can be calculated from 
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MO Z w h + w 

AC 0 = 1.2 • A ' ^ • (III-A1-8) 

b b 

where A is the aspect ratio of the wings, 

z w is the vertical distance the wing-root quarter-chord point is below 
the fuselage center line, 
b is the wing span, 

h is the average fuselage height at the wing root, and 

w is the average fuselage width at the wing root. 

The parameters 2 ^, h, and. w were estimated from drawings in [15] to have 

values of 1 foot, 3.36 feet, and 4.38 feet respectively. Therefore 

-.in 1 4.38 + 3.36 

AC 0 = 1.2 • (6.2 y ft 

*8 33.83 33.83 

= 0.0202 (III-A1-9) 

3. Influence of the tail fin. 

When the aerodynamic center of the vertical tail fin is appreciably 
offset from the roll axis (x-axis of the aircraft), the side force on the 
vertical surface (side force due to sideslip) may produce a significant 
contribution to the roll moment. For hover, this contribution was assumed 
to be negligible. 

4. Contribution of wing planform. 

The contribution of wing planform on C is found from the lower 
figure of Figure B.11.2, page 489 of [2]. The lower figure is for taper 
ratios (A) of 0.50. The taper ratio of the X-14 is 0.48. For A equal to 
6.2 and Ay^ (sweepforward of the wings at the quarter-chord) of -4.6° 

(negative for sweepforward), the value of -(C ) /C. is approximately 

W L 


= -0.0275 • C L . 


0.0275. Therefore 


(lII-Al-10) 
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5. Influence of dihedral. 

The influence of dihedral on C. is found from Figure 11.3, page 490 

3 * 

of [2]. The top curve gives C 0 for zero sweepback. This can be convert- 

3 

ed to C for the dihedral angle involved by using the formulas 

p 

( r 1 

v 8,^ ‘ partial span r 

C £ = • {r r T > (lII-Al-11 ) 

3 3 r (r ) 

1 'full span r 


where C is given by the equation 

r A + 4 • cos(A) 

C = (C ) A=Q , (lII-Al-12) 

B r (A + 4) • cos ( a) h 

A is the aspect ratio, and 

Tj is the dihedral of the partial span. 

For the X-14 r T is zero. Therefore, for A equal to 6.2 and a taper 

ratio of 0.50, the top figure gives a value of (C^ ) a= q of -0.0162. 

3 r 

Then 1 


6.2 + 4 * cos(-4.6°) 

= — * (-0.0162) 

e r (6.2 + 4) • cos(-4.6°) 

“ -0.0162 (III-A1-13) 

and 

C. = -0.0162 (2 - 0) 

6 

= -0.0324 {l 1 1— AT -14) 



The total of all contributions on C 0 is 

3 

C. = 0.032 • C. + 0.0202 - 0.0275 • C, - 0.0324 

C = 0.0045 • C, - 0.0122 . 

3 


*** 


(lII-Al-15) 
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APPENDIX 2 


Calculation of the Diverter Vane Efficiency Factor. 


The comparison of results for the hover mode (no disturbances) between 
the simulated and the actual results showed a discrepancy ih the engine 
speeds. Simulation results gave an engine speed at 90.02 percent of maxi- 
mum speed whereas actual results gave the engine speed at 98 to 99 percent 
of the maximum speed. It was then concluded that the efficiency factor for 
the diverter vane had been omitted from the model. 

Reference [1] calculated this factor by simply dividing the theoreti- 
cal speed by the actual speed. Hence 


90.02 

r = = 0.9139 . 

98.50 


(III-A2-1) 


However, the equations of the model calculate the engine' speed from 
the engine thrusts since the thrust was solved for explicitly as one of the 
control parameters. Therefore, if the efficiency factor was based on the 
thrusts, the value was 

4194 

% 0.80654 . (III-A2-2) 

a T 5200 


The thrusts 4194 and 5200 were obtained from Figure 29, which shows 
the total thrust as a function of engine speed. Figure 29 was derived from 
Figure 30, which shows the thrust curves of each engine. 

The model uses this efficiency factor in that it takes the calculated 
net thrust (obtained from the control vector) and multiplies it by the in- 
verse of K~ (1.23986). This thrust is then sent to a subroutine which 
a T 

calculates the engine speed from the total gross thrust. The equations are 

Voss = ’• Z3986 • Vt (III-A2-3) 


for the gross thrust. 
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fl = 4.52 x]0" 9 -(T~4370) 3 -9.31x10“ 6 -(T-4370) 2 +].25x10' 2 *(T-4370)+92 (III-A2-4) 

for the engine speed of gross thrusts greater than 4370 pounds, and 

0 - (T - 3305)/88.75 + 80 (III-.A2-5) 

for the engine speed of gross thrusts greater than 3305 pounds, but 
less than 4370 pounds. 

The exit velocity, V., of the engine exhausts was 2000 fps. The 

J 

diverter vane efficiency also reduces this value to E,~ *V., or 1613 fps. 

Oy J 

This is based on 


T * m • V. 

gross j 

at the engine exhausts. The net thrust is lg ross *?g so that 

Cl I I-A2-6 ) 

T gross * ~ ^net ^ ^ % T ’ 

The mass flow rate, m, is constant. Therefore 

(III-A2-7) 

t, = m * V. 
net jo 

where Vj 0 is 1613 fps. 

(EII-A2-8) 
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APPENDIX 3 

Modification of AC^ and AC Mo of the Aircraft Model. 

The aerodynamic coefficients AC^ andAC^ o were defined initially by 

[16] to be functions of U » relative longitudinal velocity, and T , thrust 

0 c 


coefficient, which 

, in turn, was 

a function of T . , the net thrust, 
net 

The 

equations for aC^ 

and aC m were 
Mo 





AC l 

= T . f(U ) 
c v o' 

and 

(III-A3-1) 


AC Mo 

= T • f^U ) 
c o' 


(III-A3-2) 

where 

T = T 
c 

+ /(l/2.p.U 2 . 
net ' H o 

s) . 

(III -A3- 3} 


Therefore, T ^ must be known in order to calculate AC. and AC M . 

net L Mo 

However, the net thrust contained, among other things, the term AC L . This 
posed a problem in that T^ et could not be calculated without first knowing 
aC^, and also, aC^ and aC Mq were to be calculated before T net was to be 
calculated. The procedure through the program was: 

i) Define initial conditions; x(0) = 

ii) Calculate the aerodynamic coefficients for x.(0). 

iii) Calculate the control vector for the initial conditions; 
a(0) = tt(x( 0)). 

iv) Introduce the disturbances w(l). 

v) Calculate the new values for the state variable vector; 

x(l) = f(x(0),a(0),w(l)). 

vi) And so on. 

T . was not calculated until (iii) whereas AC, and AC W/ , were calculated 
net L MO 

at (ii). One way to solve this dilemma was to estimate T ngt at (1). This 
would then require that T net be estimated for each different set of wind 
conditions. A second solution was found through algebraic manipulation of 
the equations involving aC^ and AC^. These equations were, for AC^, 
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forward velocity, u = vr-wq-g«sin(e)+(T / m )-{sin(a)- CO s(X)- 

U rel /V jet^ + ( 1/ 2 ) * ( 1/ m )-p - U p * 5* { { C^+AC^) * (III-A3-4) 

sin(a)-C D * cos (a)} , and 

vertical velocity, A = u*q-vp+g*cos(8)«cos(<j))-(T net /m)*{cos(a)* 

c °s(X)-W rel /V jet }.(l/2).(Vm)-p-u2-S-((C Lpo+ fin-A 3 - 5 ) 

AC^)‘COs(a)+C D *sin(a)} . 

F ° r AC Mo* 

pitch rate, q = { ( I z -I x )/I y } • p- r+( 1^/ I y ) • (r z -p 2 }-2- { I e / 1 ) -Q-r+ 

(T net /I y ) ' ((W rer x r U rer z l )/V jet +cos(>) ’ (z 2 ‘ sin(S)+ 
x 2 .cos(3)} + (M q /I y ).q + (H { /y-y(l/2).(l/I y ).p.u2. {m. A 3. 6) 

S ‘ C ’ (C Hopo +AC Mo t( C «,pA> •“-( 1 /2) • ( 1 /U o ) • C C Hq ' 

6 

By substituting equation (A3-3) into equations (A3-2) and (A3-1), and then 
substituting (A3-1) into equations (A3-4) and (A3-5) and then (A3-2) into 
equation (A3-6), the equations for u, w, and q were 

u = — +(1/2) •(1/m) - P -u|-S-{{C Lpo +T net /{l/2.p.U^*S)*f(U o )>*sin(a)‘([II-A3-7) 
w = ...-(l/2).(l/m). p .U^.S.{{C Lpo +T net /(V2.p*U^.S) *f(U Q ) }‘Cos(a) • ([II-A3-8} 
q = ...+(l/2).(l/y • P *U^.S-c.{C Mo p o +{T net /(l/2*p-U^-S) -f 1 {U Q } }+ — CTII-A3-9) 
The quantity (l/2)»p»U^»S cancelled and the resulting expressions, 

( T ne t* AC L* sin ( a ^/ m and ^net^Slo’^^y* were incor P° rated into the net 
thrust terms of the equations. The resulting equations were 

n “ *{*’ ,+ ^^*sin(a) }+(l/2) • ( 1 /m } •p*U^*S»{C^ po *sin(a)-* • » 


w = 


&U-A3-10) 

• •*(^ ne -{./ ni ) * { • * * + AC^*cos(a) }-(l/2) «(l/m) •p*U^ , S* {C^ po *cos(a) +• • • 

&II-A3-11) 


q “ - (T net /I y ) -(- +AC M o- c } t - +( 1 / 2 )-( 1 / I v)^' U o- S - c -( C Mopo + 


4 II-A3- 12) 


where AC l and AC^ were no longer functions of T c and U Q , but U Q only, 
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that is ac l - f(U 0 ) and AC = f ! (u >. Therefore, the estimation of T , 

u net 

was not needed now. 
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APPENDIX 4 


Equations of the System and Terms of the Equations. 

The equations of the twelfth-order system are presented in this ap- 
pendix. 



U = vr-wq-g*sin{0)+(T net /m)*{s1n(a)-cos(X)-U re1 /V jet +AC L *s1n(a)}+(l/2)*O/m)*p-U^S'{C L •s1n(a)-C D *cos(o)} 


= w*p-u*r+g-sin{<|))*cos(0)+(T net /m)*(sin(X)-V rel /V j . et )+(l/2)*O/m)*p-D^*S-C^ 


(III-A4-1 ) 
(III-A4-2) 


w « u-q-vp+g*cos(<J>)‘cos(e)-(T net /tn)-{cos(cf)*cos(X)-W rel /V jet +AC L ‘COs(a)}-(l/2)*(l/m)-p»U^S*{C Lpo *cos(a)+C D -sin(a)} 

(III-A4-3) 

P = < I x -I 2 /(I x ,I z- I x 2 >}-{( { V I z } / I x ),c l -r+ < I xz /I x )- P' t l + < T net /I x) - {t z l ,v rel )/v jet- sin ( x )- {z 2 +1 r cos ( 9 »l + (V I x ) ' p+ 

(L « x /I x )- V (1/2),<1/I x>’' :> '^ -S - b -(V C M a - V C 2 6 ‘ e W I x2 /I )( }-{(a x -y / I z ) - p -q-(I X z/I z )-q-r«-(I e /I 2 )-!!-q+ 
(T net /I z ) 'f s,n(x) ' {x 2- l r s<n(5)} -( ;< r V re1 )/V Jet J+(N r /I z ) ' r+(N « /I 2 )-S z +(1/2,-(1/I z ) ' 0-i ^' S ' b ‘ (C no +C n« " S r + - 


V‘ 5 a>» 

3 


(III-A4-4) 


i “ ({I z‘ I y >/I y ) ' p ' r+(I xz /I y ) ’ (r2 ' p2, ' 2 ' (I e /! y ) ‘ !! ' r+<T net /I y ) ' <(W rer )t r U r e r z l )/V jet +c< ’ stx),{z 2 ‘ s1,,(5)+x 2 -cos(5)}+ 


c - A C Mo }+(M q /I y )-q + (M (S /I y )-V {1/2) -( 1 / I y ) -P ,U o- S - c -< C Hopo t( W AC fb ) ‘ a+(1/2,-(1/U o ),C ‘ C Mq ,, ' +C m ’ 5 e) 


( III-A4-5) 


r * ^ I x‘ I z /(I x ,I z' I xz^*^ I xz /I z^^ {I y" I z }/I x ) * q * r+(I xz /I x ),p ’ q+(T net /I x ) *^ z l’ V rel )/V jet“ s ' 1n(x) ‘ {z 2 +£ l‘ cos(5)} ^ 
(L„/I x )-P+(L 5 /I x )-V (1/2 >-< 1/I x>- |5, ' J 0 ,S - b -( C 20 +C 2S -V C 11 -S)Wn x -I >/I z )-p-q-(I xz /I z )-q-r. + 2-(I e /I z ).n-q + 


ro 

to 

4*> 



'wy 


'{sin(\)'{x 2 -H 1 -s1n(S)}-(x 1 'V rel )/Vj et }+(N ( ./I z )T+(N 




5 z +(l/2) 




S-b 


t C no +( W 


V 


C nS *V» 

a 

0 = q * cos (cf> ) - r*sin($) 

<i> = P + q*sin(<t>)*tan(9) + r* cos (<}>)• tan (0 ) 

= (r*cos(<f>) + q*sin(<f>))/cos(0) 
z = v*sin(<j))*cos(0) + w*cos(<i>)*cos(9 ) - u*sin(0) 

x = u*cos(9)*cos(^) +■ v*{sin(4»)*sin(0) 4 cos(ip) - cos + w*{cos(<j>)*sin{0)*cos(\j;) + sin(<}>)*sin(iiO} 
y = u«cos(e)»sin{ifi) + v*{sin(4>)'sin(0)«sin(^) + cos(4>)»cos(^)} + w*{cos(4>)*sin(0)*sin(i|;) - s1n(<j>)»cos0fi)} 
where a = sin” 1 (W rel /U Q ) 

e = s1n * 1 < ¥ re1 / V 

U = (U 2 , + W 2 ,) 1/2 
o ' rel rel' 

U = (U 2 , + V 2 , + W 2 ,) 1/2 
o v rel rel rel ; 

U rel = u - UW - UG 

V re1 = v - VV1 - VG 

W rel = w - WW - WG 

The equations for T net > 5> X, <$ x> <5^, and & z are presented in Appendix 5. 


(III-A4-6) 

(III-A4-7) 

(III-M-8) 

(III-A4-9) 

(III-A4-10) 

(I1I-A4-11) 

(III-A4-12) 

(II I-A4- 13} 

( III-A4-14) 

(III-A4-15) 

(III-A4-16) 

( III-A4-17) 

( III-A4-18) 

( III-A4-19) 
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APPENDIX 5 


Equations for the Control Vector. 


The equations for the six control parameters were derived in [1] and 
will be restated here. The equations were derived by setting the x vector 
equal to the 0 vector and solving for the control vector, u. 

The het thrust equation was 

T net = {-E + (E 2 + D f ) 1/2 } /d (EII-A5-1) 

where E = (A-U re1 + &'V re1 + c * w re i )/ V j et " AC L *(A-sin(a) - Ocos(a}} , 

D “ 1 - ^ V !et - AC L + 2-AC L -(U rel -sin(a) - , 

F = A 2 + B 2 + C 2 , 

A = m*{vr-wq-g*sin(8)}+(l/2)»p*U 2 *S»{C L p 0 *sin(a)-C D *cos(a)} , 

B = m*{w*p-u*r+g*sin4)*cos(e)}+(l/2)*p*U^*S*Cy , and 

C = m*{yq-v*p+g*cos(<j))*cos(0)}-(l/2)*p«U 2 *S*{C^p O ‘COs{a)+Cp*sin(a)} . 


The diverter vane angles were 

X = sin '^ V rel /V jet " B/T net } and (III-A5-2) 

o = s i n" 1 { { U re1 / Vj et - A/T net - AC L *sin(ct)}/cos{X)} . (III-A5-3) 

The reaction control angles and the aerodynamic control angles were 
s x = = l( I 2 -I y )-qT-T„ et -{(V re ,-z 1 )/V jet -sin(A)-{z 2 +x 1 -cos(5)}}- 

Lp‘p-(l/2)*p*D^*S*b*(Cj t0 +C J j -fsH xz -p-q}/{L s +(l/2)-p- (III-A5-4) 

B x 

Tg.s-b-c^ } , 


6=5 
y e 


{(I K -l z )-pT+I xz -(p 2 -r Z ) + 2.I e -B.r- T Bet .{tt re i-* 1 -U r6l .z 1 )/ 

V J . e t+c o s(A)*{22*sin(a)+x 2 * c0 s{a)}+ c *A c M o}-M q *q-{l/2)*p* 

U f' S * c -<W ( W C^O/aMl/y-c-C^,} ' 


( III-A5-5) 
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{% +(V2)-P*uJ-S.c-C M 5 } , and 

y e 

6 z * 5 r - f(V I x ),p ‘ q+I x Z -‘'- r - 2-I e- a - , < + W (V r er x 1 )/v jet- sina) ’ 

{x,-H,-Sln(8)>}-N I ,T-(l/2)‘P - U^‘S-b-(C no +C r|S •6 a )} / (Nj + ( II1-A5-6) 

a 2 

(V2)-p-LT|-S-b-C n<s } . 

The control parameters had certain bounds or saturation points. 

They were 

T mav - 5515 pounds , 
gross max r 

o = 0° to 70° , 

X = -25° to +25° when used, otherwise 0° , 

V 6 y’ V 5 a = ' 20 ° t0 +20 ° » 

S 0 = -25° to +15° , and 

S r = -30° to +30° . 
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C This PROGRAM SIMULATES a yTOL AIRCRAFT <IN HOVER) which is 

-C SUBJEC-TEO— TO A -gust- F-tELO*'- THE PROGRAM- READS* THE' INITIAC " 

C WIND CONDI TIONS-'U Si f HEADWIND MAGNITUDE " GSMA&CU* V* OR W)* 

-C THE MAXIMUM MAGNITUDE OF THE COMPONENT GUST - OENQ M * THE PERIOD 

C OF THE GUST {DFNOM*iO--PERlOD3lO SECONDS, DEN0 Hs25«-PERI0D=5 

G -SECONDS* OENOM*:1^" , *"PERIODs1 SECONDS ) - AND JGUST* AN IDENTIFIER 

C TO SPEC tfy which cumpunent gust is being CONSIDERED <1* LONGI- 

t rtJDlNAL""GUST ?> RtJSl TIVE 'LATERAL' GUST" 3* NEf.ATl VT'CATERAL"' GUST 

C 4, POSITIVE VFRriCAL GUST 5* NEGATIVE VERTICAL GUST), THE 

(T PROUPA A THEN T‘AkE5' ThCSE PARAMETERS AND ITERATES UNTIL - THE '' *' 

C MAXIMUM MAGNITUDE Of _fHE_GUST THE AIRCRAFT CAN WITHSTAND IS FOUND, 


_ REAL MASS* Ix*IY* lZ*IXZ*LONE*IE.NG*LOXlX*LPl X* HOYT Y* MQl Y# NOZ! Z* NR I /* 

1 NONO I M* I.AMPDA, vl I C 1 0# 1) * YY ( 11 * U) 1 > * T J mEU 0 1 > * LS IAR< 1 0 D , NS f A R( 1 0 1 ) 
2*P[)URC iuU.Unorf I01)*R0MTC ti)l)*YA(6* 101 )*Y8(8* l 01 ) * THR STGC 1 0 1 > 

■ “COMMON/ESOC UM/N * XI * Y I ( y? ) * X* Y ( 99 ) * TH * H* KSTP* 0 Y C 99 ) 

C0MM0N/FUNC1 /G* THRUST c 10 1 ). MASS* SIGMA (1 01 ) * L A MBPA C t 0 1 ) * 

TUO* w"AREfl* CLPn*DELCL*ALPHA/CO* VREl.*UOHAR*CY* IX* IZ/IXZ* I Y* 
2 ZTWU*LUNE* 70 NE*LPIX»LUXIX* 0 E|.XC 101 )* SPAN* WREL* CLO* CL DEL A* 

30ELA* CL1F.TA* IEN G ■ OMEGA * X TWn* X ONE* NR IZ * NDZ I Z» DEL? C 1 0 t )* CND* 

4 CNDELR*')ElR,CNDELA,MDYIY*DELY( 1 U 1 >*NGIY*CHORD*CMUPO* 

5 DELCM 0 * CMALPO* DELCHA* CHD* CMOELE* DELE* BETA* J, UREL* V JET* RHO ' ‘ 

6*CPA|.PH,CL4LPH 

DATA G* HASS* V JET* RHD* WARFA * SPAN* CHORD* I X* I Y* 1 1* 1 X 2 * ZONE* IEMg*XONE* 

1 XTWO*ZTwO* L ONE* LDX I X* LPi X* KDYI Y» MOI Y*DEL CMA*CMOELE* NDZI Z*NRlZ*t;My* 

2 ANGVFL»LA-nDA*C 0 ALPH*CLALPHX 32 ,lZ 4 # 130 . 3537 * 1 6 1 3 . » 0 . 002 3 Z 8 » J 82 , 6 9 

333 . 83 ' 5 . 56»2 34 (■#* 34 Ou **54 00 *> 180 «' 0*58 3 # 6 . 0 * '‘ 0 * 0133 * 0 . 916 * 0 . 16 /* 

4 0 . o 63 |*'o* 45 *o« ft i 5 o*“ 0 » 15 * 0 , 576 # «o ,OlZS* 0 *OJ 66 * w O*? 0 * " 11 . 4 * 1 720 ,* 

51 01 * 0 , 0 * 0 . 0 * 0 . 0 / 

DATA -Vl*WI*Pl*OI*R!*THETAI* PH II* PS 11/8*0,0/- 

loo read c 5 * io*eno« 9 ouo) UW 1 * GSMaGU, gSMAC.V* GSM ftGW* DENOM* JGUS t 

10 FOPPATCSF 10 . 2 * J 5 > - - -- - 

Ul s AFJS(iJ.*I)/ 3 t 0 

ITB »0 

60 H=0 • 1 0 

I T E — 1 

TIMEC 1 )a 0,0 ----- - 

_ TlM£X=nMEC n 

J *1 

L *2 

V = VI ______ 

r=p i 
RbRI 

— THE T A = ThETA I - - 

PHlsPHlI ' 

P5-f*PS H 

UGsO.O __ _ 



olooo j I O oo wj I I i '■ l 1 ' « «J jfJ 


WG^O.O 
GO TO 20 
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t n t r n ou c r iw 'o r r$r^$sr-rtiu ci 7 r f i e l or “try snp~ is wired/' 

RLPLACt TtlE APPROPRIATE ELATION 8Y# FOR EXAMPLE# QKeGSMAGW. 

30 PIE»3#I<M59 

JIKaJ-I ' “ 

GU = 0«5*<lSMAGU*( i 40-C0ST PIE* J I K/OENOM>> 

GV*0.5*OSMAGV*ti , O-COS CPlC* JIK/DENOH J 3 ' 

GW=0,5*GSUAGW*(I ,0-COS£Plt*JIK/OENOM3J 


ROTATION OF The gust VfLGCITy flELO FROM THE ATMOSPHERE-FI X£0 
REFERENCE FRAME INTO THE AIRCRAFT REFERENCE FRAmE, 


UGs&U+COSCPS r*) + COSCTHETA J + GV*$IN(1>$I )*C0$<TH£TA)-(iW+SlN(TH£TA} ' 

VG s GV*tcns(PSI >*CUS(pHl)+S INC PS D + SIMC THETA )*SlNT PHI >W,U*( Sl« 

■ ’ 1 T p SJ >*COS<PHl>-cOs(PsI J*SlN<THETA)*SlNCPHl)) + G^ t COStTMETA)*slN 
2(PHH 

WG*GW*CuS{PHl ) + CUS(THETA} + flU*(COS<P4h*SlNcT«ETA)*CtJSCPHl3*S!NfpSI 
1 )*SiNf PH I ) >-GV* ( COS(pSl>*SlN(pHl )-$!»< PSl )*S I M< THETA ) * COST PH U ) 

- I UW»lJ 1 /I*COS(PSI)*COS(THErA)- - — - - -- 

VWsUWl*(CU$(P5l J+SINCIUETA)*SINCPHI J-SfNCPSl)*COS(PHI)j 
HW-OWI*{SIN(ThET A)*COS fPHJ )*C0SIPSI >+Sl«CPSl )*S1N<PHI }) 
UREL 3 U-UW«IJG 

VHFL«=V-\/W-VG -- - • - - — 

WREL^H-NW-WO 

' ■"'■oa = bQRT<IJRFL**24W«EL'**2X J 

U0RARsSciRr(URt!,**2+ VREL**2 + WREL**2J 

ALPHA = ARSTNCV)RFL/00) 

BETAaAHslNCVREL/UORARJ 

CAlcULATIUN^OF^ t he aerodynam ic COEFFICIENTS _ 

IffUO.LT.Sl*) CLPO=0,6l +0, 002764 *UO 
IFOJ'UGE.M . ) ClP0=0.75 

0f:LCL s 2,4O3009F*'05* jn* + 2-2.?48398E*07*O0**3-l . 3321 49 £-03*00 
CO*u.li+l./(?»l52H»i*uO+l2.5> 

IFC aBSOIE TAKLT.O.l B9) CY = - 1 . B5l *8E I A 

rFCARSC>TETA) t GF.O,l 89, AND iABSt OETAJ.LT; 0.4U)CY=C *0,35 + 0; 393* 

1 cars ( beta )-0 , 1 49 n*lA3S(OLT4)/tiETA) 
IFfARSfyErA).GF.0,Aifl,AND«ARSCRErAJ.LT*t»22) CY-OO. 26-0. 362* 
1(AB5CBLTA)~0, AiG))*CAOsTRETA)/3ErA) 

TKAijS(tjE| A), (IE, 1,22) CYs-G,55*(AB5< BETA) /SETA) - - 

IF(AMS(3ETA).LT.y,l69) CtUs-0 .0B3 + BF.TA 

I[ ( msft'LTAl.GF.O* )69*ANP*AHSf 8ETA3*LT»0*4ri ) CLO* ( "0 . 0 i A + 0 , 0066* 

IC AH', (BETA > m 0 9 i 69) ) + CA3sCBETA)23ErA) 

IFC ARStsLTAj.GL ,0,4/1 ,AN0 ( ARSC3ETA),LT,!,395) CLO*< -0,0 1 2-0 ,00665* 
) (AIJSC BETA 3-0.471 ))*CABS(HETA J/BETA ) 

If(Ad$<tfLTA>»GF»l*3*5) CLO=-0 . 02* UtlSC BFT A )/ BETA ) 

CI.dFLA- 0 ,00J*CaS(PETA)**2 

Cl Bf TA-0 ,00A5*CLP0-0 ,0122' — “ 

1 i too. LI. 34.) C M UP 0 *- 0 . 6+0 « 00?37*l)0 

If f Ilf), :jf> 34,0* AND. UD.LT.S1, ) CMUPU C *O.S2+0, 01 A8*(U0'3A , J 
iFtim.GL.S! ,0. AND.On.LT. 101 . ) CMDP0--0. 27+0. 001 579* 

1 con-si . ) - - 

IKiin,t,E. 101 . ) CM0PU=-0.!9 

fjj L<M0"A. V003SF.-07*00**3-3, 16A62F-0S*U0**2+1 »S29505E-‘03 + U0 

IFECM.LT » 3 4.) CMALPO 5 * A , 1 + 0 » 0 A 1 ttS + UU 

irtoa.GE. 3 4.0. AHDfUO.LT. 5l.> t'ULP0 = -2 , 7 + 0. 1 007* ( (JO-3 4 , ) 



s 

IFCU0.GE.51 « ) CMALPO»M.O 241. 

If(ABS<BETA),L T . 0.1905) CN0«0. 1 1 3*BETA ' ' 

IFCABSCgETA ) ,GE .0,l905tANU.ABS(BETA).LT.0.590) CN0».022#CABS(BETA > 
1 /BTTA) ' 

I F< AFlS<oETA) .GE. 0.594. AND. AfiS< BETA) »LT. 1.395) CNO®< 0. 02?+0. 1 3*8* 
1(ABSC0 EtA)-O, 59A))*<A8S<PETA)/8ETA) 
lF(AflS(H£TA},GE, 1.395) CNU« 0 , 1 3*C ABS (BETA ) /BET A ) 

’CNDl’l RaO.00 1 18*C03(BE f A )**? ’ ~ “ 

C WDEl.A«-OtO OO 1 4*C0 5 < 8 ET A U * 2 

c - 

C INCLUDE the LIFT AND DRAG COEFFICIENTS WITH RESpEcT TO THE 

"CT' “ AWGLE "Of ‘Attack'at THIS" POINT' ’ 

C IFCUO.LT-. 34.0) CDALPHsO.O 

~C “ IFCUO.L r/34,0) CLALPH*5.3-n,o3235*UO " " ’ * ” * 

_C_ lFCUn,f»K.34.0) COAl p H®-O.J_ 

C 7 F C Ufl . G E » 3 4 . 0 , A N U . U 0 • i. T i 5 1 . 0 ) CL ALPH «4 # 2"0 • 0 6 1 1 * £ On " 3 4 • 0 ) 

C IFCJ J 0.f.E,51.0) f,l.AlPH = 3.2 

C INCLUDE ALPHAfl lOEGKEES (OP WHATEVER ANGLE) AT THIS POINT. 

"C " ALPriA=AL p HA+0,i92 " ’ ' 

C 

lFCnHEl(."Eo^O;orA"ND.K.E'Ci.rT GO TO' 40“ 

IFC IDENT.Ert. 1 ) GO TO 50 

ynu,i >aur " - ‘ - — 

YIIC2.i)=VI 

ni{3#n=wr 

mu.i >*pi 

nrr5Tr) . B(jr — 

Yii(6/i )=ri 

mC7>M*THETAI - - - — 

y 1 1 c 8/ 1 )sphi r 

. --yil( 0 ,l )=PSI I • ■ — - -- 

Y 1 I C 1 0 * 1 ) s* l 00 » 0 

00 - 6 0 -IaWH> 

iu=j+i 

• Yj(i)*rn(f») > - - - 

Y‘YUJ/1 )sril (I, 1 ) 

- - - 60 CONTINUE - - - - 

YI(il>*0.0 

yj( 1?)« 0 ,0 

YY( 1 # 1 )sT IME ( 1 ) 

- N = 12- - ----- - - “ 

XI»0.0 

KODE-1 

c EsnoFci calculates the state vector, the output of fsopeg 

C ARE T HF VCD'S. " - - - - — 

C 

- - 50 CALI ESOOEOCKODE. 1 ) - - — 

K0DE=3 

JlaJ-1 

JLASTaO 

TIMET J)»TIME< JI)+H 

- -- -Y YC-1 . J >«T I HE( J ) — 

TIMt'X^TlMECJ) 


C STORAGE OF THE STATE VECTOR IN THE ARRAY YY 



oio o 


00 70 Ialfio 242. 

I K *» X + 1 ' ---- - - 

miK#j)Bru> 

TO-C'OH'TTNTJr ' 

U*¥U) _ ___ _____ 

W«Y(3) 

PsYTA) “ ‘ 

i 3 s Y ( 5 J 

RsYCrt 
THET Ab Y(7 ) 

-p H I = Yt8) - - • — • - - 

PSIsYCV) 

ITC*ITC + l- - 

NNsL-ITC 

-THRUST(t)“THRUST<N^7 — - - 

SJGMACL)=$IGMACNN) 

0ELX(L)s0E L X(NN> • - ... 

DELYCUaDLLYCNN) 

DELZCDbOELZCNN) 

c 

~ C -STORAGE or THE * CON TROL" VF C TOrt ' I H THE ARRAY YB ' ’ ~ 

C 

YB( UUbTJMEcL) - - - 

YBC2#l.)s*r,tRUSTCL) 

YBf 3 # L)=T H RUS T CL)*1 .23986 .... 

YBC4#L JaSIGHACL) 

YBf6;i'^sOELX(L ) ----- 

YB(7>L)a()ELY(L) 

YB( 6* L )sOELZ( L ) 

CALCULATION Of" ANGULAR ACCELERATIONS (EyUl V a LENT Tn THE 
CONTROL EFFORTS), 

c - - - -- 

LSTARCL )a( lY-rZ)*Q*R/IX + IX7*P*Q/IX + rHRUST(L)/IX*(Z0NE*/REL/YJFT 

I»STN<LAiiROA(L))*(ZTmO+LONE + COS(SIGMA(L))))+LPl X*P+Lf)Xl X*OELX(L ) + 

20,5/IX*RH0*unBAR**2*WAKEA*SPAN*(CLO+CLDELA*DELA+CLRETA*RETA) 

NSTAR<L) = C TX-lY)*P*G/lZ-lXZ*a*rt/lZ + 2,*lENG*DMEGA*G/lZ 

HTHRltST (I. } /I Z*{ SIN (LAMBDA Cl ) )*( XT WO -LONE *S I N ( S I GM A ( l ) ) ) -X0*'t 

2*VRE-L/V.j£T) + NRlZ*R + NnZl Z*UFlZ(U + 0.S/IZ*RH0*U0BAR* + 2*WArEA 

3*$PAN*CcNO+cNDFlR*OELrt+CN0FlA*OtLA) 

- Pdot(l>*ix*i7/ci>*iz-ixz**?)*clstap(lJ + i x 2*nstar(l>/ix) 
ROOT(L)=tX*lZ/(IX*lZ-lXZ**?)*ClXZ*LSTAR(L)/lZ+NSTAH(L)) 

- OOOT(L>={ I?-IX)*P*R/I Y+IXZ*<R**?“P**2>/I Y-2.* IENG*0MEGA*R/IY 

UTHRIJST(L) / I Y*( < ^REL*XONr-UREL*ZONF.)/vjET+COSCLAMBnA(L)) 

2 * ( Z-? * S I N ( $ l G M a ( l ) ) + XTWO*COS(S|GMA(L >) ) +0EL CMO *CHORD )+HO Y JY * 

3DELY(L )+MQl Y*G + 0,?>/T Y*.RHn*UG**2*wAR£A*CH0R0* (CMOPO + CCMALPO + 

. 40f LC mA)*ALPHA+0,S*CH0R0*CHR*Q/U0+CMDELE*DELE) - 

YAC 1/L )aTIMt(L) 

C STORAGE, of the angular acceleration^ in THE array YA, yA (5 ) 

— e and~yat«t"store -theta sr twottate variables# longitudinal 

C AND LATERAL DISPLACEMENT RESPECT I VF-i. ?« 

Y A ( ?* I ) sPDOT ( L ) 

YAC 3«L JsQOnTCl * " _ - - _ . 

YAC4#L)sROOT C) 

YAC5#L)=Y(U) — “ 

YA(6,L)oY f W) 

PP0TA*P[)0 



o o o o ( ioo 00 , 00 , ;ooo | ( | o al 


OOOTA*QOUTCL) 
ROOT AbRDOT ( L ) 
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iFCApscTFiETA),GT,o.4363A) ttf *2 

J^F(AfiS(PHI),GT.0.4 36 34)_ ITF52 
iFTAHS'fPOOTAJ.GTt 1.361) ITF = 2 
IFCABSf9D0TA),GT, 0,698) ITF = 2 
IFC ArS{ ROOTA) » G T * 0 # 332 j ITF=2' 

_ L*L + 1 _ _ 

IFC iTF.E«.?.OR.OtGE. iOf ) GO TO 1000 
IDENTsl 

KaK + 1 " 

TTc“ fs T he ”'c o' u nt e rT o“ $ i m GTa t eT h e time delay"’of the pilot” 
ifcitc.lt,3) go ro 30 

K«K»1 

calculation nr the control vector, the first time through the 
PROGRAM# "T H IS VECTOR is CALCULATED IMMEDIATELY AFTER THE 
AERODYNAMIC coefficients ARE CALCULATED. after that# This VECTOR 
“'IS CALCULATEU WHENEVER I TC = 3 , 

-4(T-A=MASS'*( V*P-H*O“6*STN(THErA)3+0.5*RH0*U0**2*WAREA*( < CLPO + CLALPN* 
l AL P HA )*$I N( ALPHA )-( CD* COALP H* ALPHA)* COSC ALPHA) ) 

7BsMASS*(W*P*U + R + ff*‘STN‘(pHn*C0S'( THETA ) ) +0 , 5 *RHO *U0R‘AR *#2 
1*WAREA*C Y 

'C s MASS*(U*9-V*P + G*CaS(PHl)*COS(THETA))-O.5*RHO*U0**'2 ' " ‘ " 

l*WARFA*((CLPO+CLALPH*ALPHA)*COS(ALi , HA) + (CO + CDALPH*ALPHA)*SlN 

2( AtPfiA ) ) ---- - ~~ - - ' 

AAsliFLCL 

O*T-.-U0hA«**2/vUET*-* 2“AA* + ?-»2.*AA*( SlN( ALPHA )*UREL*"COS 

1 (ALPHA)*WREL)/VJET 

E=CA*UREL+8*VREL+C*WREL)/VJET-AA*(A*SIN( alpha)-c 

1 *C0S( ALPHA ) ) 

F»A**2 + B**2 + C»*2 - -- 

THRUST (K)a(-E+SfcRTC£* *2+0+ F))/Q 

THRSTGCKJsTHRUSKK )*1. 23906 

THRS7= THRSTG(K) 

CALL ENgRPM( THRST#OmEgPC>UHEGA ) - 

— ' — IF- LAMBDA- is- TO BE CALCULATED* insfkt - --- -- - — 

LAMBDA tK) = ARSIN( VREL/VJET-R/THRUSTC K) ) 


SIGHA(K) = ARSIN( < UREl/ V. jET- A/T HRUST C K ) * A A *S I N ( ALPHA ) ) 

— 1-/C0S<LAMBDA(K) ) ) ■' ’ - ----- . -- 

0ELX(K)*((IZ-IY)*Q*R«THRUST(K)*(VRT L*20N,:/VJET-SIN(I.AMBDA(K)) 

1 + <ZTW0 + l0NE*C0S( SIgmA (K ) ) ) )"l PI X* P* I X-0,-)<RH0 + U0BAR + *2 + wAREA ~ 

2 + SPAN*(cLr)+CLBETA* 1 3LTA)“IxZ*P*Q)/(Ll;Xly^IX + 0, , ;»*RHO*U08AR**2* 

3WAPEA + SPAN + CL0ELA) — - 

DELY(K)s((lX-I7)*P*R + (p*-*2-R**2)*lX/+2 , * I E NG * OMEGA +R -THRUST ( K ) 

_ _1^<-CWRFL-X0NE-UREL*Z0NE )/VvIET + COS(LAhHDA(.<) )*CZTWO*SIN(SIGMA(K) ) 

2 + XTwd*CnSCSlGMA<K) ) )+CH0RD + 0ELCM0)"M0I Y*N* T Y-0.5*RH0*U0**2*WAREA 

3*CHORD*(CMOPO + (CHALpn + DtLCMA)*ALPHA + 0,5*CHORO*CMQ*9/UO)) / - - 

A{M0YIY*TY + 0,5*RHU*UU + *2*WAREA*CH0fU;<-CMDF IF ) 

-DELZTKJaT ( I Y-IXJ*P*Q+lXZ*0*R-2.*IFN«i*0MFuA*(3-THRUST(K) + <-VREL + X0N£ 

l/VJET + SINCLAMBDACK) )*(XTW0-L0NF*Sl.'i(SlG*5A(K) ) ) ) -NR I Z *R* I Z-0 , *>*RHO 
2*UOBAR**2*WAREA*SPAN*CCNO + CNI)CLA*OF LA> ) /( Nt)/I/*IZ+0,5*RH0 
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3*U0BAR**2*WAR£A*SPAi-J*CNUE LR) 

-Y B (I,K)s riM E(K) ~ ' ‘ ‘ 

YB(2»K)aTHRUST ( K ) 

YBO,K)arTriR5TGCK) 

YB(4,K)*SIGMA(K ) 

YB ( 6, K ) sDELX ( K ) 

YB(7,K)*0ELY(K) 

YB< 8,K HOELZTK ) - - ‘ - ' 

C 

C CffF^»'rPtJP-*SATURA-T-I-ON-Or'COHT L ROt--pARAMrT£RS. 

c 

IFCI)FIX(K).GT,20.0) DELX( KH2Q.0 - - ■ 

IF(OELX(K).LT,»20.0) DEL XCK H-20 , 0 

IFCOE.LV<K>,GT,?0.0> DELYCKH20.0 

IF COEIY(K) .LT.-20.0) DEL Y ( K ) s-20 . 0 

IF(DE12(K}.GT.20»0) OELZ ( K )®20 « 0 - — — 

IF([)FIZ(K),LT»*"2Q«Q> OELZ(K) = -20,0 
DtLAsOELX(K) ------ 

DLLE=DLLY(K) 

— — - -OKI K*rElZ<K) - --- • 

IFCOFLE.Gr.15,0) DELE-15.0 

IF<S1GMA(K). GT.l. 22173) -SI6MAfK-)sl.-22173 " ‘ 

IFtSlGMA(K).LT.O.O) SiGMA(K) = 0.0 

-JFfTHRSTGCK ).GT,5!>15. ) THRUSH K ) = 5 b 15 . /1 .23996 

I T C 5 0 

K=*K+J - - 

GO TO 30 

C ITERATION PROCEDURE plNOS THE RANGE A TO A + 20 FIRST, THEN IT 

C FINDS The RANGE a TO A + 2, AND FINALLY THE EXACT VALUF < 

IT CAN r,Q TU EITHER GREATER OR LESSER VALUES FROM ThC INITIAL 
•VALDES OF THE GUST MAGNITUDES. 


rooo-iFCTTF'.r f i.n 'GnrT0'99a 
IFUTB.EQ.l) GO TO 996 
IFCJgDST.EO.I) GSMAGU“GSMAGli4 10.0 
IF(JGUST,E0,2,0R, JOUST ,EQ.3) GSM AG V*GSMA G V+ 1 0 . 0 

“ ' IFCJGUST.EQ.A.OR. JGUSl.E0*5) GSMAGWaGSMAGW+10.0 

I T A* 5 

GtJ~T'(T 80 

996 IF< JGUST.E9. 1 > GSMAGU'GSMAGU+l . 0 

- -~IF( Jr,uST.E0.2,0R.'JGUST.EQ.3) G$MAG9=GSMAGV+1 , 0 

IFIJGUST .Eq,4,0R.JGUSr,E0.5) GSMAGWaGSMAGw+t .0 

- -ITD«l - - - ” 

GO TO 80 

9.9-6— I F<ITA*E9.l) GO TO- 997 

IF < JGUST.EOl ) GSMAGUsGSMAfaU-10.0 

If(JGUST.E0.2.QR. JGIJST.ES. 3) G$mA(»V = GSMAGV-10.0 

IF C JGUS r . EO , A.OR.JGUST »EO,5) G$MAGrfsG$MAGW-i 0 ,0 

—GO^TO 80 - 

997 I F C JGUSJ.EO. I ) GSMAGU = GSMAGU**2.0 

IF-CJGU5T.E9.2-,0R.JGUS F.tQ. 3) GSMAGV»GSMA&V-2 « 0 — • 

IF(JGUST,E3, A, OR, JGUSI.E0.5) GSMAGWaGSMAGW-2,0 

|-T8*1 - - - - - - ' - 

ZFdTD.EO.i ) GO TO 1001 

- GO TD 80 - 

1001 WRITE(6«800)(( YY( I, J), Hl> 1 1 JtAST) 

800 F0RMATC1P1U 11,4) - ------ ...... 

WRI TE(6,801) C C YAC I, J), Hi, 6), HI, JLAST) 

301 FORMATC 1P6E1 1 .4) 
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WHITE (6,602)C(YB< I, J)/ I *1 / 8 ) / JM , JlAST ) 

HO? FORMAT ( 1 POEM .4 } — - — 

go ro loo 

- - 9000 S f OF — 

£N0 

FUNCTION F ( I » X jf Y ) 

c 

C FUNCTION SUBROUTINE 'FUR' ESOOEUT” " ’ 

C 

• COHitON/ESO/NOBO 

REAL Y(99) 

COM.' iON/fUNCl /G> THRUST < 101") MASS# $IC,MA ( I 01 )> LAMBDA ( 
tUU» viAREA'CLPO/DEICL/ A lPHA/CO/ VRE l*UGBAR/C Y / IX, I// IXZ/ iy, 

2ZTW)> LONt, ZUNp, t. F I X, LOX IXj> 0£LX{ 1 0 1 ) / SP AN, WRFL , CtO# CLOEt A-> 

3DEL A, CLnET A, I FNG* OME GA , XT WO/ XONE , NR I l > ND Z 1 1’ D£ LZ ( 1 0 1 ) /_CNO>__ 

- - 4 C NDF.lR>oELR»CN[)ElA,MDYlY>0iriYC10n/M0lY/CH0RD/CM0P0/“ ~ 

5DELCM0»c M ALF0»DELCMA/CH^#c M 0ELe/0ELE/8ErA/ J/UREL/VJET/RHO 
6,C0ALPH,CLALPH ----- - 

IFCI.EO. 1 ) FsY(2)*Y(6)-Y(3)*Y(5)-G*SIN(Y(T) ) fTHRUST C J )/MA$$ 

1 *( SlNCblGMAC J))*C0SCLAMi3DAC J))*UREL/VJ£T + DELCL*SINf ALPHA ) )+0.5/ 
2MAss*R f <n*un**2*WARFA*c cclpo+clalph*alpha)*sin( alpha)-(co+cdalPh* 

3 ALPHA )*'C n $ ( ALPHTFD' 

IFtl.FQ.D RETURN 

IFC1.E0.2) F=YC3)*YCA)-YC1 )*YC6)+G*SIN(Y(8) )*C0StY(7>) 

l + thkust c J)/mass*i sin( lamBdai J))-vrel/vjet)+o,5/mass + rho*uobar 

2** 2*V(AREA*CY ----- - 

IFCI.tu.2) RF TURN 

ifci .eo, o--F*Yri^*Yr5) , -Yr2)’» Y T4)-+G*cosfYC8) )*cos(YCtn ‘ 

t - THKIJS r c U)/MASS*(C03(5IGMA( J))*COSClAMBdA( J)) + WREL/YJET + DELCL* 

- - ?cOS( ALPHA ) )-0,5 /mA$S*RH0*U<H*2*wAREA*( (CLPO+CLALPH*4LPHA)*COS 

i{ ALPHA ) + ( COt CO ALPH* ALPHA )*SIN( ALPHA) ) 

IFCI.EQ.3) RETURN - * 

I F C ! » £ N , A ) F»(iX*iZ/{jX*lZ-lXZ**2))*aiY-iZ)*Y(5)*Y(6)/lX 
1 tlX/*YC A )*Y(5)/IX+THRUST( J)/IX*(Z0NE*VRFL/VJET-SIN(LAH8DA(U) )* 

2(/r v .0 + L0NE*C0S(SlGM*( J) ) ) ) + L PI X* Y £ 4 ) +LOX I X *0ELX ( J) + 0,b/IX 

3*RHO*UU*AR**2*wAR£A*SPAH*(CLO+CLDELA*dElA+CL0ETA*b£TA J 
A + IXZ/IX*C( I X - I Y )*YC4)*YC5)/IZ-IX7*Y(5)/Y(6)/|Z + 2,*IFNG 
5*UMCGA*/<5)/IZ + THRtJiT( J)/I7*(SIN(LAMBUA( J) )*( XTW0-L0NE*SIN -- 
6(SlGMA(j) ) )-XONF*VREL/VJET)+NRIZ*Yc6)4NOZIZ*[)ELZCJ)+0.5/1Z*RHO 

- - 7*UO8AR**2*wAREA*SPA,Ni*<CHP+CN0ELR*DELR+CNDELA*DELA> >) 

IF(I.EQ.*> RETURN 

IFCI.FQ.a) Fs(IZ-IXj*Y(«)*Y{6)/IY + IXZ*(YC6)**2-Y(fl)**2) -- ~ - 

1/1 Y’2.* IENG*0MF.GA*Yc 6)/IY+THRUST( J)/I Y*f (WREL*XONE-UREL 

- - - 2*Z0Mf )/vJET+COSClAHBOA( J ) ) * f Z T WO* S I N£ S I GM A £ J) )4XTwO + COS 

3(SiGMaC J) ) )+CHORD*i)fLCHO) + MD Y lY*DELYt J)+M9lY*YC5)40,5/IY*RH0*U0/*2 

A*wAREA*cHORD + (C MUpO + <CMAL POt-OELCMA) * ALPHA fO • 5*CM0R0*CM9* Y (5) /UO 

5 + CM[)F.LE>DELE> 

IFC I .E0.5) RETURN 

IF (I ,E0,4> F = ( lX*IZ)/( IX*1Z-IXZ**2)*(IXZ/IZ*C( IY-IZ)*Y(5) 

1 *Y(6)/Ix+1XZ*Y(A)*y< 5)/IX+THRUST( J)/IX+(Z0NE*VREL/VJET-SIN(LAMR0A( 
2J) )*cz rwO + LONE*CUS<SlGHA( J)))) + LPIX*Y( A) + L0XlX*0ELX( J) + 0.5/lX* 

3RH0*UGB A** **2* W ARE A *^»P AN *<-CL U+CLO£L A *OEL A + CLBET A* BE TA)) + ( I X • I Y )*- 

4Y£ 4 )*Y£b >/IZ-IXZ*Y{5 )*YM J/IZ + 2. * IFNG*(1MEG A * Y ( 5 ) /I Z + T HRUST ( J)/IZ* 
5 (SIN(EAmPi)A( J)>*< XThQ-lONE*SIN(SIGMAC J)))-X0NE*VREL/VJ£T) + NRIZ* 
6Y(6)+NUZIZ*UELZ( J) + U,5/lZ*RH0*U08AR**2*WAREA*SPAN*(CN0'tCNDELR* 

-- -- rOELrt+CN()ELA*OELA>) - - - — - -- " 

IFCI.E9.6) RETURN 

I-FM .EU, 7->~F-y<5)*C05{ YUT))-Y<^>*SIN(Y-(6')-) ■ 

IF £ I * EG » 7 ) RETURN 

--IFU ,E«.-8>- FsY< 4) + Y(5>*SlN{Y<6))*TAN( Y (7))*Y(6)*C0S(Y(8)) 
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1 *TAN( YC 7 ) ) 

IFCTTEr. 8 ) “RETURN 7 

IFCI.EQ.9) Fs(Y(6)*C0S(Y(8))4-Y(5)*SIN(Y(8)))/C0S(Y<7)) 

If^'tTE^T9T~RETtJHN 

IFU.E^.IO) FaYC2)*sINCY(8))*COS(Y(7)) + Y(3)*COS<YC8))*COS(Y(75)- 

mi") * s f m < r c n y 

IFU.EO.lO) return 

IF( IlEtti 11) F=YC 1 3+C0SCYC73 )*C0SfY( 9) T+YC2)*CSIN(Y(0))+SIN(Y(7))+ 

lCOSCY(V) )-C0S(YC8))*SIN(Y(9)))+Y(3)*(C0SCY(3))*SIN(YC7 ))*c0S(Y(9)) 

2 -rsi Nm?7) *"$ r n crrsrn 

1FU.EQ. in return 

IFU.Ee. 12) F = Y(1')*C0$(Y(7)')*SIN(Y(9)) + YC2)*(SIN(Y(8))*$INCY{7))* 

lSlN<Y(9)) + C0S(Y(8))*CCsCY{O))) + YU)*CC0S<Y(8))*SlN(YC7))*SlfJCY(9)) 

2- S INtY(8)>* C 0SCY(9))) ------ 

RETURN 

END — — 

SUBROUTINE EMGRPM<ThRST#OmEGPC#OmEGA) 


SUBROUTINE TO CALCULATE ThE ENGINE ANGULAR 
VELOCITY FROM THE- THRUST* — — 

- - A NGVE-L-at72fi-i 

IF<THRS7‘LE.4370. ) GO TO 300 
GO - T q -3f)i 

300 OMELPC*C THRST-3 305. )/88. 75 + 80* 

- - - — GO TO 302 - 

301 BA«A.52446?3£-09 

8BS9T3T30978E-06 

BC«1 .24984971-02 
' " BU*T MRS T “4 3 70 , 

OMFGPC»nA*BO+*3-BB*BD+*?+8C*ftO+92.0 

302 OMEGA*0MLGPC*ANGYEL/10O. ... — -- ... 

RETURN 

“END" “ ~ 

SUBROUTINE £ S 0 F> 6. Q ( KilDt * NRO ) 
common/esu/hobo 

COHHON/L$OCOM/N#XI# YI(99)#X#y( 99), TR# H# KSTP# DY< 99 ) 

REAL XC, YCC99)#0YPC4#99)#S(99> ' ' ~ 

DATA H0/0/,0YP/396+0/ 

DArA“MS T r>~NSTP>~iin » rtf?* tn3# in4/“o #~ o #-“r>- ?.r~ 3 # —4 t 

TR * C-l.O) 

“ NUB OS r - 

IF(N.GT,99) GO TO 7777 

- GO "TO" "( 1 000# 2000# 3003)# KODF 

1000 DO 1001 I=1#N 

YU)- -s--Y I U ) 

1 00 1 DY(I> S F(I,XI#YI) 

INI = 1 

- IN2--a--2 - 

I N 3 * 3 

“”i N4~s~4 

1050 MS TP = *STP 

- NS TP -8 0 

HO c H 

“ GO TU“4000 

2000 GO TO 1030 

” 3000~ I F "( H0T\|E7 HTGdTT{rT05Tr 

GO TO 4000 

300T RETURN- 
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4000 00 4001 la 1#N 

4001 OYPC INl , I ) a’DYCI) 

I F__{ NS fP.LE .2) GO TO 4500 

DU 4002 1 = t # N 

4002 Yc (X) =» Y ( i ) + H*f55.0*DYp( InIp P"59,0*0Y P C In2^ J)+37»0*0 Y P< IN3# I> 

1 ' * V « 0 *D YP ( 1 N 4 "# I ) ) / 2 4 » 0 " “ ' 

HSTP * M5TP-1 

IF' ( MSTp .UE'.Y) >“5iT7QTI'0ff 

DU 4003 1*1# N 

4 o 03 y cit** - yc m 

X = X ♦ H 
" GO TO 4 300 

4100 X = X + H 

DO 4101 ~I a 1VN “ ■' ‘ 

4101 DYCI) = FU# X# YC) 

DU 4102 ' 1*1 #N ‘ ” ’ 

4102 Y ( I ) = Ytl) + H*(9.0 + 0Y< I ) + 19.0*DYP( INI# I )"5 , 0*DYP( IN2 j I ) 

. p - -- + DY p(IN3#n)m*0 ‘ ... - •- - - 

MSTP * KSTP 

- T *~ s - o;o — " 

DO 4103 I* I # N 

5 1 0T3“Wvnr iT~+“TT 1 ?'30*rY{ I >^TcTDT5**2 
TR s SWRT(TR) 

go m 48T>0" " ■" 

4 500 DO 4501 I*i#N 

SIT) 'i 'H*OY(I) “ 

4501 YCU) * Y ( I ) ♦ sCD/2,0 
X“T"~H77TtJ 

DO 4502 1*1, N 

“DYtn = fci# xc# ^rtn 

SU) * Stn f 2.0*H*0YU) 

4502 yc<i> ■ yu) t (H*fmi)>/2.a — - - - - - - 

DO 4503 I * 1 # N 

DY( I) •* F( I # -XC#-Ye-) 

St I ) * SCI) + 2.0*H*0YCI) 

4503 YOU) * YU) ♦ H*DY<I> -- 

XC = X + H 

00 4504 I a ! , N — - - 

DYCI) =» FCI# XC# YC) 

-+ -h*ovu-> 

4504 Y ( I ) * Y ( I ) ♦ SCI>/6*0 
NSTP = NSTP + 1 

4 8 00 DO 4801 I*1#N“ 

4801 DYCI) * FU# X# Y) 

1 a - jt) 4 

I N 4 s 1 n3 

IN2 = INI 

GO TO 3001 

7 777-MR I Tg~U # - 77 7 ) 

CALL EXIT 

777 FORMAT { 1 ril# 9H-T00-HUC#) - — - — 

END 
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APPENDIX 7 

Inclusion of C La and In the Aircraft Model. 

C L and Cp are to be Included In equations (A4-1) and (A4-3) in 
Appendix 4. The altered equations are 

u = ••• + 0/2) ♦ p • Uj; • S • {(C^ + C La • a) • sln(a) - 

( c p + C Da ♦ a) • cos (a)} • (1 /m) (III-A7-1) 

and 

w = ••• - (1/2) * p * U p * S • {(C Lpo + C La • a) * cos (a) + 

( C 0 + C Da • a) • sin(a)} «(l/m) . ( III-A7-2) 

These two coefficients are also included in the terms A and C 
of the thrust equation (A5-1) of Appendix 5. These terms are altered 
in the same manner as above with the exception that (1/m) does not appear. 

A = *** + (1/2) * p • Ug * S * {(C Lpo + C La * a) • sin(a) - 

(C D + C Da • a) * cos (a)} . ( III-A7-3) 

C = ••• - ( 1 / 2 ) * P • • S * {C Lpo + C La * a) • cos (a) + 

(Cp + Cp a * a) • sin(a)} . 


( III-A7-4) 
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APPENDIX 8 


Determination of Gust RMS Intensity. 

The longitudinal gust will be used in this example. The value of 
from WINDY was 52 fps. The gust field was moving at a speed equal to 
the headwind, 10 fps and the gust frequency was 0.10 cps. The displace- 
ment of the aircraft was 23.9 feet in the positive x-di recti on. The dis- 
placement was in the opposite direction to the movement of the gust field. 
Therefore , 

2d m = (10) • (1/0.10) + 23.9 = 100 + 23.9 = 123.9 feet. 

Also L u = 854 feet which implies that 

d m /L = (123.9)/(2 • 854) = 0.0724 . 

Using the probabilities of occurrence of 10, 1, and 0.10 percent of exceed- 
ing v m (52 fps) for a specific d m /L (0.0724), the ratios (v^cr) for each 
probability can be found from Figure 12. 


v g 

0.10 + m 0.62 

0.01 + 1% 0.95 

0.001 *0.1* 1.25 


The corresponding intensities are 

= 52/0.62 = 84.0 fps rms , 
o.| = 52/0.95 = 54.7 fps rms , and 
Oq 1 = 52/1.25 =41.6 fps rms . 
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